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Preface 


This  thesis  contains  elements  of  theory,  design,  and  experimentation.  You  will 
find  that  it  is  organized  in  a  rather  atypical  manner.  In  writing  the  paper,  I  tried  to 
communicate  my  thought  processes  as  they  occurred,  for  a  report  is  supposed  to  reflect 
the  work  that  was  done.  There  are  several  places  in  the  report  where  it  appears  that  1  got 
sidetracked  because  peripheral  excursions  sometimes  span  several  pages.  Please  bear 
with  these  parts  of  the  paper  because  I  feel  they  have  a  necessary  place  in  helping  show 
you  why  I  did  what  I  did. 
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give  continual  encouragement,  make  several  “dine-and-dash”  dinners,  clean  up  after 
them,  and  be  a  source  of  inspiration.  Holly  has  been  all  of  these  for  me  and  I  am  most 
grateful.  Thank  you.  Holly,  for  your  unwavering  faith  in  me. 

Dr.  Peter  J.  Torvik,  my  thesis  advisor,  has  been  extremely  helpful  to  me.  He  is 
truly  a  master  of  this  subject,  and  his  timely  advice  and  valuable  counsel  saved  me 
countless  hours  pursuing  academic  dead  ends.  I  am  grateful  for  his  time  and  for  taking  an 
interest  in  me. 

Any  student  who  uses  a  laboratory  knows  the  value  of  knowledgeable  and 
experienced  technicians.  I  am  thankful  to  Mr.  Jay  Anderson  for  providing  me  with 
training  on  how  to  use  the  modal  analysis  hardware  and  software.  Mr.  Andy  Mills  and 
Mr.  Maik  Derriso  are  appreciated  for  all  of  the  footwork  they  did  for  me,  hunting  down 
parts  and  supplies  I  would  have  never  found  on  my  own. 

Finally,  the  men  of  the  AFTT  Model  Fabrication  Shop,  Mr.  Joe  Hofele  and  Mr, 
Jack  Tiffany  in  particular,  are  some  of  the  finest  craftsmen  I  have  ever  met.  Their 
expertise  and  suggestions  resulted  in  providing  me  with  a  test  apparatus  and  damping 
device  of  exceedingly  high  quality.  Thank  you  for  giving  me  a  superior  effort. 


Mark  A.  Scharpen 


Table  of  Contents 


Preface  .  u 

Table  of  Contents . iii 

List  of  Figures .  iv 

List  of  Tables  .  vii 

List  of  Symbols .  viii 

Abstract  .  xi 

I.  Introduction .  1 

Background  .  1 

Problem  .  2 

Review  of  the  Literature .  3 

Scope  .  4 

Assumptions  .  5 

Approach  .  6 

Sequence  of  Presentation .  7 

II.  Improving  on  a  Conventional  Constrained  Layer  Damping  Treatment  .  8 

III.  Design  of  the  Test  Article .  1 1 

Characterization  of  the  Truss  .  12 

Summary  of  Equipment  Used  .  12 

Summary  of  Truss  Dynamics  .  13 

IV.  Conventional  Constrained  Layer  Treatment  -  Analysis  and  Design  .  15 

V.  Improved  Damping  Treatment  .  25 

Displacement  Equations  for  a  Three  Layer  Sandwich  Beam  in 

Direct  Shear . 29 

Effect  of  Lever  Geometry  on  Bending  Analysis .  35 

The  Differential  Equation  of  Transverse  Displacement  for  a  Three 

Layer  Sandwich  Beam .  38 

Boundary  Conditions  for  the  Sandwich  Beam .  39 

Solution  of  the  Differential  Equation  .  42 

Calculation  of  Strain  Energy .  53 

Calculation  of  the  Mode  1  Loss  Factor .  58 

VI.  Experiment  and  Results .  63 

Apparatus .  63 


III 


Experimental  Procedure .  67 

Results .  68 

Observations  of  Results .  77 

Conclusions  Based  on  the  Experimental  Results .  78 

VII.  Recommendations  for  Further  Research .  87 

VIII.  Conclusion .  90 

Appendix  A:  Characterization  of  the  Test  Article  and  Finite  Element  Model .  91 

Appendix  B:  Analysis  of  a  Five-Layer  Sandwich  Beam  Using  the  Variational 

Calculus .  92 

Appendix  C:  Comparison  of  Experimental  and  Theoretical  Modal  Loss  Factors 
(Modes  1  and  2)  for  Various  Constrained  Layer  Segment  Lengths 
on  a  6061-T6  A1  1 1"  Beam .  104 

Appendix  D;  Derivation  of  the  Differential  Equation  of  Transverse  Displacement 

of  a  Three  Layer  Sandwich  Beam .  107 

Appendix  E:  Development  of  the  Volume  Spedfic  Energy  Dissipation  per 

Cycle  D  .  112 

Appendix  F:  Determination  of  Shear  Strain  Due  To  Bending  yt .  1 14 

Appendix  G:  Determination  of  Improved  Treatment  Parameters  h,  s,  t,  and  the 

Cross-Section .  1 17 

Appendix  H:  Matlab  Function  File  LOSS .  120 

Appendix  I:  Raw  Data  -  Conventional  and  Sandwich  Beam  Damping  Treatments  122 

Bibliography .  124 

Vita . 126 


IV 


List  of  Figures 


Rgure  Page 

1.1  Sandwich  Beam .  2 

1 .2  Deformation  of  Structure .  3 

13  Viscoelastic  Sandwich  Beam  With  Combined  Shear  Force  and  Moment .  4 

3. 1  Test  Article  With  Dimensions .  1 1 

4. 1  Depiction  of  Constraining  Layer  Segment  With  Free  Body  Diagram  of 

Slice  dbc .  16 

4.2  Theoretical  Loss  Factor  vs.  Segment  Length  for  032005 ISD  112 .  21 

4.3  1  r  Beam  Experimental  Data  With  Theoretical  Curves .  22 

5. 1  Illustration  of  Lever  Tip  Movement .  26 

5.2  Two  Lever  Concept .  26 

53  Depiction  of  Unloaded  and  Loaded  Geometry .  27 

5.4  Exploded  View  of  Sandwich  Beam .  27 

5.5  Typical  Sandwich  Beam  and  a  Differential  Slice  dx .  29 

5.6  Shear  Strain  of  Viscoelastic  Layer  due  to  Axial  Loading .  34 

5.7  Eccentrically  Loaded  Beam-Colunm .  35 

5.8  Tensile  and  Compressive  Cycle  Displacements  38 

5.9  Deformed  Geometry  of  Beams  Pinned  at  Supports  (comer  right  angles 

maintained) .  47 

5. 10  Displacement  of  the  Comparison  Function  Term  as  a  Function  of  Vt 

Normalized  t^a  Sine  Curve .  52 

5. 1 1  Illustration  of  Effects  of  Neglecting  Highly  Damped  Terms .  57 

5. 12  Depiction  of  Finite  Element  Static  Load .  58 

5. 13  Damping  Trend  of  Dominant  Term  for  Varying  Thickness  t2 .  61 

6. 1  Electronic  Equipment  Setup .  64 

6.2  Close-up  of  Truss  Damping  Device .  64 


V 


65 


63  Close-up  View  of  23"  Sandwich  Beam 

6.4  Top  View  of  5.4"  Sandwich  Plates .  65 

6.5  Detail  of  Long  Lever  and  Truss  Fastening  Plate .  65 

6.6  Load  Distribution  Plate  for  Short  Lever-SB  Plate  Interface .  65 

6.7  Short  Lever  (not  to  scale) .  66 

6.8  Long  Lever  (not  to  scale) .  66 

6.9  Distribution  Plate  Rivet  Patterns  with  Dimensions .  66 

6. 10  Loss  Factors,  Undamped  T russ .  70 

6.1 1  Loss  Factors,  Improved  Treatment,  t2  =  0.005" .  72 

6. 1 2  Loss  Factors,  Improved  T reatment,  t2  =  0.0 15" .  73 

6.13  Loss  Factors,  Improved  T reatment,  /2  =  0.030" .  75 

6. 14  Loss  Factors,  Conventional  Treatment,  t2  =  0.005" .  76 

6. 15  Tip  Displacement  vs.  Lever  Rod  Axis  Vertical  Displacement .  79 

6. 16  Force-Displacement  Relationship  -  Sandwich  Beam  \  J2  -  0.005" .  81 

6. 17  Force-Displacement  Relationship  -  Sandwich  Beam  2,  t2  =  0.005" .  81 

6.18  Force-Displacement  Relationship  -  Sandwich  Beam  1,  h  =  0.030” .  82 

6. 19  Force-Displacement  Relationship  -  Sandwich  Beam  2,  t2  =  0.030" .  82 

B.  1  Differential  Slice  of  a  Layered  Beam .  92 

C.  1  Undamped  Experimental  Loss  Factors  for  11"  Beam,  Modes  1  and  2 .  104 

D.  1  Differential  Slice  d!x  of  a  Three-Layer  Sandwich  Beam .  107 


VI 


List  of  Tables 


Table  Page 

3.1  Experimental  Modes  and  Frequencies .  14 

3.2  Finite  Bement  Model  Modes  and  Frequencies .  14 

4. 1  Example  of  Relationship  Between  Li  and  /2  for  Nominal  Parameters  .  23 

5.1  Theoretical  Damping  Values  for  Improved  Treatment  .  62 

6.1  Data  Set  1  -  Damping  Tests,  Undamped  Truss .  70 

6.2  Data  Set  2  -  Damping  Tests,  Improved  Treatment,  t2  =  0.005"  .  71 

63  Data  Set  3  -  Damping  Tests,  Improved  Treatment,  t2  =  0.015"  .  71 

6.4  Data  Set  4  -  Damping  Tests,  Improved  Treatment,  t2  =  0.030"  .  74 

63  Data  Set  5  -  Damping  Tests,  Conventional  Treatment,  /2  =  0.005"  .  74 

6.6  Sandwich  Beam  Force/Displacement  Experimental  Data .  80 

6.7  Modified  Finite  Element  Model  Results .  83 

6.8  Frequencies,  Loss  Factors,  and  Effective  Moments  of  Inertia  of  Test  Beams..  84 

6.9  Loss  Factor  Calculations  Using  Tripled  Flexural  Rigidity  D  and  Old  FE 

Model  Data .  84 

6. 10  Loss  Factor  Calculations  Using  Tripled  Flexural  Rigidity  D,  Modified 

Finite  Element  Model  Data,  and  Experimentally  Determined  Loads  P .  85 

6. 1 1  Comparison  of  Theoretical  and  Experimental  Damping  Values, 

Conventional  Treatment  .  86 

A.l  Comparison  of  Finite  Element  Model  and  Experimental  Modal  Frequencies..  91 

C.l  3M ISD  112  Material  Property  Data  . 104 

C.2  Undamped  Modal  Loss  Factors .  104 

C.3  Experimental  Data,  11"  Beam  Experiment .  106 

1.1  Damping  Tests,  Improved  Treatment,  /2  =  0.005" .  122 

1.2  Damping  Tests,  Improved  Treatment,  t2  =  0.015" .  122 

13  Damping  Tests,  Improved  Treatment,  /2  =  0.030"  .  123 

1.4  Damping  Tests,  Conventional  Treatment,  r2  =  0.005" .  123 


VII 


List  of  Symbols 


Symbol  Definition 

Oi  -  comparison  function  coefficients 
at,  bi,  Ci,  di,  ei,  ki  -  dummy  coefficients  for  Mindlin  and  Goodman's  method 
b  -  width  of  a  sandwich  beam  (in.) 
b\  -  wide  sandwich  beam  width  (in.) 
b2  -  slim  sandwich  beam  width  (in.) 
dx  -  differential  slice  of  a  structure  or  beam  (in.) 


fi  -  time-dependent  moment  (lb.in.) 

Fg^  1 

g-  shear  parameter  =  J - (in.*i) 


gi .  comparison  function  for  the  ith  boundary  condition 
h  -  long  lever  length  (in.) 

qo  •  distributed  load  (lb7in.) 

qi .  coefficients  of  Appendix  B 
s  -  short  lever  length  (in.) 

U  -  thickness  of  the  i***  layer  (in.) 

Uo  -  longitudinal  displacement  of  the  substrate  at  the  interface  of  the  viscoelastic 
layer  (in.) 

Ui  -  longitudinal  displacement  of  the  ith  layer  (in.) 


w  - 
Bo- 


transverse  displacement  (in.) 
a  parameter  defined  as 


Bi  -  operator  for  boundary  conditions  in  Chapter  5 
C  -  coefficient  of  the  modal  expansion  in  Chapter  5 
D  -  volume  specific  energy  dissipated  per  cycle  (lb.in.) 


Ds  -  energy  dissipated  in  a  system  per  cycle  of  oscillation  (lb.  in.) 

Ej  -  Y oung's  modulus  of  the  i^  elastic  layer 

FE  -  finite  element 

(T  -  viscoelastic  layer  complex  shear  modulus  (psi) 

G'  -  viscoelastic  layer  real  shear  modulus  (psi) 

G"  -  viscoelastic  layer  loss  shear  modulus  (psi) 

Gf  -  shear  modulus  of  the  i*^  layer  (psi) 

Gn  -  coefficient  of  spatial  solution  orthogonalized  with  respect  to  MIg{(x)] 

G/u*  -  coefficient  of  spatial  solution  orthogonalized  with  respect  to  L(gi(x)] 

Hz-  hertz  (se<r*) 

/  -  imaginaiy  part  of  the  square  root  of  qg  in  Appendix  B 
L\  -  characteristic  length  of  a  constraining  layer  segment  in  Chapter  4  (in.) 

£f  jc]  -  stiffness  operator 

-  mass  operator 

Mf,  -  maximum  moment  applied  to  the  ends  of  a  sandwich  beam  during  a 
vibration  cycle 

NM)  -  generalized  force  for  the  n***  generalized  coordinate 
P  -  load  (lb.) 

Po  -  maximum  load  due  to  maximum  truss  deflection  for  a  given  cycle  of 
vibration  (lb.) 

R  -  real  part  of  the  square  root  of  qg  in  Appendix  B 
SB  -  sandwich  beam 

TJt)  -  assumed  temporal  solution  of  the  sixth-order  differential  equation 
U  -  elastic  strain  energy  (lb.  in.) 

Unom  -  nominal  strain  energy  equal  to  the  maximum  strain  energy  in  an  element  of 
a  constraining  layer  if  it  was  strained  by  e^flb.in.) 


IX 


Y-  shear  strain  (rad) 

Y  -  rate  of  change  of  shear  strain  (sec*l) 
yi,  -  shear  strain  due  to  bending  (rad) 

-  shear  strain  due  to  axial  loading  (rad) 
r]  -  loss  factor 
?)5  -  system  loss  factor 


Vv- 

X- 


viscoelastic  matenal  loss  factor 


V  -  poisson's  ratio 
0,^  -  angles  (rad) 

(lb.  sec* 


<7-  normal  stress  (psi) 

X-  shear  stress  (psi) 

(o-  frequency  (sec**) 

(OsB  -  mode  1  frequency  of  the  sandwich  beam(sec**) 
(Or  -  mode  1  frequency  of  the  truss  (sec**) 

^  -  damping  ratio  ( r;  =  20 

A  •  disp|{K;ement  (in.) 

— -  -  i***  derivative  (also  given  by  “primes”) 


X 


Abstract 


This  study  was  an  attempt  to  conceptualize,  design,  build,  and  test  a  constrained 
layer  damping  device  that  uses  mechanical  advantage  to  increase  the  damping  of  transient 
vibration  of  a  five-bay  planar  truss.  The  aluminum  truss  was  clamped  to  a  fixed  support 
and  simulated  a  solar  array.  A  conventional  constrained  layer  treatment  was  designed, 
built,  and  tested  to  serve  as  a  benchmaric  for  comparison.  For  the  conventional  treatment, 
Plunkett  and  Lee's  segmented  constraining  layer  technique  was  used. 

The  improved  device,  two  three-layer  sandwich  beams  spanning  the  two  bays 
nearest  the  built-in  end,  were  supported  several  inches  away  from  the  truss  plane  by  a 
series  of  levers  that  were  designed  to  impart  a  combined  shear  and  bending  load  to  the 
beam  ends.  Mindlin  and  Goodman's  method  was  used  to  solve  a  sixth-order, 
homogeneous  differential  equation  of  motion  with  time-dependent  boundary  conditions. 

Experimental  results  showed  that  the  sandwich  beam  device  did  not  exceed  the  loss 
factors  of  the  conventional  treatment  for  any  of  the  first  four  out-of-plane  bending  modes. 
However,  measurable  levels  of  system  damping  were  achieved  for  all  modes  of  interest, 
particularly  at  the  mode  whose  frequency  coincided  with  that  of  the  sandwich  beams. 
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PASSIVE  DAMPING  OF  A  SOLAR  ARRAY 


/.  Introduction 


Background 

Developments  in  propulsion  technology  have  now  made  it  possible  to  send  large 
structures  into  space.  The  field  of  large,  flexible  space  structures  is  one  of  growing  interest 
(the  NASA  space  station  Freedom  program  is  one  example  of  the  need  for  research  in  this 
area).  For  the  purposes  of  this  paper,  large  can  be  arbitrarily  defined  as  anything  larger  than 
the  space  shuttle.  Space  structures  are  typically  designed  to  be  as  light  as  possible,  so  an 
increase  in  mass  most  likely  means  an  increase  in  »ze.  This  increase  in  size  Imngs  with  it  a 
host  of  new  challenges  for  design  engineers.  A  primary  area  of  concern  is  the  control  of 
vibration.  Increasing  size  and  mass  creates  vibration  modes  which  can  be  very  low  in 
frequency.  The  seriousness  of  this  problem  is  amplified  by  ever  more  demanding  mission 
requirements  such  as  optical  pointing  accuracies,  spacecraft  shape  control  (of  a  parabolic 
mirror,  for  example),  and  settling  times  after  slewing  maneuvers.  Furthermore,  low  natural 
frequencies,  in  general,  cause  relatively  large  deflections,  which  give  rise  to  high  maximum 
stresses.  Consequently,  vibration  suppression  is  a  major  area  of  concern. 

By  virtue  of  their  placement  and  function  on  a  space  structure,  solar  arrays  typically  are 
problem  areas  for  vibration  control.  They  are  usually  long  and  thin,  having  low  natural 
frequencies  and  sometimes  high  modal  densities.  The  shape  of  a  solar  array  lends  itself  to 
control  by  means  of  constrained-layer  damping.  The  technology  of  constrained  layer 
damping  is  well -developed,  and  its  use  in  space  structures  is  desirable  because  passive 
damping  is  relatively  simple,  reliable,  lightweight,  and  does  not  require  power.  In  addition,  it 
is  inherently  stable  and  is  an  excellent  complement  to  active  damping.  Hnally,  passive 
damping  is  a  good  post-fabrication  fix  for  vibration  problems  that  were  not  anticipated  in 
the  design  phase.  These  characteristics  show  that  passive  damping  is  a  reasonable,  practical 
means  of  controlling  resonance  in  structures. 

For  structures  that  are  already  in  service,  the  effectiveness  of  constrained  layer  damping 
treatments  is  limited.  Ideally,  the  maximum  effectiveness  can  be  achieved  by  redesigning 
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the  structure  so  that  damping  treatments  can  be  placed  in  the  most  advantageous  locations. 
Obviously,  this  is  not  always  an  option  because,  for  a  complex  structure,  it  is  cost- 
prohibitive.  The  inability  to  optimally  locate  the  damping  material  on  a  given  structure 
provided  the  impetus  for  this  research. 

Problem 


Problem  Statement:  With  passive  means  and  by  using  an  optimized  conventional 
constrained-layer  damping  treatment  as  a  benchmark,  design,  build,  and  test  an 
alternative  damping  treatment  that  provides  more  energy  dissipation  per  unit  mass  than 
the  conventional  treatment.  Specifically,  the  loss  factor  of  the  improved  treatment  must 
exceed  that  of  the  classical  treatment  for  the  lowest  six  vibration  modes. 

A  truss  structure  was  fabricated  to  serve  as  the  test  article.  Operationally,  the  truss 
supports  photovoltaic  solar  cells  that,  for  the  purpose  of  this  research,  were  massless,  rigid, 
and  incapable  of  dissipating  energy.  Naturally,  modifications  to  the  truss  cannot  interfere 
with  the  operation  of  the  solar  cells. 

Since  additional  payload  can  cost  up  to  several  thousand  dollars  per  pound,  the  best  way 
of  damping  a  structure  may  not  be  the  simplest  or  the  cheapest.  A  reliable,  easily 
manufacturable  design  that  significantly  reduces  the  amount  of  added  weight  can  be  worth 
the  added  complexity. 

As  stated  earlier,  the  damping  material  in  a  conventional  constrained-layer  treatment  is 
not  located  optimally.  From  basic  strength  of  materials  theory,  the  maximum  shear  forces 
occur  on  the  neutral  surface  of  a  rectangular,  prismatic  beam.  Hence,  a  sandwich 
configuration  (Figure  1.1)  is  optimal. 


viscoelastic  layer 

constraining  layers 
(load-bearing  members) 


Figure  1.1  Sandwich  Beam 


Designing  an  optimal  damping  treatment  requires  an  examination  of  how  strain  energy  is 
distributed  throughout  a  structure.  Since  strain  energy  is  a  function  of  stress  and  strain,  it 
will  vary  throughout  a  member  in  bending.  In  general,  an  optimum  damping  treatment  is 
tailored  to  exploit  areas  of  high  strain  energy.  For  example,  the  strain  energy  in  a  cantilever 
beam  in  its  fundamental  mode  is  a  maximum  at  the  built-in  end. 
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Torvik  and  Lazan  [1]  have  investigated  the  possibility  of  moving  the  damping  treatment 
further  away  from  the  neutral  surface  of  the  beam  as  a  means  of  gaining  mechanical 
advantage  to  increase  the  shear  deformations  in  the  viscoelastic.  The  author  conducted 
research  in  this  area  with  emphasis  on  a  different  mechanism;  bending  of  the  viscoelastic 
layer.  See  Figure  1.2  below. 


constrained  layer  sandwich  beam  • 


structure 

Undeformed  Structure 


M 


Deformed  Structure 

Figure  1.2  Deformation  of  Structure 


The  device  is  very  simple:  when  the  structure  bends,  the  relative  rotation  of  the  levers  causes 
a  displacement  at  their  ends,  which  in  turn  causes  the  sandwich  beam  to  undergo  direct  and 
bending  shear.  Any  dimension  in  the  design  can  be  optimized  for  a  given  material  and/or 
frequency.  It  was  presumed  that  the  added  bonus  of  shear  due  to  moment  would 
significantly  increase  the  effectiveness  of  the  damping  treatment. 


Review  of  the  Literature 

The  application  of  viscoelastic  damping  treatments  is  well  documented.  Torvik's  tutorial 
work  [2]  contains  not  only  an  extensive  investigation  of  various  analysis  and  design 
techniques,  but  also  a  thorough  look  at  the  historical  developments  in  constrained  layer 
damping  over  the  years  from  1957  to  1980.  Soovere  and  Drake  [3,4]  have  compiled 
perhaps  the  most  comprehensive  work  ever  written  on  material  damping.  Their  Damping 
Design  Guide  is  full  of  design  procedures,  case  histories,  and  analysis  techniques. 

The  development  of  the  theory  of  sandwich  beams  occurred  over  a  period  of  about 
twenty  years.  E.  M.  Kerwin  [5]  is  regarded  as  the  first  person  to  derive  an  expression  for 
damping  in  a  constrained  layer  configuration.  Ungar  [6]  developed  energy  expressions  of 
Kerwin's  analysis  in  1962.  EHTaranto  [7]  was  credited  with  deriving  a  differential  equation 
of  motion  for  a  three  layer  sandwich  beam  in  1965.  The  general  boundary  conditions  for 
DiTaranto's  equation  were  developed  by  Mead  and  Markus  [8]  in  1969.  Van  and  Dowell 
[9]  followed  in  1974  by  generalizing  the  analysis  to  asymmetric  sandwich  beams.  Rao 
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[10,1 1],  among  others,  added  the  effects  of  shear  layer  extension  and  inertia  in  order  to 
develop  an  analysis  technique  for  short  sandwich  beams.  The  author  sought  to  apply  this 
body  of  knowledge  to  a  modified  geometric  configuration  that  increases  the  effectiveness  of 
the  damping  layer. 

Much  research  has  been  conducted  in  finding  more  efficient  ways  of  applying  damping 
treatments  to  structures.  Whittier  [12]  investigated  various  configurations  involving  spacers 
in  1958,  as  did  Torvik  and  Lazan  [1]  in  the  mid-1960's.  In  1965,  Metherell  and  Sokol  [13] 
examined  the  effects  of  anchoring  the  ends  in  various  ways.  In  1970,  Plunkett  and  Lee  [14] 
explored  the  effect  of  optimizing  the  length  of  the  constraining  layers  in  an  effort  to  increase 
the  strain  in  the  damping  layers.  It  was  the  author's  intention  to  investigate  the  effect  of 
pulling  on  the  ends  of  the  beam  shown  in  Figure  13,  and  optimize  this  configuration  for  the 
solar  array  truss. 


Figure  1 .3  Viscoelastic  Sandwich  Beam  With  Combined  Shear  Force  and  Moment 

Research  in  this  area  of  increasing  mechanical  advantage  ior  gains  in  damping  loss  factor 
has  generally  been  fruitful.  Some  configurations,  although  innovative  in  theory,  have  been 
difficult  to  fabricate  because  of  material  cost  and  odd  geometry.  One  of  the  challenges  of 
this  study  was  to  design  a  practical  damping  treatment  that  was  easy  to  make  and  relatively 
simple  to  analyze.  The  reason  why  simplicity  is  desirable  is  because  this  treatment  was  not 
only  designed,  but  built  and  tested  as  well. 

In  the  past  fifteen  years,  the  study  of  viscoelastic  material  behavior  has  been  extended  by 
Torvik  and  Bagley  [16]  to  include  the  use  of  fractional-order  derivative  models  that  predict 
elastomer  material  response.  The  author  did  not  use  this  analysis  technique. 


Scope 

The  experimental  nature  of  this  project  means  that  a  large  portion  of  time  was  spent  in 
testing.  Therefore,  no  deliberate  attempt  was  made  to  reach  a  theoretical  breakthrough 
regarding  the  analysis  of  sandwich  beams. 
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This  project  was  an  interdisciplinary  study  with  many  variables.  Consequently,  some 
areas  of  the  design  process  were  treated  more  rigorously  than  others.  Since  the  problem 
was  defined  arbitrarily,  there  is  no  real  reason  for  some  of  the  constraints  to  exist  other  than 
that  they  helped  define  a  practical  apparatus  to  experiment  with. 

There  are,  of  course,  many  ways  to  use  mechanical  advantage  to  increase  damping.  The 
proposed  sandwich  beam  was  chosen  for  investigation  primarily  because  it  is  simple.  Other 
configurations  can  probably  be  used  effectively  as  well,  but  even  a  simple  device  like  the 
sandwich  beam  can  present  difficult  analysis  issues  when  unusual  geometries  are  present. 

A  logical  approach  was  used  to  reduce  the  problem  to  a  woikable  level  of  complexity. 
The  five-bay  truss  was  designed  for  easy  handling  and  predictable  finite  element  analysis. 
The  specified  range  of  problem  frequencies  was  arbitrarily  defined  to  be  low.  Considering 
only  a  small  number  of  modes  to  be  dangerous  in  a  given  vibration  problem  is  a  common 
way  of  reducing  complexity.  Time  constraints  prohibited  more  than  one  new  damping 
design  to  be  built  and  tested.  Impact  testing  was  chosen  for  collecting  data  because  it  is  a 
quick,  simple  and  accurate  way  of  determining  natural  frequencies  and  modes. 

The  basic  goal  of  the  author's  research  was  to  show  that  greater  damping  efficiency  can 
be  achieved  by  judiciously  choosing  where  and  how  to  apply  a  damping  treatment  to  a 
structure.  Only  the  proposed  treatment  and  a  conventional  damping  treatment  were  tested. 

The  effectiveness  of  the  damping  treatments  was  judged  by  the  modal  loss  factors. 
Impact  testing  allowed  several  test  runs  to  be  averaged.  Averaging  served  to  minimize  the 
effect  of  noise  and  was  less  susceptible  to  error  than  a  single  test.  Five  to  ten  test  runs 
provided  sufficiently  reliable  data.  In  addition,  every  effort  was  made  to  insure  that  three 
significant  figures  of  accuracy  could  be  obtained  for  all  measurements.  Fortunately, 
accurate  calibration  of  the  test  equipment  was  done  and  this  level  of  precision  was  possible. 


Assumptions 

In  order  to  reduce  the  problem  to  a  manageable  size,  some  simplifying  assumptions  were 
made  for  the  analysis.  In  addition,  although  the  operational  environment  of  the  solar  array 
will  be  a  given  earth  orbit  (assume  the  proposed  orbit  is  geostationary),  practical 
considerations  precluded  the  use  of  the  following  design  constraints  such  as  near-vacuum 
pressures,  large,  rapid  temperature  swings,  steady  state  temperatures  well  below  CC,  and 
zero-g  effects. 
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The  first  assumption  was  that  the  truss  could  not  be  scrapped  and  redesigned.  Replacing 
the  structure  is  cost-prohibitive.  Next,  the  modification  must  have  been  producible  from 
**017  the  shelf”  materials.  Composite  materials  were  not  used  for  the  constraining  layers. 
The  assumed  operating  temperature  for  design  purposes  was  the  laboratory  temperature 
(60-85‘’F). 

The  analysis  and  design  of  the  sandwich  beams  proceeded  with  the  following 
assumptions:  Young's  modulus  of  the  viscoelastic  material  is  small  compared  to  that  of  the 
structure  and  constraining  layers;  material  thicknesses  remain  constant  throughout  all 
deformations;  the  bonds  between  the  different  layers  are  perfect  (so  that  forces  could  be 
considered  to  be  distributed  over  the  entire  area);  and  published  viscoelastic  material 
property  data  are  accurate. 


Approach 

The  research  project  was  divided  into  the  following  major  steps: 

1. )  Design  a  test  article  that  is  simple  enough  to  characterize  quickly. 

2. )  Calibrate  the  test  equipment 

3. )  Conduct  a  coarse-mesh  finite  element  (FE)  analysis.  Refine  the  model  and  analyze 
again. 

4. )  Characterize  the  actual  structure  when  received.  Find  frequencies,  mode  shapes, 
and  modal  loss  factors. 

5. )  Establish  test  conditions  in  order  to  obtain  design  constraints  (testing  temperature 
is  one  example  of  a  test  condition). 

6. )  Design  a  conventional  constrained-layer  treatment.  Explore  different  techniques 
and  optimize  the  design. 

7. )  Use  the  theoretical  results  obtained  from  the  conventional  treatment  as  a  benchmaik 
to  establish  design  specifications  for  the  improved  treatment.  Design  the  improved 
treatment  with  the  constraint  that  it  must  be  able  to  be  removed. 

8. )  Produce  design  drawings  for  the  AFTT  Model  Fabrication  Shop  and  have  both 
treatments  built 

9. )  Install  the  improved  treatment  first,  and  then  characterize  the  improved  damped 
system. 

10. )  Remove  improved  treatment  and  apply  the  conventional  treatment 

11. )  Characterize  the  system  again. 
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12. )  Write  a  lab  report  that  includes  a  comparison  of  the  both  sets  of  results.  Discuss 
uncertainty,  acknowledge  unfairness  in  the  comparison,  and  propose  further  improvements. 

13. )  Write  the  report. 


Sequence  of  Presentation 

This  project  is  a  bold  undertaking  whose  approach  and  schedule  had  no  precedent 
Therefore,  the  order  of  presentation  is  similar  to  the  order  of  the  research.  Since  no  design 
was  predetermined  at  the  outset  of  the  analysis,  a  chapter  was  needed  to  include  some 
general  ideas  on  improving  a  conventional  constrained  layer  treatment  Chapter  three 
explains  how  the  test  article  was  designed  and  characterized  (modes  and  frequencies 
determined).  The  next  chapter  includes  the  analysis  and  design  of  the  conventional 
damping  treatment  In  (Tiapter  five,  a  theory  is  developed  for  analyzing  the  proposed 
improved  damping  treatment  a  sandwich  beam  subjected  to  time-dependent  boundary 
conditions.  Expressions  for  the  shear  strain  and  loss  factor  follow  in  this  section.  Finally, 
the  design  of  the  improved  treatment,  and  construction  and  testing  of  both  treatments  is 
discussed. 
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//.  Improving  On  A  Conventional  Constrained  Layer  Damping  Treatment 


In  researching  the  many  different  applications  of  constrained  layer  damfnng,  it  became 
evident  that  many  treatments  possess  the  same  problem:  uniform  treatments  are  applied 
over  distances  where  strains  change,  sometimes  dramatically.  Designing  variable-thickness 
layers  for  uniform  stress  is  possible,  but  fabricating  such  a  treatment  is  a  non-trivial  matter. 
Given  the  constraints  regarding  the  manufacture  of  constrained  layer  damping  treatments, 
compromises  must  be  made  when  choosing  lengths  and  thicknesses.  This  unfortunate  fact 
warrants  a  careful  look  into  how  constrained  layer  damping  can  be  made  more  effective. 

The  amount  of  energy  dissipated  per  cycle  of  vibration  per  unit  volume  is 


D-itrydt  (2.1) 

where 

r-  shear  stress, 
y  *  rate  of  shear  strain,  and 
(O  -  natural  frequency  of  the  vibrating  structure  . 

The  natural  frequency  of  a  given  structure  will  change  a  relatively  small  amount  with  the 
addition  of  mass  and/or  stiffness.  It  is  a  parameter  of  which  we  have  relatively  little  control 
unless  major  changes  are  made  to  the  structure.  The  shear  stress  is  a  forced  condition 
based  on  geometric  deformations.  On  the  other  hand,  the  rate  of  shear  strain  is  based  on  the 
rate  of  change  of  the  structure's  shape.  The  distribution  of  strain,  and  thus  strain  energy,  is 
certainly  not  uniform  for  any  structure  unless  it  is  exceedingly  basic  (a  simply  supported 
beam  with  circular  curvature  is  one  important  example).  Therefore,  the  value  of  shear  strain 
in  a  structure  depends  a  great  deal  on  where  we  look.  It  is  possible  to  exploit  large  shear 
strains  by  carefully  choosing  where  and  how  to  install  a  damping  device.  If  we  assume  that 
the  operating  environment  of  the  structure  cannot  be  changed,  we  then  look  to  design  a 
damping  treatment  with  two  important  characteristics:  the  ability  to  generate  high  rates  of 
shear  strain,  and  the  ability  to  evenly  distribute  high  strains  over  significant  areas.  In  this 
way,  we  can  obtain  maximum  effectiveness  from  a  uniform  coverage.  Clearly,  the 
mechanism  to  employ  is  strain  amplification  by  mechanical  advantage. 

Any  given  structure  has  many  vibration  modes  and  corresponding  natural  frequencies. 
It  is  important  to  identify  structural  modes  that  a.)  have  the  possibility  of  being  excited  (by 
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movement  of  the  structure,  rotation  of  an  on-board  motor,  or  some  other  type  of  forcing 
function),  and  b.)  have  the  possibility  of  causing  damage  if  a  resonance  condition  goes 
unchecked.  For  most  metals,  very  high  frequency  vibrations  are  dissipated  by  natural 
acoustic  radiation.  In  space,  this  range  of  frequencies  is  a  problem  because  there  is  no 
medium  in  which  to  transfer  energy.  For  this  paper,  high  frequencies  will  not  be  considered 
a  problem.  In  general,  frequencies  too  low  to  acoustically  radiate  yet  too  high  to  visually 
observe  are  typically  dismissed  because  every  material  has,  to  some  extent,  a  material  loss 
factor  (albeit  small)  which  causes  energy  dissipation  in  a  short  time  simply  because  the 
period  of  oscillation  is  short.  An  important  exception  to  this  case  is  the  operation  of  a 
slightly  unbalanced  motor  maintaining  a  given  rate  of  rotation.  Barring  this  exception,  the 
frequencies  typically  considered  to  be  of  concern  are  the  lowest  ones.  There  are  three  major 
reasons  for  this: 

1. )  Movement  of  a  structure  is  a  common  occurrence  and  tends  to  excite  low  modes, 

2. )  Low  frequency  modes  produce  the  largest  displacements  and  stresses. 

3. )  By  virtue  of  their  relatively  long  periods  of  vibration,  low  modes  have  the  longest 

settling  times.  This  becomes  a  problem  when  performance  specifications  require  a 
short  settling  time. 

These  facts  emphasize  the  need  to  control  low  modes. 

The  geometry  of  a  complex  structure  will  most  likely  have  one  or  more  regions  of 
relatively  high  strains.  For  this  paper,  the  structure  of  interest  is  a  two-dimensional  truss 
that  simulates  a  solar  array.  Because  the  energy-collecting  cells  it  holds  are  most  likely 
brittle,  the  array  is  designed  to  prohibit  movement  in  the  plane  (shear  and  bending). 
However,  out-of-plane  bending  and  torsion  modes  generally  have  low  frequencies  and 
relatively  large  displacements.  The  movement  of  the  stracture  in  these  modes  allows  several 
ways  of  using  mechanical  advantage  to  increase  damping.  A  few  novel  ideas  were 
proposed,  but  some  of  them  would  require  a  considerable  support  structure.  Most  of  the 
design  effort  in  this  project  centered  around  trying  to  exploit  a  large  distance  change 
between  two  points  during  a  typical  structural  deformation. 

Generating  bending  in  some  kind  of  sandwich  beam-like  addition  seemed  to  offer  the 
most  versatile  design  possibilities.  It  then  became  apparent  that,  since  this  "improved” 
treatment  would  most  likely  be  attached  at  discrete  points,  the  most  fundamental  concern 
was  getting  as  much  displacement  or  rotation  out  of  the  attachment  points  as  possible.  With 
a  possible  solution  at  hand,  the  focus  became,  “What  will  this  improved  treatment  be 
compared  to?”  At  this  point,  the  course  of  action  for  the  project  was  defined.  The  research 
would  proceed  through  the  following  phases: 
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1. )  Design  of  a  realistic  test  article.  It  needed  to  be  easy  to  make,  relatively  simple  to 
analyze,  and  easy  to  modify. 

2. )  Second,  two  unknown  parameters  had  to  be  determined  as  constraints  for  an 
improved  design:  maximum  weight  addition  and  a  target  system  loss  factor.  The  only 
reasonable  way  to  find  these  two  numbers  was  to  design  a  conventional  constrained  layer 
treatment 

3. )  Once  the  weight  and  loss  factor  are  determined,  conceptualize,  design,  and  analyze 
an  improved  damping  treatment 

4. )  Finally,  the  two  designs  need  to  be  built  and  tested.  At  that  point,  results  can  be 
compared  and  any  testing  tnases  acknowledged. 

Once  the  general  sequence  of  research  was  determined,  it  was  time  to  proceed  to  the 
design  of  the  test  article;  the  solar  array  truss. 
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HI.  Design  of  the  Test  Article 


The  inspiration  for  the  test  article  originated  from  the  observation  of  the  Passive  and 
Active  Control  of  Space  Structures  (PACOSS)  Representative  System  Article  (RSA)  data 
[16].  Martin  Marietta's  data  on  the  low  frequency  modal  identification  showed  that  159  out 
of  the  structure's  lowest  210  modes  (all  below  10  Hz)  were  solar  array  modes.  Since  the 
thrust  of  the  proposed  research  was  the  design  and  testing  of  a  damping  treatment, 
designing  a  structure  with  this  kind  cS  modal  density  would  needlessly  complicate  matters. 
Since  the  number  of  system  modes  (tecreases  geometrically  with  decreasing  complexity,  the 
best  way  to  avoid  high  modal  density  is  to  make  it  with  fewer  members.  In  addition,  most 
arrays  are  attached  to  a  space  structure  with  a  single  shaft  which  allows  rotation  about  their 
longitudinal  axis  for  the  purpose  of  pointing  the  array  with  a  minimum  amount  of  energy. 
Converting  the  single  support  point  to  two  attached  beams  serves  to  “shiff*  the  torsional 
vibration  modes  to  higher  frequencies. 

Since  the  test  article  is  supposed  to  be  a  preexisting  structure,  every  effort  was  made  to 
make  it  as  realistic  as  possible.  However,  some  desirable  features  needed  to  be  designed  in 
to  make  laboratory  work  easier.  For  example,  low  frequencies  are  visually  observable,  and 
being  able  to  actually  see  the  various  modes  helps  give  one  a  ‘Teel”  for  what  is  happening. 
Additionally,  an  abundant  material  such  as  aluminum  makes  procurement  and  fabrication 
easy.  Furthermore,  some  FE  analysis  was  planned,  so  a  high  level  of  idealization  was 
desired.  Fully  clamped  ends,  rectangular,  prismatic  cross-sections,  and  welded  connections 
were  convenient  attributes  that  allowed  easy  finite  element  modeling. 

Using  these  considerations,  the  final  design  for  the  test  article  was  determined.  The 
article  is  a  two-dimensional,  five  bay  truss  whose  dimensions  are  listed  on  the  following 
diagram.  For  an  oblique  close-up  view  of  the  clamped  area,  the  reader  is  referred  to  Chapter 
6. 


ir 


Figure  3.1  Test  Article  With  Dimensions 
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The  long  members  are  continuous  beams  whose  cross-section  is  detailed  on  the 
figure.  The  cross-members  are  made  of  the  same  cross-section  and  are  butt  jointed  and 
welded  to  the  long  members.  The  weld  bead  extends  a  small  distance  (about  1/16")  into  the 
long  members  and  was  ground  smooth  to  their  edge.  The  cantilever  connection  points 
consist  of  two  pipe  clamps  made  from  0.500*  plate  aluminum.  They  were  secured  to  a 
6.75"  outer  diameter,  6.25"  inner  diameter  steel  tube.  There  is  a  thin  aluminum  strip  (.040*) 
that  lines  the  inside  of  each  clamp  in  order  to  insure  a  very  snug  fit.  Finally,  the  securing 
plates  depicted  in  the  figure  are  the  same  thickness  as  the  truss  cross-section  and  have  the 
purpose  of  stiffening  up  the  clamped  end.  The  pipe  only  had  about  ten  inches  of  clearance 
for  clamping.  Thus  the  reason  for  the  complicated  connection  points.  The  pipe  is  itself  a 
cross-member  in  one  of  PACOSS's  three  suspension  supports.  These  supports  are  about 
20  feet  high  and  are  secured  to  the  cement  floor  with  several  large  bolts.  In  addition,  each 
support  has  a  few  hundred  pounds  of  lead  on  their  bottoms  to  further  secure  the  supports. 
The  transmission  of  energy  through  the  supports  to  the  truss  is  only  possible  at  high 
frequencies  (  greater  than  about  200  Hz)  and  can  be  accomplished  by  striking  the  support 
with  a  metal  object.  However,  no  amount  of  human  force  applied  to  the  support  caused  the 
truss  to  move  in  any  way. 

From  this  point  on,  the  solar  array  will  be  referred  to  as  the  "'truss." 


Characterization  of  the  Truss 

Once  the  truss  was  received  and  mounted,  its  modal  frequencies  needed  to  be 
determined.  First,  the  equipment  needed  to  be  calibrated  over  the  range  of  frequencies  of 
interest.  The  most  important  consideration  in  calibrating  the  equipment  was  making  sure 
the  trials  were  repeatable.  The  goal  of  the  testing  was  to  determine  the  ratio  of  output 
voltage  (acceleration)  to  input  voltage  (impact  force)  across  the  frequency  spectrum.  This 
transfer  function  would  be  used  to  determine  modal  frequencies  and  damping  ratios. 


Summary  of  Equipment  Used 

The  accelerometer,  a  Kistler  PeizoBEAM^’*  (model  8632A5),  was  calibrated  with 
respect  to  a  standard  accelerometer  calibrated  by  the  National  Bureau  of  Standards.  The 
output  voltage  was  calibrated  at  100  Hz  and  measured  99  mV/g.  The  force  transducer 
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(Kistler  model  9712A5)  was  calibrated  at  840  mV/lbf.  The  computer  software  used  was 
Tektronix  TekSTAR™  version  2.2.  This  software  has  an  averaging  feature  whereby  several 
response  measurements  can  be  averaged  to  smooth  out  data.  Simple  “boxcar'’  windowing 
was  used  to  obtain  data.  The  software  also  includes  a  feature  where  the  user  can  choose  the 
number  of  samples  per  trial,  up  to  a  maximum  of  4096.  Choosing  the  frequency  band  ci  0- 
20  Hz  and  miming  the  trial  for  4096  samples  made  for  a  71. 1  second  trial.  The  effect  of  the 
60  Hz  noise  spikes  was  kept  to  a  negligible  level  by  choosing  the  maximum  number  of 
samples  per  trial.  The  reason  why  this  noise  had  to  be  kept  to  a  minimum  was  because  the 
spikes  caused  the  polynomial  curve  fitter  (a  feature  of  the  SMS-StarStmct™  software 
package  for  manipulating  data)  to  form  less  abmpt  modal  peaks  (an  indication  of  higher 
damping). 


Summary  of  Truss  Dynamics 

Once  the  equipment  was  calibrated,  the  accelerometer  was  mounted  to  the  truss.  In 
order  to  fully  characterize  the  dynamic  response  of  the  tmss,  the  accelerometer  had  to  be 
mounted  at  more  than  one  location.  For  example,  the  second  out-of-piane  bending  mode 
has  a  node  at  about  0.8L.  This  location  happened  to  coincide  with  a  cross  member,  making 
it  possible  to  excite  the  fundamental  mode  almost  exclusively.  The  tip  of  the  truss  was  an 
excellent  accelerometer  location  for  recording  all  modes,  and  a  large  number  of  excitation- 
response  configurations  had  the  accelerometer  at  the  truss  tip. 

For  each  data  point,  ten  trials  were  averaged  and  the  corresponding  spectral  transfer 
function  was  saved  as  a  TekSTAR  data  file.  Then  the  modal  peaks  were  curve  fit  using  the 
SMS  StarStruct  software.  The  output  was  given  in  percent  damping,  so  the  value  had  to  be 
doubled  to  obtain  the  modal  loss  factor  (according  to  rf=  2(D> 

Some  modes  were  quite  close  to  one  another  in  frequency,  and  determining  the  mode 
type  was  done  by  moving  the  accelerometer  to  a  different  location.  For  example,  placing  the 
accelerometer  on  the  longitudinal  centerline  of  the  truss  should  yield  a  response  devoid  of 
torsional  modes.  Then,  by  moving  the  accelerometer  off  the  centerline,  all  of  the  new  peaks 
ought  to  represent  torsional  modes.  A  more  complete  characterization  that  includes  the 
undamped  loss  factors  appears  in  Appendix  A,  but  the  frequencies  of  the  modes  of  interest 
are  listed  below. 
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Mode 

Frequency  (Hz) 

Mode  Type 

1 

1.06 

2 

6 

BerKling  2 

3 

7.4 

Torsion  1 

4 

16.5 

Bending  3 

5 

22 

In-Rane  Bending  1 

6 

22.5 

Torsion  2 

7 

31.8 

Bending  4 

Table  3.1  Experimental  Modes  and  Frequencies 


A  finite  element  model  was  constructed  for  comparison.  An  accurate  FE  model  was 
needed  because  analysis  of  the  distribution  of  strain  energy  would  be  used  in  designing  the 
improved  damping  treatment.  Since  the  structure  was  designed  to  be  easily  idealized, 
modeling  was  quick  and  simple.  The  input  dimensions  were  accurate  to  four  significant 
figures  on  cross-sectional  area  and  the  area  moments  of  inertia,  and  to  within  three  figures 
on  distances.  The  modal  analysis  on  the  FE  model  produced  the  following  results: 


Mode 

Frequency  (Hz) 

Mode  Type 

1 

2 

&94 

Bending  2 

3 

7.90 

Torsion  1 

4 

16.5 

Bending  3 

5 

23.2 

In-Rane  Bending  1 

6 

24.2 

Torsion  2 

7 

31.7 

Bending  4 

Table  3.2  Rnite  Element  Model  Modes  and  Frequencies 


These  results  compare  favorably  to  the  actual  structure,  which  suggests  that  the  irrodel 
is  a  good  one.  In  addition,  an  accurate  model  allows  one  to  iterate  through  a  number  of 
proposed  improvements  quickly,  and  thus  saving  a  great  deal  of  time.  Without  a  usable 
model,  the  improved  treatment  would  have  had  to  have  been  made  adjustable  in  order  to 
iterate  to  an  optimum  configuration. 

Once  the  test  article  was  measured  and  its  modal  frequencies  determined,  the  next  step 
was  the  design  of  the  conventional  treatment 
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IV.  Conventional  Constrained  Layer  Treatment  -  Analysis  and  Design 


The  purpose  of  analyzing  and  designing  a  conventional  constrained  layer  damping 
treatment  was  to  obtain  two  key  design  constraints  for  the  improved  treatment:  the 
maximum  weight  and  the  target  modal  loss  factors. 

An  extensive  review  of  the  earlier  literature  (1970  and  previous)  was  conducted  in 
order  to  find  an  analysis  technique  suited  to  the  problem.  Most  analyses  on  this  subject 
can  be  derived  from  E.  M.  Kerwin's  1959  paper  [5].  Two  methods  of  analysis  were 
chosen  for  this  project:  Torvik's  variational  development  [17],  and  Plunkett  and  Lee’s 
approach  which  assumes  constant  strain  in  the  substrate  [14].  Torvik's  1969  paper,  using 
the  calculus  of  variations,  is  a  good  teaching  tool  because  of  its  elegance  and  because  the 
boundary  conditions  are  a  natural  by-product  of  the  mathematics.  However,  after  finding 
the  expression  for  the  loss  factor  and  substituting  in  the  truss's  dimensions,  it  was 
determined  that  a  constraining  layer  length  of  71"  (the  length  of  the  long  truss  members) 
was  too  long.  For  this  reason,  the  method  of  Plunkett  and  Lee  was  chosen  as  the 
theoretical  basis  for  the  conventional  design.  The  development  of  the  loss  factor  by  the 
variational  technique  appears  in  Appendix  B. 

The  theory  of  Plunkett  and  Lee  is  repeated  here  for  convenience.  They 
hypothesized  that  there  exists  an  optimum  length  constraining  layer  segment  which  gives 
maximum  damping  for  a  given  viscoelastic  material  and  a  given  frequency.  Their 
explanation  of  this  hypothesis  is  the  following: 

If  the  constraining  layer  is  very  long,  the  shear  stress  near  the  ends  induces  the  same 
axial  strain  in  it  as  in  the  basic  structure;  thus  there  is  no  shear  in  the  viscoelastic 
layer  away  from  the  ends  and  the  damping  is  small.  If  the  length  of  each  element  of 
the  constraining  layer  is  very  short,  it  exerts  no  constraint  on  the  underlying 
viscoelastic  layer,  there  is  little  shear  strain,  and  the  damping  is  small.  At  some 
finite  value  for  the  lengths  of  the  elements  of  the  constraining  layer,  the  damping  is 
a  maximum.  [14:  151] 

Plunkett  and  Lee’s  theory  is  based  on  the  following  assumptions: 

1. )  All  energy  dissipation  in  the  viscoelastic  layer  is  attributable  to  shear  strain. 

2. )  The  constrained  and  constraining  layers  have  small  thicknesses  with  respect  to 
the  substrate.  This  assumption  a.)  makes  bending  stiffnesses  in  the  constraining  layer 
negligible  with  respect  to  the  structure,  and  b.)  allows  us  to  analyze  the  constraining 
layer  as  if  it  were  in  tension  only,  and  the  viscoelastic  layer  as  if  it  were  in  shear  only. 
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3.)  The  viscoelastic  material  behaves  linearly  throughout  the  range  of  operating 
strains  and  can  be  characterized  by  its  shear  modulus,  the  complex  value 

G* -G'  +  iG"-G'(l+in,)  (4.1) 


where 


G'- 

G"- 

n»- 


the  elastic  shear  modulus , 
the  loss  shear  modulus,  and 

G” 

the  material  loss  factor  - — . 

G' 


4. )  The  constraining  layer  and  the  substrate  are  elastic  and  dissipate  no  energy. 
This  is  a  fairly  good  assumption  since  aluminum  has  a  structural  damping  factor  of  Tfs  % 
0.0005. 

5. )  Poisson  effects  are  negligible. 

6. )  The  axial  strain  is  uniform  at  the  interface  of  the  substrate  and  the  viscoelastic 
layer.  For  a  beam,  this  assumption  implies  circular  curvature.  Even  in  deformations 
where  this  assumption  is  not  explicitly  the  case,  the  strain  in  the  substrate  will  approach  a 
constant  value  as  the  segment  length  gets  small. 

7. )  The  shear  strain  in  the  viscoelastic  layer  is  constant  through  the  thickness. 

8. )  The  normal  strain  in  the  constraining  layer  is  assumed  uniform  through  the 
thickness. 

9. )  Young's  moduli  of  the  constraining  layer  and  substrate  are  much  greater  the 
shear  modulus  of  the  viscoelastic  layer. 

Suppose  we  have  a  given  constraining  layer  segment  like  the  one  in  figure  4.1. 


dr 


Figure  4.1  Depiction  of  Constraining  Layer  Segment  With 
Free  Body  Diagram  of  Slice  dx 


The  equilibrium  equation  for  the  differential  element  dx  is 
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(4.2) 


-Ij-T 


Using  assumption  7,  the  constitutive  relations  for  the  constraining  layer  and  the 
viscoelastic  layer  ate 


a  -  Ee.  ■  E — 
dx 


X  -Gy  —  G 


respectively.  Assuming  constant  strain  in  the  substrate, 


Substituting  eqs.  (43),  (4.4),  and  (4.5)  into  (43)  gives 


d^U  m(  U  —U 

vX  I 


UjE,  <?V 

For  future  convenience,  use  Plunkett's  notation  to  define 


Bi.m 


The  boundary  conditions  for  any  single  constraining  layer  segment  are 


The  general  solution  of  this  equation  is 
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u  (ar)-c^  jr-4 


sinh 

U/. 

rzri 

cosh 

2Er 


Now,  using  assumption  7  again,  the  energy  dissipated  per  cycle  per  unit  width  and 
length  is 


dD  - 


(4.10) 


where 


Using  the  relation 


u  -Up 
U 


(4.11) 


gives 


dD-m,G'^U„  x-b: 


sinhi 


sinh 

[  X  j 

coshj 

'a.] 

(4.12) 


Integrating  D  over  the  length  of  the  constraining  layer  gives  the  energy  dissipated 
per  cycle  per  unit  width. 


D  «  ^jtBg  tjEjBg 


sinh  -^cos-^  sin— -sin  -^sin—  cos— 

U  2j  2  [b:  2}  2 

- - Tl  — rTi - IL - 1 - 

cosh  -^cos-^  -fcos  -^sin-^ 


(4.13) 


b:  2 


b:  2 


We  now  define  a  nominal  energy  (/no^  which  equals  the  maximum  strain  energy  in 
an  element  of  the  constraining  layer  if  it  was  strained  by  as 
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(4.14) 


U. 


We  can  nondimensionalize  D  by  dividing  by  Unom>  Define  a  dimensionless  loss 
coefficient  rjias 


(4.15) 


This  value  is  found  to  be  a  maximum  at  Li/Bq  =  3.28.  Therefore,  the  optimum  length  Li 
=  3.28Bo. 

Plunkett  and  Lee  used  a  cantilever  beam  coated  on  both  sides  (the  same 
configuration  used  here)  as  an  application  of  their  theory.  Their  analysis  of  the  cantilever 
beam  follows. 

Consider  the  deformations  of  a  cantilever  beam  subjected  to  sinusoidal  flexural 
oscillation.  The  strain  at  die  interface  is 


e 


O 


2  dx‘ 


(4.16) 


where 


/,  -  thickness  of  the  substrate 

curvature  of  beam 
dx^ 


Assuming  the  damping  treatment  is  applied  from  0  to  L,  the  energy  dissipated  per 
cycle  in  the  system  Ds  is 


D. 


(4.17) 


where  the  2  comes  from  treating  both  sides  of  the  beam.  Substituting, 
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The  maximum  strain  energy  in  the  system  is 


The  system  loss  factor  rjs  is  then 


(4.18) 


(4.19) 


(4.20) 


Assuming  a  given  segment  length  is  much  less  than  a  wavelength  and  the  treatment 
is  applied  over  the  entire  length  of  the  beam,  the  integrals  are  identically  equal  at  all 
locations  jc.  For  £[  =  £3, 


sinh 

(  j  \ 

-Icos— 

[b:  'll 

.  IJv  .  1 

Sin— -Sin 

2  1 

-^sin— 

.B.  2) 

COS— 

2 

cosh 

4co,i 

.t-cos 

l».  2, 

1 

(4.21) 


A  parametric  study  was  conducted  on  a  spreadsheet  to  find  the  optimum 
configuration  of  thicknesses  and  segment  length.  Since  G*  is  dependent  on  frequency, 
each  mode  will  exhibit  a  different  optimum  length  of  segment.  Any  given  length,  then, 
will  have  different  loss  factors  for  each  mode.  A  graph  comparing  the  theoretical  loss 
factor  and  segment  length  was  constructed  for  a  number  of  candidate  configurations. 
Figure  4.2  is  one  example  of  such  a  graph.  The  curves  for  the  higher  modes  are  less 
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Loss  Factor  Ela 


accurate  because  the  assumption  that  the  segment  length  is  much  shorter  than  the 
wavelength  is  no  longer  valid. 


Figure  4.2  Theoretical  Loss  Factor  vs.  Segment  Length  for  032005  IS0 112 


Before  a  configuration  was  chosen,  a  short  laboratory  experiment  was  conducted  to 
verify  whether  or  not  this  theory  would  match  the  results  of  an  actual  experiment  Three 
identical  1 1"  A1 6061-T6  beams  with  cross-section  0.126"  x  0.75"  were  used  to  perform  a 
number  of  tests,  each  test  having  a  different  segment  length.  Each  test  consisted  of 
collecting  data  on  several  sets  of  trials.  The  damping  material  was  a  0.002"  thick  3M ISD 
1 12  viscoelastic  layer  with  a  0.005"  thick  type  1 100  dead  soft  aluminum  constraining 
layer.  Each  segment  length  was  tested  by  means  of  conducting  ten  imp^t  trials  averaged 
into  a  single  frequency  response  function  (FRfO.  which  was  then  curve  fit  for  the 
damping  ratio.  The  damping  ratio  (in  percent  damping)  was  then  multiplied  by  two  to 
obtain  the  experimental  system  loss  factor.  Each  set  of  ten  impact  trials  produced  two 
data  points  (for  modes  1  and  2),  which  were  plotted  in  Figure  43.  The  details  of  this 
“mini-experiment”  and  the  tabulated  data  can  be  found  in  Appendix  C. 
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Figure  4.3  11"  Beam  Experimental  Data  With  Theoretical  Curves 

Mantena,  Gibson,  and  Hwang  [18]  ran  a  similar  experiment  with  3M  ISD  1 12  and 
concluded  that  the  dead  soft  aluminum  backing  (theoretically  assumed  to  have  negligible 
energy  dissipation)  had  a  significant  material  loss  factor.  3M  does  not  test  ISD  112  with 
the  aluminum  backing,  so  loss  data  was  not  available.  If  the  material  loss  factor  is  indeed 
significant,  the  experimental  loss  factor  would  be  higher,  particularly  for  long  segment 
lengths  (where  strain  in  the  constraining  layer  is  close  to  strain  in  the  substrate  over  a 
considerable  distance).  This  phenomenon  would  explain  why  the  long  segment  lengths 
gave  considerably  higher  damping  than  the  theory.  Based  on  the  data  above,  it  was 
concluded  that  Plunkett  and  Lee's  theory  is  verifiable.  Furthermore,  by  eliminating  as 
much  uncertainty  as  is  reasonably  possible,  obtaining  repeatable  data  is  possible. 

Once  the  preceding  experiment  was  finished,  a  configuration  was  chosen.  Several 
factors  were  considered  in  the  decision.  Within  reason,  the  potential  loss  factor  is  limited 
only  by  the  constraining  layer  thickness.  Hence,  a  restriction  was  placed  on  the  weight  of 
the  entire  treatment.  A  maximum  weight  of  one  pound  (about  V4  of  the  truss's  weight) 
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was  arbitrarily  chosen.  The  limit  on  the  constraining  layer  thickness,  then,  is  about  V32" 
(0.032").  Aluminum  of  this  thickness  is  produced  in  sheet  form.  The  added  uncertainty 
of  using  the  dead  soft  aluminum  backing  precluded  its  use  (even  though  it  would  be 
advantageous  to  have  another  energy-dissipating  component).  Therefore,  bare  ISD  112 
was  selected.  Table  4.1  shows  the  change  in  optimum  length  for  varying  t2  (while  all 
other  parameters  are  held  constant).  As  observed  from  Table  4. 1,  if  the  segment  length  is 
allowed  to  vary  with  varying  t2,  then  the  loss  factor  is  only  iruUrectly  dependent  on  (2-  It 
would  be  advantageous,  then,  to  select  the  thinnest  t2  because  it  is  the  lightest.  In 
addition,  the  table  shows  how  the  segment  length  decreases  with  decreasing  t2- 
Intuitively,  since  the  distribution  of  strain  along  a  cantilever  beam  varies,  decreasing  the 
segment  length  (to  a  reasonably  small  len^h)  should  give  an  overall  higher  damping 
factor.  If  the  strain  decreases  monotonically  (as  in  mode  1  vibration),  shorter  segments 
will  better  exploit  the  concentrations  of  strain  and  give  rise  to  higher  damping  because 
energy  dissipation  is  proportional  to  the  square  of  the  shear  strain.  Besides,  one  of  the 
initial  assumptions  was  that  the  segment  length  is  much  smaller  than  a  wavelength. 
Theoretically,  then,  the  optimum  value  of  L\  is  probably  small.  Unfortunately,  a 
configuration  with  a  small  segment  length  would  require  a  viscoelastic  layer  of  vanishing 
thickness.  The  following  are  nominal  pararameters  used  to  construct  Table  4.1  (t;i  was 
determined  from  the  other  parameters): 


tl  =0.195",  £=107psi,  G  =  23psi,  ijv  =  0.61,  w  =  i'l/Bo  =  3.28,  m  =0-7495 


b 

»?svs 

0.002 

0.032 

15.987 

0.0625 

0.1174551 

0.004 

0.032 

22.608 

0.125 

0.1174551 

0.005 

0.032 

25.277 

0.1563 

0.1174551 

0.007 

0.032 

29.908 

0.2188 

0.1174551 

0.009 

0.032 

33.912 

0.2813 

0.1174551 

0.01 

0.032 

35.747 

0.3125 

0.1174551 

0.012 

0.032 

39.159 

0.375 

0.1174551 

0.015 

0.032 

43.781 

0.4688 

0.1174551 

Table  4.1  Example  of  Relationship  Between  Li  and  ^for  Nominal  Parameters 


Another  unfortunate  fact  is  that  the  contact  surfaces  sandwiching  the  viscoelastic 
must  have  some  nonzero  value  of  average  surface  roughness  to  enable  adequate  bonding 
of  the  ISD  1 12.  The  problem  of  surface  roughness  could  not  be  ignored  because  its  value 
on  both  the  truss  and  the  constraining  layer  is  on  the  order  of  the  thinnest  available  t2 
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(0.002"  is  3M's  thinnest  layer).  This  roughness  was  deliberately  exacerbated  by  sanding 
the  surfaces  to  insure  a  good  bond.  Since  the  shear  stresses  are  so  sensitive  in  this  region 
of  values  of  t2,  it  was  decided  that  selecting  /2  =  0.002”  would  be  unwise.  However,  this 
effect  can  be  diminished  markedly  by  going  to  the  next  thinnest  available  layer. 
Consequently,  a  viscoelastic  layer  thickness  (2  =  0.005”  was  chosen. 

According  to  Figure  4.2,  the  optimum  segment  length  differs  for  each  mode.  The 
mode  1  optimum  segment  length  of  253  in.  would  be  a  poor  choice  because  the 
effectiveness  of  that  length  at  higher  modes  is  low.  As  the  graph  indicates,  a  much  better 
choice  of  length  is  somewhere  between  the  peaks  for  modes  1  and  2.  In  order  to  obtain 
higher  damping  in  the  higher  modes,  a  length  of  18”  was  chosen.  As  you  can  see,  this 
length  is  closer  to  the  mode  2  peak,  but  the  relative  insensitivity  of  the  mode  1  curve  to 
changes  in  length  in  this  region  justifies  the  choice.  Therefore,  the  final  configuration  for 
the  conventional  treatment  is  the  following: 

Constraining  Layer  Thickness  ti  =0.032” 

Viscoelastic  Layer  Thickness  t2  =  0.(X)5" 

Constraining  Layer  Length  Li  =  18.0” 

Based  on  Rgure  4.2,  the  expected  values  of  rj  for  modes  1  and  2  are; 

wei  ght  =  0.96  lb. 

Mode  1:  rj  =  0.11 
Mode  2:  ty  =  0.14 

These  expected  values  became  inequality  design  constraints  for  the  improved 
damping  treatment. 
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V.  Improved  Damping  Treatment 


The  major  goal  of  devising  an  improved  damping  treatment  for  a  cantilevered 
structure  is  exploiting  regions  of  high  strain  energy.  For  low  frequency  vibration  modes, 
this  region  is  the  built-in  end.  Somehow,  a  damping  treatment  must  be  designed  so  that 
strain  energy  is  “channeled”  into  an  attached  device.  Unless  this  device  contacts  the 
structure  over  a  large  area  (like  a  conventional  treatment),  it  will  have  to  be  secured  to  the 
structure  at  discrete  points.  This  observation  leads  one  to  consider  which  point  properties 
are  important  The  two  properties  which  are  useful  in  this  case  are  displacements  and 
forces.  If  we  assume  that  mode  1  and  possibly  mode  2  are  problem  frequencies,  then  we 
can  reduce  the  problem  to  a  reasonable  level  of  complexity.  At  this  point  it  appeared 
that  the  goal  of  the  thesis  statement  (obtaining  higher  modal  loss  factors  of  the  first  six 
vibration  modes)  made  for  a  design  problem  whose  solution  would  be  too  complex.  For 
the  improved  treatment,  it  was  decided  to  design  such  that  only  mode  1  would  be 
considered.  This  concession  is  not  an  admission  of  defeat:  rather,  it  is  a  practical 
reduction  of  the  problem  in  order  to  obtain  a  solution.  By  designing  to  damp  mode  1 
vibration  most  effectively,  modes  2  through  6  may  in  fact  have  higher  damping  as  well. 

It  should  be  noted  here  that  mode  3  is  the  fundamental  torsion  mode.  This  mode 
will  not  be  taken  into  consideration  for  the  design.  However,  it  will  be  tested  and 
compared  with  the  conventional  treatment. 

As  stated  in  Chapter  1 ,  a  sandwich  beam  offers  the  advantage  of  locating  the 
viscoelastic  material  in  the  region  of  the  beam  cross-section  having  the  maximum  shear 
stress  (the  neutral  surface).  Additionally,  it  is  possible  to  generate  large  bending 
deflections  with  relatively  small  forces.  In  a  conventional  treatment,  the  major  cause  of 
shear  strain  in  the  viscoelastic  layer  is  attributable  to  the  extension  of  the  substrate  at  the 
interface,  whereas  the  causes  of  shear  strain  in  a  sandwich  beam  are  both  the  longitudinal 
extension  of  the  constraining  layers  and  the  relative  motion  of  the  top  and  bottom 
interfaces  of  the  viscoelastic  layer  due  to  bending.  Considering  these  facts,  it  seems 
necessary  then  to  connect  this  sandwich  beam  in  such  a  way  that  these  purportedly  large 
bending  deflections  will  have  enough  physical  space  to  occur.  Therefore,  some  sort  of 
support  structure  will  be  needed  to  hold  the  sandwich  beam  away  from  the  truss. 

With  a  weight  constraint  on  the  design,  minimizing  the  weight  of  any  “non¬ 
damping”  parts  is  of  crucial  importance.  However,  in  order  to  generate  bending  in  a 
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sandwich  beam,  forces  from  the  truss  (which  will  be  felt  as  a  combination  of  axial  forces 
and  end  moments  by  the  sandwich  beam)  will  have  to  be  transmitted  through  the  support 
structure  to  the  sandwich  beam.  Therefore,  we  need  a  lightweight  support  structure  that 
can  resist  the  applied  forces  as  rigidly  as  possible,  for  any  deformation  in  the  suppoit 
structure  will  likely  cause  a  decrease  in  damping. 

The  next  question  is  one  of  how  to  generate  large  bending  deflections  in  a  sandwich 
beam.  If  the  sandwich  beam  needs  to  feel  end  moments,  then  we  must  look  to  optimize 
the  product  of  force  and  distance.  The  following  configuration  was  considered  to  be  the 
most  practical  for  the  given  problem. 

Consider  a  rigid  rod  positioned  vertically  with  a  pivot  on  the  bottom  (the  lower 
lever  of  Figure  5.2  illustrates  this  setup).  When  the  bottom  point  rotates  a  small  amount 
0,  all  points  on  the  lever  will  rotate,  and  all  points  except  the  pivot  will  displace  some 
amount  according  to  A  -  xd.  Now,  place  the  pivot  somewhere  on  the  truss  (like  the 
location  depicted  in  Figure  5.1).  Assuming  the  truss  will  vibrate  in  a  given  way 
regardless  of  what  is  connected  to  it,  we  can  theoretically  generate  as  much  displacement 
A  as  we  want.  Unfortunately,  several  factors  limit  the  practical  maximum  length  h.  If  a 
resisting  force  P  is  applied  at  the  tip  of  the  lever,  the  angle  6  will  be  reduced  (for  a  static 
case,  this  reduction  in  0  can  be  found  by  summing  moments  about  the  pivot). 
Furthermore,  the  resisting  force  P  will  cause  some  bending  deflection  in  the  lever  to 
further  reduce  the  displacement.  Finally,  the  longer  the  lever,  the  heftier  it  must  be  to 
rigidly  transmit  the  desired  displacements  to  the  sandwich  beam,  and  in  consequence  less 
weight  is  available  for  the  damper.  Assuming  some  optimum  length  can  be  found, 
consider  how  this  ‘^displacement-generating”  lever  can  be  connected  to  the  truss. 


sc 

K 


Figure  5.2  Two  Lever  Concept 
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The  cross-members  of  the  truss  provide  an  adequate  securing  location  for  the  levers. 
However,  an  increase  in  distance  between  two  points  does  not  solve  the  problem  of  how 
to  generate  bending.  Consequently,  a  second  set  of  levers  was  proposed  for  placement  at 
the  lever  tips  (refer  again  to  figure  5.2).  The  additional  lever  is  then  rigidly  attached 
perpendicular  to  the  sandwich  beam  end.  Figure  53  below  illustrates  this  configuration. 

As  point  o  (connected  to  the  truss)  rotates  by  an  angle  Q,  point  v  displaces  some 
distance  A  =  hd.  This  distance  is  approximately  equal  to  (for  small  ^).  By  choosing 
the  ratio  of  /1/5,  it  is  possible  to  control  the  ratio  of  Care  must  be  taken  to  avoid  a 
configuration  where  ^  becomes  too  large  (it  will  be  shown  later  that  the  limits  of  the 
truss's  deflection  will  prevent  ^  from  becoming  too  large). 

In  Figure  53,  the  connection  of  levers  h  and  5  is  a  pin  joint,  and  the  short  levers  and 
sandwich  beam  maintain  a  right  angle  orientation  with  respect  to  each  other  throughout 
the  vibration  cycle .  In  addition,  assume  the  pin  joint  is  frictionless. 

Placing  this  setup  on  the  truss,  we  have  the  following  situation. 


In  addition  to  bending,  it  is  possible  to  generate  direct  shear  strains  in  the 
viscoelastic  layer.  Consider  the  following  diagram  of  a  sandwich  beam. 


Figure  5.4  Exploded  View  of  Sandwich  Beam 

(VI«oo«lMtlc  layw  tMcknew  greatly  exaggaratad) 

If  the  left  lever  is  connected  to  the  bottom  sandwich  plate  and  the  right  lever  connected  to 
the  top  plate,  the  load  P  will  become  a  shear  load  that  the  viscoelastic  layer  must  resist 
for  equilibrium. 

By  equipollency  and  assuming  small  ^  the  direct  loads  P  produce  corresponding 
moments  (equal  to  the  product  of  P  and  .v)  at  the  end  of  the  sandwich  beam.  Therefore, 
the  sandwich  beam  is  in  a  state  of  combined  loading  (axial  and  bending).  We  will 
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assume  that  the  displacement  equation  for  any  inchvidual  layer  of  the  sandwich  beam  is 
the  superposition  of  displacements  due  to  each  separate  loading.  This  assumption  will  be 
shown  to  be  valid  (for  small  deflections)  in  spite  of  beam-column  effects  that  occur 
during  the  compressive  half-cycle  of  oscillation.  Since  the  levers  will  move  back  and 
forth  when  the  truss  vibrates,  both  loadings  will  be  sinusoidal  in  time  according  to 


where 


P{t)  «  P^sixKOj-t  and  M(t)  -  Pj  sino),-/ 


(5.1) 


P„  -  the  maximum  load  P  due  to  maximiun  truss  deflection  for  a  given  cycle 

(Oj.  -  the  mode  1  natural  frequency  of  the  truss 

We  will  assume  the  natural  frequency  of  the  truss  is  minimally  affected  by  the 
addition  of  the  damping  device.  For  this  reason,  changes  in  the  viscoelastic  material 
properties  are  considered  to  be  negligible.  Although  the  direct  load  is  time-dependent, 
the  natural  frequencies  of  longitudinal  vibration  can  be  safely  assumed  to  be  much  greater 
than  any  truss  frequency  of  interest,  so  a  static  loading  analysis  is  in  order.  However,  the 
natural  frequencies  of  transverse  vibration  cannot  generally  be  assumed  to  be  much 
greater  than  <ur>  so  a  dynamic  analysis  is  required. 

The  sequence  of  the  sandwich  beam  analysis  will  be  as  follows: 

-  develop  the  displacement  equations  for  a  three  layer  sandwich  beam  in  direct 

shear 

-  determine  the  effect  of  the  lever  geometry  on  the  dynamic  bending  analysis 

-  statically  analyze  an  eccentrically  loaded  beam-column 

-  discuss  asynunetric  phenomena 

-  develop  the  homogeneous  differential  equation  of  motion  for  a  three-layer 

sandwich  beam  (most  of  this  development  is  in  appendix  D) 

-  establish  boundary  conditions  for  the  solution  of  the  equation  of  motion 

-  solve  the  differential  equation  for  the  displacement  w{x,t) 

-  determine  shear  strains  and  system  loss  factors. 
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Displacement  Equations  for  a  Three  Layer  Sandwich  Beam  In  Direct  Shear 


Consider  a  slice  dx  of  the  beam  of  unit  width  illustrated  in  figure  5.5.  Assume  there 
is  no  shearing  strain  in  the  elastic  members  and  that  the  viscoelastic  layer  has  no  axial 
strain.  Assume  also  that  the  shear  strain  is  distributed  uniformly  through  the 
viscoelastic  thickness  and  the  axial  stress  is  uniformly  distributed  through  the  thickness 


dx 


P*—  x'-*  •  '  '  •  i  i  ^  ‘  ‘  ‘  •  >■! 

(thickiMMM  exaggeratad  lor  etarily) 

a,  - - 


Figure  5.5  Typical  Sandwich  Beam  and  a  Differential  Slice  dx 


of  the  constraining  layers.  We  will  prescribe  the  constraining  layers  to  be  made  of  the 
same  material  and  have  the  same  thickness.  Finally  assume  that  Young's  modulus  of  the 
elastic  layers  is  much  greater  than  the  shear  modulus  of  the  viscoelastic. 

Looking  at  the  differential  slice  depicted  in  figure  5.5,  we  can  sum  the  forces  in  the 
upper  and  lower  constraining  layers  to  obtain 

r,  +  T<&  -  0  and  tj-^^dc-rdc  -0  (5.2) 

'  dx  ^  dx 


respectively.  The  shear  strain  in  the  viscoelastic  layer  y  (assuming  y  is  small)  is 


y 


(53) 


Using  the  stress-strain  constitutive  relations. 


T  -  Gy  -G 


and 


dx 


Substitution  of  (5.4)  into  (5.2)  gives 


(5.4) 
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t. 


0  and 


o2  ^ 


-0 


(5^ 


Rearranging  (5.5),  we  obtain 


£!fi 

dbc* 


G 


(“l-“3)-0 


and 


G 

^3^2 


(«l-“3)-0 


(5.6) 


For  convenience,  let 


and,  since  r3  =  /i,let 


__G _ G _ g 

’  2 


(5.7) 


(5.8) 


The  set  of  two  homogeneous,  coupled  differential  equations  is  then 

“3 

Subtracting  the  second  equation  from  the  first  yields 

(«i  -"3  )-«{“i-“3)-0 

Let  A  =  Ml  -  M3  and  A"  =  mi"  -  M3"  to  simplify  (5.10)  to 


A"-gA-0 


The  solution  to  (5. 1 1 )  is 


A(jc)  -  A  sinhAjr  +  B  coshAr 


(5.9) 


(5.10) 


(5.11) 


(5.12) 
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where  A  aud  S  are  constants  and 


A-V?  (5.13) 

The  boundary  conditions  for  this  problem  are 

A(0)  -  A(L)  -  e  (5.14) 

where  e  is  some  difference  in  the  end  displacements  with  respect  to  each  other.  This 
quantity  will  be  determined  later.  Solving  the  left  boundary  condition, 

A(0)-fi-e  =>  A(jc)  -  A  sinhAjc -t- e  cosh  Ac  (5.15) 


Using  the  right  boundary  condition. 


A(L)*-  AsinhAX.-f  e  coshAL  >  e 


A-  e 


(1  -coshAL) 


sinhAL 


-etanh 


AL 


(5.16) 


The  expression  for  A  becomes 

A(jc)  -  f^-tanh  -^sinhAx  +  coshAxj  (5.17) 

Since  A  =  «i  -  «3, 

«,''-^A-0  (5.18) 

2 

We  can  arbitrarily  set  mi(0)  to  zero.  This  defines  the  boundary  conditions  for  the  solution 
ofu]  as 

«,(0)-0  and  a,'(L)-0  (5.19) 

Substituting  (5.17)  into  (5.18)  gives 
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^  2 


-tanh  — sinh  Ar  +  cosh  Ax 
2 


(5.20) 


Integration  of  Uiis  expression  yields 


^  .ii 

^  2A 


-  tanh~cosh  Ax  +  sinh  Ax 
2 


+  C, 


Using  the  right  boundary  condition, 


MO-o  -  f 


AL 


-tanh — cosh  AL  +  sinh  AL 
2 

XL 

-tanh — coshAL-i-sinhAf. 
2 


+  C,  -0 


Therefore, 


(5.21) 


(5.22) 


^  2A 


-  tanh— cosh  Ax  +  sinh  Ax 

_gg 

AL 

-  tanh — cosh  XL  +  sinh  AL 

2  J 

2k 

L  2 

(5.23) 


Integrating  once  again,  we  obtain 


XL 

-  tanh — sinh  Ax  +  cosh  Ax 
2 


2A 


-  tanh — cosh  AL  +  sinh  AL 
2 


x  +  C, 
(5.24) 


The  left  boundary  condition  gives 

u,(0)-0 


-  2A- 


(5.25) 


Consequently,  the  expression  for  mi  is 


M,(X) 


2A- 


-  tanh  —sinh  Ax  +  cosh  Ax 
2 


£g 

2A 


AL 

-tanh — cosh  AL  + sinh  AL 
2 


X- 


gg 

2A' 

(5.26) 


To  find  M3,  we  use  the  relation  M3  =  mi  -  A  to  get 
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U^ix)  ~u^{x)~  ^ 


-tanh  —sinh  Ax  +  cosh  Ax 

-11 

-tanh— cosh  AL  +  sinh  AL 

[2k^  1 

2  J 

‘’2A 

2 

ge 

{5.rj) 


We  now  need  to  find  the  unknown  quantity  e  in  terms  of  known  quantities. 
Differentiating  U3(x)  gives 


Mj  (jc)-  -As^-tanh^coshAjc  +  sinhAxj 
At  this  point,  we  prescribe  right  boundary  condition  usXL)  to  be 


'  P 

u,  (L) - 

^  AE 


where  A  =  /i  for  a  unit  width. 

Substitution  into  i43'(Jt)  gives  the  value  of  e  as 


AAE 


-tanh  ~coshAL  +  sinhAL 
2 


Substitution  of  e  into  the  expressions  for  ui(x:)  and  u^ix)  gives 

JgL  . 

I  ■ 


1- 


u^(x) 


-  tanh  sinh  Jgx  +  cosh  Jgx  I 


UgAE 


u^{x) 


UgAE] 


-tanh^^  cos  •f  sin  hJgL 

I  -tanh  ^^^smh^x  +  cosh-Jgx  \ 
-tanh  ^^^cosh^L  +  sinh^L 


+  1 


Px 

2AE 


Px 

2AE 


(5.28) 


(5.29) 


(530) 


(531) 


(532) 


What  we  are  really  interested  in,  however,  is  the  magnitude  of  u  i  -  U3.  This 
expression  is 
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Ic  Layer  Shear  Strain  (rad) 


(533) 


u^ix)-U3ix) 


p 

-  tanh^^^sinh.^x  +  co^Jgx 

-tanhJ^^cosh.^L  -t-sinh.^L 

For  the  given  set  of  parameters,  the  difference  of  ui  -  is  illustrated  in  figure  5.6. 


P=  2  lb.  £=  1 E7  psi,  G=  50  psi,  fi  =  0.032r.  ^  =  O.OOy.  L=  14*.  and  =  0.173  in* 


Figure  5.6  Shear  Strain  of  Viscoelastic  Layer  due  to  Axial  Loading 


A  quick  substitution  into  equations  (5.9)  verifies  that  these  expressions  are  indeed 
solutions.  Equations  (5.31-2)  represent  the  midpoints  of  the  constraining  layers  for  all 
time  (assuming  a  static  analysis  is  appropriate).  From  (5.3),  we  know  the  shear  strain  is 
the  quotient  of  U3  -  u  i  and  t2.  This  value  is  the  average  shear  strain  through  the 
thickness.  Intuitively,  this  concept  is  contradictory  since  we  know  that,  at  either  end,  one 
of  the  constraining  layers  has  maximum  stress  while  the  other  has  zero  stress.  This  is 
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considered  to  be  a  local  end  effect  that  dies  out  a  very  short  distance  into  the  beam 
according  to  St.  Venant's  principle.  Loosely  speaking,  the  shear  distribution  through  the 
thickness  will  self-equilibrate. 


Effect  of  Lever  Geometry  on  Bending  Analysis 

The  determination  of  the  lever  lengths  h  and  s  is  complicated  by  the  presence  of  an 
asymmetric  bending  phenomenon  between  the  tensile  and  compressive  halves  of  a 
vibration  cycle.  Gere  and  Timoshenko  [19:  567-8]  discussed  the  effect  of  eccentric 
loading  on  the  displacement  of  a  homogeneous  beam-column.  Their  development  for  the 
compressive  half  of  the  cycle  is  presented  here. 

Consider  an  eccentrically  loaded  beam-column  like  the  one  in  figure  5.7. 


Figure  5.7  Eccentrically  Loaded  Beam-Column 


Assume  the  distance  s  remains  constant  in  the  c-direction  and  that  the  levers  are  rigid. 
The  moment  at  any  location  x  along  the  beam  is 

M~P{s  +  w)  (534) 


From  strength  of  materials  theory. 


EIw"  -  -M  -  -P(s  +  w) 
..  P  -P 


w  + — w 
El 


El 


with  boundary  conditions  w(0)  -  h<L)  -  0.  If  we  let 


El 


(535) 


(536) 


the  homogeneous  solution  to  (5.35)  is 
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w{x)^  -  A  cos  kx+  B  sin  ^ 


(537) 


To  find  the  particular  solution,  assume  ~Cs+  D,  whose  derivatives  are 


w.  -  C  and  w.  -  0 

p  p 


(538) 


Substitution  of  this  form  into  the  differential  equation  yields  the  algebraic  expression 


P  P 

— (Cs+D)--- -s 
El  '  El 


C--1  ,  D-0 


(539) 


Therefore, 


(5.40) 


Using  the  boundary  conditions  to  find  the  constants. 


h'(0)-A-j»0  =>  A-j 
vv(L)-s  coskL-¥B  sinkL-smO 

(5.41) 

„  (l -cos  kL\  kL 

=>  — : -  -Stan  — 

[  sin  kL  )  2 

(5.42) 

The  expression  for  w(Jc)  is  then 

>v(x)  -  s^tan  ^sin  kx  +  cos  fcc  -  ij 

(5.43) 

The  maximum  deflection  occurs  at  x  =  ^2 .  where 

"’(f) 

(5.44) 

If  P  was  uncontrolled  and  allowed  to  grow  without  bound,  w(x)  would  grow  to  infinity  as 
*^/2  approaches  */2.  Fortunately,  the  distance  over  which  P  acts  is  controlled  by  the 
limited  deflections  of  the  truss.  Therefore,  we  need  not  be  concerned  about  exceeding  the 
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yield  stress  of  the  constraining  layers  at  jc  =  ^/2  and  about  the  effect  of  this  highly 
nonlinear  phenomenon  on  the  assumed  linear  analysis.  Therefore,  in  choosing  h  and  s, 
the  check  must  be  made  that  ^l2  is  not  close  to  */2. 

Consider  the  other  half  of  the  cycle,  where  P  is  pulling  the  lever  tips  apart.  The 
moment  M  becomes  M  =  P{s-W)  since  the  beam  is  deflecting  toward  the  load  center. 
The  differential  equation  then  becomes 


w" - w  «  -P{s  -  w) 

El 


The  solution  to  this  equation  is 


w(jc)»  s 


I  • 

tanh — SI 

V  2 


sinhkt  -coshfcc  +  l 


The  maximum  displacement,  at  jc  =  A/2 ,  is  given  by 


(5.45) 


(5.46) 


(5.47) 


The  value  of  1-sech  *^/2  approaches  I  zs  k  gets  large.  This  phenomenon  is  a  problem 
because  stresses  will  grow  high  at  the  comers  formed  by  the  short  lever  and  the  sandwich 
beam.  Figure  5.8  illustrates  the  difference  between  the  secant  and  hyperbolic  secant 
functions  as  their  arguments  become  large.  For  comparison,  the  strength  of  materials 
solution  is  included.  This  solution,  based  on  a  (1  -  cos  0)  relationship,  which  is  what  a 
linear  finite  element  solver  uses.  Conveniently,  over  a  significant  range  of  values  of  ^^/2. 
the  ordinate  of  the  strength  of  materials  solution  is  close  to  the  midpoint  between  the  two 
other  functions.  This  fact  essentially  allows  us  to  greatly  simplify  the  analysis  to  a 
symmetric  loading  throughout  the  vibration  cycle  if  we  are  willing  to  accept  a  small  error. 
Since  the  strain  energy  is  proportional  to  the  square  of  the  curvature  (which  is 
proportional  to  the  moment),  this  error  will  be  squared,  which  could  produce  significant 
error  at  extreme  deflections  in  the  cycle.  If  the  value  of  is  monitored  throughout  the 
design,  we  can  safely  neglect  the  effects  of  an  asymmetric  bending  cycle. 

The  equations  of  the  curves  in  this  illustration  are  the  following: 


wi^2)  For  Tension  =  5(sec(*^/2)-l) 


(5.48) 
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w(L/2) 


Comivession  =s(I-secIi(*^/2)) 
Circular  Curvature  =  s{  l-cos(*^/2)) 
5  =  1  for  all  curves 


(5.49) 

(5.50) 


kL/2 

Figure  5.8  Maximum  Tensile  and  Compressive  Cyde  Displacements  [w{^^] 


The  Differential  Equation  o/Trattsverse  Displacement  for  a  Three-Layer  Sandwich  Beam 

Mead  and  Markus  [8]  used  DiTaranto’s  derivation  of  the  equation  of  longitudinal 
displacement  of  a  three-layer  sandwich  beam  [7]  as  a  basis  for  developing  the  equation  of 
transverse  displacement  They  converted  the  longitudinal  variable  ^  to  the  transverse 
variable  w  and  arrived  at  a  sixth-order  differential  equation  similar  to  DiTaranto's.  Their 
development  is  given  in  Appendix  D,  but  the  missing  steps  from  their  analysis  are 
rigorously  developed.  The  resulting  equation,  for  no  external  forces  applied  in  the 
interior,  is 
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0 


(531) 


For  this  equation. 


8 


(pOfd^w 


2G 

Et,t,  ’ 


2D 


,  and  D  >■ 


6(1 -v^) 


(5.52) 


Boundary  Conditions  for  the  Sandwich  Beam 

Since  the  equation  of  motion  is  of  order  six,  we  need  six  boundary  conditions  to 
find  w{x,t).  Four  boundaiy  conditions  (two  displacement  and  two  force)  come  from  the 
transversely  vibrating  composite  beam.  The  other  two  (displacement)  boundary 
conditions  come  from  the  restraint  of  longitudinal  motion  of  the  constraining  layers.  The 
displacement  boundary  conditions  are  all  homogeneous,  but  the  forced  boundary 
conditions  are  nonzero  and  time-dependent.  Solutions  to  problems  of  this  type  are  not 
generally  obtainable  by  using  the  classical  method  of  separating  the  spatial  and  temporal 
variables.  Therefore,  in  order  to  solve  this  equation,  we  must  convert  the 
nonhomogeneous  boundary  conditions  to  homogeneous  boundary  conditions  through  a 
transformation  of  variables.  The  effect  of  this  transformation  is  that  the  homogeneous 
equation  of  motion  becomes  nonhomogeneous.  The  resulting  forced  vibration  problem 
can  then  be  solved  by  modal  expansion.  Mindlin  and  Goodman  [20]  developed  a 
procedure  for  solving  these  types  of  problems  by  using  what  they  call  “elementary 
mathematical  techniques”  as  opposed  to  integral  transform  methods  that  use  the  Laplace 
and  Fourier  transforms. 

Their  solution  method  involves  the  construction  of  a  comparison  function  that 
effectively  homogenizes  the  time-dependent  boundary  conditions.  Then,  modal  analysis 
is  employed  to  find  the  forced  response  solution  to  the  nonhomogeneous  equation  of 
motion. 

Mead  and  Markus  [8]  discussed  the  possibility  of  entirely  real  solutions  to  the 
sixth-order  differential  equation.  They  reported  that  if  a  beam  is  simply  supported  at  its 
ends  and  the  faceplate  displacements  are  unrestrained  at  both  ends,  then 

w.  -  Af  sin  A,jc  (533) 


is  a  possible  solution. 
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Mead  and  Markus  [8]  developed  expressions  for  a  number  of  different  boundary 
conditions.  For  this  analysis  it  is  assumed  that  the  rotary  inertia  of  the  small  levers  is 
zero.  This  is  clearly  unacceptable  at  high  frequencies,  but  at  low  frequencies  the 
associated  error  will  be  small.  In  addition,  axial  loading  will  be  neglected  as  it  has 
already  been  accounted  for.  Assume  the  changes  in  the  lever  angles  cause  a  negligible 
transverse  displacement.  This  assumption  gives  two  displacement  boundary  conditions: 

<5-54) 


Two  more  boundary  conditions  come  from  the  free  edges  of  the  plates: 


dUj 

^X 

X.L  ^ 

0 


Symmetry  gives 


dUj  _  ^u^ 

dx  dx 


=  0 


(5.55) 


We  can  equate  these  boundary  conditions  in  u  to  boundary  conditions  in  w  by  using 
equation  (D.23)  (found  in  Appendix  D). 


and 


du, 

D 

d*w  yd‘w  (pt)  d^w 

dx 

-6  + 

dx*  dx-  ~  D  dt\ 

(5.56) 


(pO  d‘Y/ 
D  dt- 


(5.57) 


To  find  the  moment  boundary  condition,  consider  the  summation  of  moments  M/  of  each 
individual  layer. 


d'w 

D. — and  M, 
dx~ 


(5.58) 
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M2  results  from  the  forces  Pi  and  Pi  which  are  opposed  and  offset  by  distance 


fj  1+il- 


W 


‘2// 


-P. 


Ti 


\\\ 


// 


( 


'■‘V 

\  *5 


~Pt 


y£J!!L- (P^) 

~  5*“  ^  dx^~  D  dt^ 


D 

g 


d*w  (pt)  d^w 


(5J9) 


Summing  A#/  gives  the  moment  as 


M.® 

g 


obc" 


obc* 


D  *■ 


(5.60) 


The  boundary  conditions  are  then 


\f(0)  • -M^sincOjt  and  M{L)  -  M,, sin (Ojt 


(5.61) 


where  Mb  is  the  maximum  external  moment  imposed  by  the  vibrating  truss  at  the  ends  for 
any  given  cycle.  One  important  simplifying  assumption  is  that  the  moments  prescribed  at 
each  end  are  equal  in  magnitude  and  opposite  in  sense  for  all  time. 

Mindlin  and  Goodman  used  a  convenient  operator  notation  that  is  suitable  here. 
Define  the  boundary  condition  operators  B,  to  be 


^-^4-1 


(5.62) 


where 


-- 

g 


d*  ..  ipt)  2 

dx  dx  D 


(5.63) 
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and 


-01 


2 


-^6 


(pt)  j 

■TT  +  «^“rT+  n  " 
dx  dx  D 


(5.64) 


The  boundary  conditions  can  be  expressed  in  summation  notation  as 

BHO,/)] -/,(/)  .  i-1,  2.  3 
£5[»v(L./)]-/,(0  ,  i-4v  5^  6 


In  unabridged  form, 

fl;[vi<a/)]-0 

fl3[vt<aO]-  0 


fl4N^0]-0 
B5[w(L,/)]-  M^sino>7-r 
£t[w(L.O]-0  (565) 


where  o>r  =  the  frequency  of  vibration  of  the  external  force.  2^ro  initial  conditions  are 
assumed. 


Solution  of  the  Differential  Equation 

Since  the  boundary  conditions  are  non-homogeneous,  we  will  transform  the 
variable  w  into  v  and  force  the  boundary  conditions  in  v  to  be  homogeneous.  To  do  this 
transformation,  we  assume  a  solution  of  the  form 

w{x,  t)  -  v{x,t)  +  2  gfx)fit)  (5.66) 

I 


We  must  also  transform  the  boundary  conditions  into  v.  Substituting  (5.66)  into  (5.65), 


B,[w(0,r)]  -  fl,[v(0,r)]+  V  //r)fll[g,(x)] 

pi 

ft[w(0.0]  -  fl3[KO,0]+  ^/;(Oft[g/x)] 
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J‘t 

B,[w(L,0]  -  Bj;v(L,0]+2^j(t)B,[gj(x)] 

j-i 

ajw{i.<)]  -  2  w] 

>-• 

>-• 

(5.67) 

We  know  from  (5.65)  that  fi  =  -f5  =  -Mb  sin  mit  and  f\=h-h  =f6  =  0*  Since  the 
purpose  of  the  variable  transformation  is  to  homogenize  the  boundary  conditions  in  v,  we 
must  choose  gj^x)  so  that  Bi[v(x,0]  =  0.  Therefore,  we  rearrange  (5.67)  to  get 

fll[K0,()].  fliH<o./)]-{A(/)i([&(0)]+/,(/)a;[fe(0)]}  -o 
i^[  Ka  0]  -  a,[»>(o,  t)]  -  {/,(»)^(o)] + /,(()4*5(0)]}  -  0 
«,[Ka/)]-  fl,[>*<ft<)]-{A«)^[&(0)]+/,{()ii(&(0)}}  .0 

B4[><t./)]  -  B.[H<i.()]-{/2«)B4[&(0]+  /s«)B4[fo(t)S  -0 

Bjv<£,l)]  -  BiH<£.<)]-{A(»)8,[fe(i)]+ /,(()B.[*,(i)]}  -0 

B[  KAO]  -  BJivC  AO]  -  {/jWBiUiC  A)]+ /,(04[*4(i)]}  -  0 

(5.68) 

In  order  to  satisfy  these  relationships. 

«{«/<»]-«,  .  <  -l.  z  3 

B[*;(A)].«,  .  <-4  5  6 

(5.69) 

The  initial  conditions  expressed  in  v  are 

v(.*;0)-0-2/i(0)g,(A:)-0 

i-\ 

v(jtO)  =  0  -  2//(0)ft(*') »  0iTMb[g2ix)  -  «5('^)) 

i-1 

(5.70) 

since 

(OjM^  cos(w^(0)  - 
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g2(x)  and  g5(x)  must  be  determined.  We  assume  a  polynomial  form  of  gi,  and  since 
at  least  one  of  the  operators  is  of  order  zero,  let 

gjix)  +  bjX  +  CjJC^  +  djX^  +  +  AjO:’  (5.71) 

To  solve  for  each  operator  to  construct  an  algebraic  equation  that  is  either 

equal  to  zero  or  unity  according  to  (5.69). 

^b2(0)]-0  =»  Oj-O  (5.72) 

^[ft(0)]  - 1  -  7  *  »’)5-  ^ 

^^^e2  +  2Z)(l+K)c2 -1 

g 

fi^[&(0)]  -  0  =>  ~  24e2  +  2gy£^  -  0  (5.74) 


£2  and  C2  are  now  determined  as 


(5.75) 


Next,  we  solve  the  boundary  conditions  at  L  by  considering  e2  and  C2  to  be 
constants. 

B,[g,(L)]-0 

Lb.  +  l^d.  +  Zf ifc,  -  -  Z?c,  -  L*e.  (5.76) 

Bs[griL)]^0 

+  |-6g(l  +  JOf-  -  d.  +  120L  -  20g(l  +  Y)V  *2 

-  |2g(l+K)  +  ^a>=L"  c,+  +12g(l+m"-24je,  (577) 

4[«2(^)]-0 
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(5.78) 


Inserting  (5.80)  back  into  the  g2  boundary  conditions  will  verify  that  they  are  indeed 
satisfied.  It  remains  to  find  gsix).  We  again  assume  the  polynomial  form 

gj(jc)  •a2+bgX  +  +  dgX^  +  e^x*  +  k^x^  (5.81) 

^[«s(0)]-0  =*‘  0,-0 

ftk5(0)]-0  =>  -2Ae,  +  2gi\  +  Y)c,~0  (5.82) 

fi3[«5(0)]-0  =>  -24€5  +  2grc5-0 


Clearly,  05  =  65  =  0.  The  boundary  conditions  at  L  are  the  following: 

«.[&(«] -0 

Lb,  +  ed,  +L%  -0 


(5.83) 


d£!l(o'‘L\b,  +  \-6g{l+  Y)L-^-^(o-L^]d.  +  fl20L-20g(l  +  K)Z?  - 

Id  J  L  D  J  I  D  J 


1^(0- L]bs  +  \-6gYL-^a)^L^  U  +  [l20L-  lOgYL^  -  0 

D  J  L  D  \  [  D  J 


(5.84) 

(5.85) 
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Solving  (S.83-SS)  simultaneously. 


eQL-lgYJl  ^  gYB-6 
360D  '  36DL 


and 


gY 

120DI 


(5.86) 


Therefore, 


€OL-7gY£^^gY£-6^,  gY  ^5 
360D  36DL  120DL 


(5.87) 


This  function  also  satisfies  the  boundary  conditions.  At  this  point,  g2  and  gs  will 
homogenize  the  boundary  conditions. 

The  next  step  is  to  solve  the  nonhomogeneous  equation  of  motion  in  v  with  the 
boundary  and  initial  conditions  expressed  in  v.  Meirovitch's  text  [21:  303-5]  contains  an 
excellent  summary  of  the  following  solution  to  the  nonhomogeneous  differential 
equation.  We  shall  express  the  differential  equation  in  his  operator  notation  as 

I[v(x,t)]  +  M(x)’^v(x,t)  -  Fix,t)  -  (0f[&(^)]+  /i(0Af(x)ft(x)]  (5.88) 


where 


— ^-gd+K)-:^ 


MO 


,  i-2 

,  /-5 

0  otherwise 


First,  the  homogeneous  eigenvalue  problem  in  v  is  solved.  Now  that  the  boundary 
conditions  are  homogeneous,  we  can  separate  the  spatial  and  temporal  variables  by 
assuming  a  solution  of  the  form 


v(j;0-2A,(x)r.(0  (5.89) 

ii^i 
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The  eigenvalue  problem  becomes 


Ilv(x)]  -  fo^l^vix)]  -  0  (5.90) 

Since  we  are  prescribing  a  moment  at  the  ends,  assume  the  mode  shapes 

v,(x)  -  V,  sin^  (5.91) 

Let  us  pause  and  consider  the  validity  of  this  assumption.  If  we  investigate  the  dynamic 
motion  of  the  fundamental  mode,  it  is,  in  fact,  asymmetric.  Imagine  pressing  down  on 
the  center  of  the  left  sandwich  beam  depicted  in  Figure  5.8.  This  half  of  the  vibration 
cycle  resembles  very  much  a  beam  simply  supported  at  its  tips.  However,  if  we  pull  up  in 
the  middle,  this  deformation  is  better  described  by  a  fixed-fixed  beam.  In  fact,  it  will  be 
even  stiffer  than  a  fixed-fixed  beam  because  as  deformation  causes  the  levers  to  move 
inward,  the  comers  will  begin  rotating  opposite  the  motion  of  the  beam  center.  Clearly, 
by  assuming  this  mode  shape,  we  are  describing  only  part  of  the  cycle  partly  accurately. 
If  we  are  concerned  with  determining  the  worst-case  displacement  scenario,  equation 
(5.91)  is  probably  the  most  conservative  full-cycle  mode  shape  to  assume.  However, 
depending  on  the  geometry,  damping  factors  may  or  may  not  be  adversely  affected  by  the 
limited  displacements  in  the  stiffer  half-cycle.  If  the  dynamic  part  of  the  solution 
becomes  the  dominant  term  (as  in  frequencies  near  sandwich  beam  resonances),  then  a 
more  accurate  investigation  need  be  made  with  regard  to  what  mode  shapes  to  assume. 


Figure  5.9  Deformed  Geometry  of  Beams  Pinned  at  Supports 
(comer  right  angles  maintBined) 


If  we  assume  (5.91),  then 

.  nxx  (  mY  .  nnx 


d*  .  mcc  /  rm\  *  .  mix 


47 


(5.92) 


.  ruK 
— rsm— — 
L 


.  ma 


Substitution  of  these  derivatives  into  (5.51)  yields 


This  expression  for  <0^  is  the  same  result  obtained  by  Torvik  [17]. 

Now  we  must  normalize  the  coefficients  Vn  with  respect  to  M{x)  of  the 
eigenfunctions  v,^x).  This  step  turns  out  to  be  necessary  in  order  for  the  units  of  the 
“dynamic”  part  of  the  solution  to  match  the  units  of  the  part  of  the  solution. 

Since  M(x)  is  a  differential  operator, 

(5.95) 


where 


.  rnnx  .  mtx 

and  v.-V.sin  — 

f  ,,  .  mmc  ,,  {pt){  (rm\'  .  rmx  . 

(-(t) 


dx  »l 


(5.96) 


For  m*n,  the  value  of  the  integral  is  zero.  For  m  =  n. 


(5.97) 


We  assumed  a  solution  of  the  form 
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(5.89) 


v(j;/)  -  |jv,(x)r,(/) 

When  this  form  is  substituted  into  the  non-homogeneous  equation  of  motion,  we  have 

|[r,(()4v.w]+  (5.98) 

Since  Vn  and  satisfy  the  eigenvalue  problem  ( 5.90),  we  can  rewrite  this  expression  as 

iJnCO  +  <»/r,(»)]A4v.(x)]  -  -|[/,«)4a(-«)]  +  /,(»)*4'’.W]]  (5.99) 

This  infinite  set  of  coupled  equations  can  be  uncoupled  by  using  orthogonality.  We 
multiply  each  side  by  v^ix)  and  integrate  from  0  to  L.  We  now  introduce  the 
orthogonalized  spatial  coefficients 

^8.  •  t^nix)^^g^{x)]dx  and  G*  -J^v,(x)L[ft(^)]<fc  (5.100) 

This  step  gives  the  uncoupled  equations 

f,(/)  +  a>;r,(/)-A^,(0  (5.101) 

where 

^<»)-|[g;/,(')-g,, /,{<)] 

The  non-zero  terms  G_, ,  G.. ,  G*  ,  and  G'  are 

lij  ^5  "2 
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G„  -  V.  <^{(-l)’[-2ci(i)  -(H  *(i2^  -"‘(:^)l 

.(20*,  -*4)[6t(i)’  -  j  .  .  24(i)’  -  12L^[  A)’] 

.«*,[t’(A)  .  ,aot(A)’  -20L»(i)’]].2c,(i)  -2(.2. 


(5.102) 


G.,  -  V  M(.l)-|(6<y,  -,.j(-z(i))  .(20*i  -i>(^)) 

+g/:  [ Zffil]  +  120l(— -20I^(i)^) 

^  \ruc/  \n3t}  \rui/  j 


(5.103) 


G;,  -  -V.s(Ur)[(-.)-(-24.,(A)-,20*,L(±)).2^,(^)]  (5.104, 


G„  -120(-l)-V.g(l  +  r)L 


a 


(5.105) 


The  length  of  these  terms  is  a  result  of  integrating  the  product  of  a  sine  function  and  a 
fifth-order  polynomial  (several  integrations  by  parts  are  needed). 

Nnit)  can  then  be  expanded  as 


A'.(»)  G.,A(»)+ c;,/,(o  .  c„7,(»)) 


(5.106) 


For  zero  initial  conditions,  the  temporal  solution  is 


T„it)  -  7“ ji A^,( T)sin t) dr 


(5.107) 


We  can  now  write  the  solution  for  v{x,t)  as 
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(5.108) 


This  is  the  relative  displacement  of  the  sandwich  beam  with  respect  to  a  translating, 
rotating  truss.  This  is  an  approximate  solution  because  the  polynomials  are  merely 
comparison  functions  and  because  the  mode  shape  assumption  made  on  the  first  term  is 
approximate.  A  cursory  look  at  this  solution  indicates  that  the  "forcing”  term  will  be 
equal  to  zero  for  all  of  the  boundary  conditions  since 


(5.110) 


Since  gi(x)  were  prescribed  to  satisfy  the  boundary  conditions  in  v,  they  will 
automatically  satisfy  the  boundary  conditions  in  w.  Additionally,  an  inspection  of  the 
comparison  functions  indicates  that  the  unit  is  indeed  a  displacement.  Strictly  speaking, 
the  unit  is  in-,  but  the  moment  is  divided  by  the  width  since  this  is  a  per  unit  width 
analysis.  Finally,  the  suspicious  lack  of  an  **^/£/-looking'’  term  is  explained  by  the 
presence  of  D  in  each  denominator.  Assuming  the  end  moments  are  equal,  the 
comparison  function  reduces  to 
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The  following  graph  shows  the  comparison  function  versus  a  sine  wave  (normalized  with 
respect  to  each  other). 


Figure  5.10  Displacement  of  the  Comparison  Function  Terms  as  a  Function 

of  ^4  Normalized  to  a  Sine  Curve 


An  investigation  of  the  convolution  integral  term  gives 


a) 


SB 


-fg  sinctfrT  sincUj^g( 


1 

( 

(Oj  sincuss/ 

oi^gsincu^r 

<»SB 

(5.111) 


The  evaluation  of  the  convolution  integral  reveals  that  the  resonance  phenomenon  is 
accounted  for.  The  oj-p-ios^  term  in  the  denominator  suggests  that  the  amplitude  of 
vibration  grows  without  bound  as  cur  approaches  (jusb-  will  not  be  the  case  since  the 
sandwich  beam  is  damped.  However,  an  investigation  of  ojsb  is  warranted  for  likely 
configurations.  The  natural  frequencies  are  calculated  from  the  equation 
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(5.112) 


f  IK. 

ipAif 

where  p  for  a  simply  supported  beam  in  mode  \  =  jt  With  a  fundamental  truss  frequency 
on  the  order  of  unity  and  a  sandwich  beam  fundamental  frequency  on  the  order  of  35  Hz, 
the  forcing  term  is  likely  to  be  small  for  this  configuration  in  mode  1  vibration.  For  other 
configurations  or,  even  for  higher  modes  of  this  configuration,  this  term  will  certainly  be 
significant. 

The  forcing  term,  evaluated  for  mode  1  vibration  (n  =  1)  is 

V,  sin— ji  A^,(T)sinmsB^  (t  -  r)  dr 
L 


(5.113) 

Now  that  we  know  the  behavior  of  the  sandwich  beams  for  the  improved  treatment 
under  sinusoidal  loading,  we  must  determine  the  (Ustribution  of  strain  energy  so  that  the 
modal  loss  factors  of  the  system  can  be  determined. 

Cakulation  of  Strain  Energy 

In  order  to  obtain  the  loss  factor  rf,  the  energy  dissipated  per  cycle  Ds  must  be 
found.  The  quantity  Ds  per  unit  width  per  unit  length  is 

D,-nG%|y|-  (5.114) 

Appendix  E  contains  the  development  of  Ds.  To  find  Ds,  y  must  first  be  determined. 
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Earlier  in  this  chapter,  we  assumed  the  total  shear  strain  to  be  the  superposition  of 
shear  strain  due  to  axial  loading  and  shear  strain  due  to  bending.  Using  equation  (533), 
the  shear  strain  due  to  the  axial  load  P  is 


1 

1 

ar 

1 

-tanh  J^^sinh  Jgoc  +  coshJgx 
2 

/j  yfgAEt2 

-  tanh  sinh  +  cosh 

2 

"Sin 

(5.115) 


We  can  obtain  Yb  from  the  development  of  the  sixth-order  differential  equation  of  motion 
for  transverse  displacement  (Appendix  D).  The  equation  for  the  shear  strain  due  to 
bending  is 


r. 


dx 


(5.116) 


We  also  know  that 


£!ffL 

dx^ 


(5.117) 


In  addition,  we  determined  the  quantities  du^/dx  and  du^fdx.  Since  we  know  the  first 
derivatives  of  w,-,  their  second  derivatives  can  be  determined.  We  can  rearrange  (5. 1 17)  to 
isolate  the  part  that  equals  yt.  The  other  side  of  the  equation,  which  is  in  terms  of  known 
quantities,  must  also  equal  yi,. 


£r,  ^'u^  d'u^ 
2G  ^  dx^  ^ 


After  differentiating  ^u^|dx  and  (^3 is  found  to  be 


(5.118) 
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1  ptna 

— 7 - r  yC,c<» - i— - X 

/,?2(l+/|/^2)j  L  D  L  <Ogg^ 


2G[gEt,t^{lA-tJt,))[i^i  L 

~*®S«r  0>r  S*®  ^SB.^  +  ^SBn 


+  (f”^V  rrf”^V]  ®^rSinQ>5a,f-a)^a.sina>r/ 

[(V  L  j  ^  \  /  j  [  0)/  -  o>j^2 

+AffcSintt»j/|-gK(6a3  +24<i4x)-^cOj.2(a,  +2a^x  +  3a^^  +4a4r’)j| 


(5.119) 


where 


C  -  U. 


-[(6i-«t>)(i)  *6*z(i)’]^.[(l2£^  -*^)(^)  -«)(i)' 


l2L-gYn 


^  “~2D 


12£> 

gy 


Appendix  F  contains  the  unabridged  development  of  The  total  shear  strain  is  then 

y-y, -n  (5.120) 

There  are  six  separate  terms  in  x: 
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sinh-Jgx  ,  sinh-Jgx  ,  cos{nx/L)  ,  x  ,  x^,  and  x^ 


To  square  and  integrate  this  complex  quantity  over  the  length  would  involve  many 
separate  integrals.  Therefore,  a  numerical  technique  is  in  order.  Before  this,  however, 
some  simplification  can  be  done. 

An  investigation  of  the  expression  for  y  indicates  that  part  of  the  shear  strain  is 
dependent  on  sinwsat .  Since  the  amplitude  of  vibration  of  the  truss  is  transient,  the 

entire  solution  is,  in  fact,  transient.  However,  due  to  the  inherently  high  damping  of  the 
sandwich  beam  component,the  terms  having  a  factor  of  sinniss/  will  become  small  much 
sooner  than  the  terms  dependent  on  sintu^-f .  The  sinni^iBr  terms  represent  the  transient 
response  of  a  sandwich  beam  due  to  an  initial  impulse  excitation.  Once  the  sandwich 
beam’s  individual  response  is  damped,  these  terms  will  be  a  small  contributor  to  system 
response.  In  effect,  after  the  initial  damping  of  the  sandwich  beam  occurs,  we  will  have, 
loosely  speaking,  a  quasi-static  vibration  condition  that  will  be  perpetuated  by  the  truss. 
Therefore,  if  we  eliminate  the  sin  -dependent  terms,  the  resulting  loss  factor  will  be 

conservative.  The  effect  on  the  solution  will  be  as  follows:  if  the  initial  amplitude  is 
unity,  the  total  system  will  damp  relatively  quickly  over  the  first  few  periods  of  (Ogg  and 
then  settle  to  the  slower  decay  of  the  sina>rf  terms.  By  eliminating  the  sincOa^  terms,  we 
will  begin  decaying  at  the  rate  of  the  sinru^-t  terms.  Therefore,  the  amplitude  of  response 
will  be  higher  for  all  time  t  than  if  those  terms  were  included.  Remember,  the  goal  of  the 
project  is  to  make  the  sinm^t  rate  of  decay  greater  than  the  rate  of  decay  for  the 
conventional  treatment.  The  following  figure  shows  decaying  oscillation  envelopes  that 
illustrate  this  concept 
The  equations  for  this  illustration  are 


improved 

sb-0.4e"^’‘ 
sum  -  e"  **  +  0.4e‘’ 
conv.».6ef^ 


(5.121) 


To  reiterate,  this  simplification  is  only  reasonable  because  (O^b  »  o>r  • 
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Arbitrary  Amplitude  Scale 


Figure  5. 1 1  Illustration  of  Effects  of  Neglecting  Highly  Damped  Terms 


Extracting  the  sinc^ssr  -  dependent  terms  from  yt  and  setting  n  =  1  leaves 


"  2G 
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A^sinctfrq 
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\"r-  -<Os~8 


(5.122) 


-gK(6i^  +  24ajjc)  +  20,^:  +  SojJC'  +  4*140:^  )| 


The  expression  for  yi  is  unchanged  since  the  analysis  was  assumed  to  be  quasi¬ 
static. 
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Calculation  of  The  Mode  1  Loss  Factor 


The  expression  for  the  loss  factor  of  the  system  is 


n- 


D 

2jtU 


_ _ 


(5.123) 


In  order  to  calculate  ri,  the  unknown  quantities  Mb  and  P  needed  to  be  determined 
because  they  impact  the  viscoelastic  strain  energy  and  damping  energy  terms.  The  finite 
element  analysis  program  I-DEAS'™  (a  product  of  Structural  Dynamics  Research 
Corporation)  was  used  to  construct  a  truss  model.  This  model  was  the  same  one  used  to 
determine  the  undamped  natural  frequencies  and  mode  shapes.  Mb  and  P  are  dependent 
on  the  geometry  and  the  loading,  so  a  parametric  study  was  conducted  to  determine  the 
dimensions  h  (long  lever  length),  s  (short  lever  length),  t  (thickness  tj),  and  the  cross- 
section  of  the  levers.  The  goal  of  the  study  was  to  minimize  the  percentage  of  total 
system  strain  energy  located  in  the  truss  elements.  The  study  was  basically  an  iterative 
process  whereby  different  configurations  of  h,  s,  t,  and  the  cross-section  were  input  to  I- 
DEAS.  The  loading  was  assumed  to  be  a  linearly  increasing  distributed  load  from  zero  at 
the  clamped  end  of  the  truss  to  a  maximum  at  the  tip,  as  the  following  diagram  illustrates. 


Figure  5.12  Depiction  of  Finite  Element  Static  Load 

Interestingly,  this  load  gives  almost  the  identical  mode  shape  of  mode  1.  This  load  was 
applied  by  way  of  point  loads  on  the  nodes  of  the  finite  element  model. 

The  static  solution  from  I-DEAS  is  almost  exactly  equal  to  the  strength  of  materials 
solution  for  a  similar  loading.  The  finite  element  program  output  data  consisted  of 
displacements,  strain  energies,  and  element  forces.  The  total  elastic  strain  energy 
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solution  for  the  undamped  truss  was  0.92  lb.  in.  from  I-DEAS,  while  the  strength  of 
materials  solution  was 


-  0.84  lb.  in.  for  Af  -  2L"  -  3L*  +  -xA 

Jo  2EI 

6V  L  1 

(5.124) 


where  qo  for  this  calculation  was  backed  out  by  using  the  I-DEAS  moment  of  35.135 
lb.  in.  This  result  was  considered  to  be  close  enough  to  verify  the  finite  element  model. 
The  element  forces  data  was  then  assumed  to  be  close  enough  to  use  in  the  theoretical 
expressions  for  D^,  lJs^  and  i)s-  Appendix  G  briefly  describes  the  procedure  for 
determining  the  dimensions  of  the  damping  device.  As  a  result  of  the  parametric  study  in 
Appendix  G,  the  chosen  configuration  reduced  the  total  elastic  strain  energy  in  the  system 
from  0.92  lb.  in.  to  0.63  lb.  in.  (based  on  the  finite  element  output  data). 

During  the  parametric  study,  it  was  discovered  that  a  wide  sandwich  beam  placed  in 
the  root  truss  bay  has  the  effect  of  transferring  a  large  portion  of  the  strain  energy  in  the 
bay  down  to  subsequent  bays.  Therefore,  a  two-bay  configuration  was  considered.  The 
penalty  for  adding  a  second  sandwich  beam  is  the  weight  of  an  extra  set  of  levers.  This 
added  weight  reduced  the  total  allowable  width  of  the  sandwich  beams.  A  few 
preliminary  estimate  configurations  indicated  diat,  not  only  did  the  total  percentage  of  the 
strain  energy  in  the  sandwich  beams  increase,  but  the  total  elastic  strain  energy  was 
significantly  decreased.  Reducing  the  total  elastic  strain  energy  is  desirable  because  the 
value  of  the  denominator  in  r)  decreases. 

The  winning  configuration,  constrained  by  the  readily  available  materials  and 
construction  capabilities  of  the  resources  on  Wright-Patterson  Air  Force  Base,  is  the 
following: 


long  lever  length  h*7in. 
short  lever  lengdi  s  -  1.25  in. 
sandwich  plate  thickness  -  0.(X32  in. 

cross-section  -0.75 in.  x  0.080  in. 
wide  sandwich  beam  width  -  54  in. 
slim  sandwich  beam  width  h^  -  5.4  in. 

Figures  6.2-9  show  the  damping  device. 
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The  thickness  of  the  viscoelastic  layer,  t2>  only  parameter  left  to  be 

determined.  t2  is  a  key  parameter  in  the  expression  of  damping  energy,  which  is 

A-  /«G"|yfdV' 

(5.125) 


Once  the  configuration  was  set,  a  curve  of  t2  vs.  integrated  shear  strain  could  be 
constructed.  One  big  problem  with  this  method  is  that  G*  is  a  complex  quantity,  which 
makes  calculation  cumbersome.  Therefore,  Matlab™,  a  commercial  matrix  manipulation 
computer  program  designed  to  handle  complex  algebra,  was  employed  for  the  purpose  of 
iterating  through  a  number  of  thicknesses  (2.  A  numerical  iterative  function  file  was 
constructed,  and  results  were  analyzed.  This  program  comprises  Appendix  H. 

The  resulting  values  of  energy  (hssipation  for  nominal  thicknesses  were  several 
times  higher  than  the  elastic  strain  energy  in  the  trass.  Because  the  viscoelastic  material 
is  so  compliant,  it  is  conceivable  that  this  amount  of  dissipation  is  possible,  but  large 
deflections  must  occur  before  these  strains  can  be  realized.  It  seems  that,  for  a  given 
thickness,  assuming  the  truss  deflection  is  unaffected  by  the  sandwich  beam  thickness, 
doubling  the  beam  displacement  will  roughly  quadruple  the  energy  dissipation.  There 
seems  to  be  no  significant  penalty  for  trying  to  get  the  largest  elastic  displacements  the 
material  will  allow. 

During  the  iteration  process,  it  was  discovered  that  a  single  term  dominates  the 
expression  for  y  (for  configurations  of  interest).  The  term  is 


where 


Et, 

2D  ^ 

2G 

gY(6a3+24a^x) 


gYL  -gY  ^  (lever  tip  load  PXshort  lever  length  s) 
\2D  '  24D  ’  ^  sandwich  beam  width  b 


This  term  is  found  in  the  comparison  function  part  of  y.  The  term  is  linear  and 
antisymmetric  about  the  center.  If  we  examine  the  x  =  0  end  of  the  sandwich  beam,  only 
the  ‘*6a3”  part  of  this  term  is  present.  For  the  given  load  magnitude  (P  =  2.42  lb.  for  a  4" 
tip  displacement),  the  shear  strain  reaches  a  minimum  value  of  about  1.05  radians  at  t2  = 
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Value  of  the  Dominant  Terni  of  the  Shear  Strain  (rad.) 


0.015”  and  then  increases  with  increasing  t2.  The  figure  below  illustrates  the  behavior  of 
the  term  with  changing  t2. 


0.001  0.01  0.1  1 

thickness  12  (In.) 


Figure  5.13  Damping  Trend  of  Dominant  Term  With  Varying  Thicknesses  ^ 

Intuitively,  this  supposed  behavior  beyond  t2  =  0.015"  must  be  in  error.  Cottle  [221 
discovered  the  theoretical  phenomenon  of  increasing  damping  out  past  the  region  where 
damping  begins  to  decrease,  as  did  the  author  while  investigating  Torvik's  five-layer 
beam  analysis.  According  to  the  equation  for  energy  dissipated  per  cycle  Ds,  an  increase 
in  volume  of  viscoelastic  material  causes  a  proportional  increase  in  damping,  but 
viscoelastic  layer  thickness  has  an  inverse  squared  relationship  with  damping  energy.  It 
was  presumed  that  this  theoretical  increase  in  shear  strain  with  increasing  thickness  could 
be  the  result  of  the  breakdown  of  one  or  more  of  the  assumptions  of  the  analysis.  It 
appears  this  phenomenon  generally  becomes  significant  as  the  viscoelastic  material 
occupies  a  significant  percentage  of  the  beam  cross-sectional  area.  Torvik's  more  general 
analysis  (Appendix  B)  does  not  make  the  simplification  of  =  (3,  so  the  shear  strain  term 
remains 
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y* 


rather  than  y„-(l  +  /,/0 


dw 

lx 


When  ti  » (/2  and  (3) ,  this  effect  will  not  be  significant  until  t2  gets  really  large  because 
ti  dominates  the  coefficient  of  ^/ax . 

The  assumption  whose  breakdown  would  probably  have  the  greatest  impact  is  the 
one  that  says  the  viscoelastic  layer  does  not  resist  direct  strains.  Clearly,  as  tz  becomes 
large  (on  the  order  of  lOr),  the  bending  stiffness  becomes  significant  with  respect  to 
£1/1^,  especially  for  higher  frequencies  where  the  elastic  shear  modulus  G2  is  high. 

Since  a  close  scrutiny  of  the  algebra  did  not  result  in  finding  a  math  error,  it  was 
decided  to  test  a  number  of  different  viscoelastic  thicknesses.  Even  if  the  magnitude  of 
the  experimental  results  was  not  in  line  with  the  theory,  it  was  hoped  that  a  trend  could  be 
observed  from  the  testing  of  three  thicknesses. 

According  to  the  Matlab  numerical  integration  file,  the  expected  values  for  the 
mode  1  loss  factors  for  various  thicknesses  were  the  following: 


0,  (b.in.) 

0^  Ob.in.) 

Ih+UsB 

LUf 

n 

^325" 

!5!35 

b.65 

4.fl4 

. 1.78 

— (5S5" 

.015" 

25.56 

9.42 

0.65 

6.78 

^5 

0.56 

.030" 

67.21 

24.83 

0.65 

17.83 

6.59 

0.58 

Table  5.1  Theoretical  Damping  Values  for  Improved  Treatment 
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VI.  ExperimerU  and  Results 


The  two  parts  of  this  project  -  theoretical  development  that  led  to  the  design  of  the 
damping  device,  and  laboratory  testing  -  were  closely  intertwined  during  all  phases  of 
this  effort.  For  example,  the  design  of  the  conventional  damping  treatment  could  not 
proceed  without  the  resonant  frequency  data  obtained  from  impacting  the  truss.  Much 
has  already  been  said  with  regard  to  the  apparatus  of  this  study.  Chapter  3  contains  a 
detailed  drawing  of  the  truss,  and  Appendix  G  describes  the  dimensions  of  the 
components  of  the  improved  damping  device.  In  addition,  the  method  of  obtaining  modal 
frequencies  and  damping  ratios  is  well  known,  and  a  detailed  overview  of  the  testing 
procedure  is  not  necessary  (Appendix  A  describes  how  the  truss  frequencies  and  loss 
factors  were  obtained).  Consequently,  most  of  the  laboratory  worie  is  included  in  this 
section  of  the  report.  First,  the  improved  damping  device  will  be  described  in  detail. 
Then  a  brief  explanation  of  the  procedure  for  how  the  data  was  collected  will  follow. 
Finally,  the  reduced  data  (modal  frequencies  and  loss  factors)  will  be  presented. 
Appendix  I  contains  the  tabulated  raw  data. 


Apparatus 

The  electronic  equipment  used  for  collecting  data  was  a  typical  setup  for  modal 
analysis.  The  schematic  diagram  below  includes  the  most  important  components. 

Both  the  force  transducer  and  the  accelerometer  were  calibrated  for  consistency 
over  the  operating  ranges  of  force  and  acceleration,  but  determining  calibration  constants 
was  not  necessaiy  since  the  output  data  was  a  transfer  function  in  the  frequency  dommn. 
The  calibration  constants  would  merely  shift  the  entire  curve  up  or  down;  it  would  not 
affect  the  shape. 

The  conventional  damping  treatment  consists  of  16  aluminum  strips  measuring 
l.(X)”  X  18.0"  X  0.032.”  The  damping  material,  0.(X)5”  thick  3M  ISD-112,  came  in  the 
form  of  a  large  roll  and  was  applied  to  one  side  of  each  aluminum  strip.  The  strips  cover 
both  sides  of  each  longitudinal  truss  member.  The  strips  on  the  tip  end  of  the  truss  had  to 
be  trimmed  to  about  16." 
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Ptezolronics 
Impact  Hammer  with 
KMer  Foioe  Tranaducar, 
model  9712A5 


KteHer  Plezahon  Coupler 
model  5124 


27VDC 
Power  Supply 


Figure  6.1.  Electronic 

Equipment  Setup 


Kisller8832A5 

Accototometar 


The  components  of  the  improved  damping  device  are  illustrated  in  the  diagrams 
below.  The  sandwich  plates  overlap  by  about  12.5". 


Figure  6.2  Close-up  View  of  Truss  Damping  Device 
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Shoft  Levers 


The  plates  were  made  separately  and 
are  affixad  Id  each  other  only  by  the 
visco^astlc  material  between  them. 


^  Viscoelastic  Material  between  Plates 

Figure  6.3  Closeup  View  of  2.3”  Sandwich  Beam 


Long  Levers 


Bearing  (5/16*  outer  dia., 
3/16*  inner  dia.) 


Lever  Cross 
Section 
.080"  X  3/4' 
Made  from 
sheet  stock 

Fastening  Plate 


Oblong  slots  drilled  out  for 
the  purpose  of  taking  up 
slack  during  assembly 


Welded  to 
fastening  plate 


Figure  6.6  Load  Distribution 
Plate  for  Short  Lever  -  SB 
Plate  Interface 


Figure  6.5  Detail  of  Long  Lever  and 
Truss  Fastening  Plate 
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JaB* 

- ^ - 

1.25* 

k 

R=5/16* 

1.0" 

Figures.?  Short  Lever 
(not  to  scale) 


Rgure  6.8  Long  Lever 
(not  to  scale) 


The  fabrication  of  the  sandwich  plates  required  a  slight  modification  from  the 
original  design.  The  design  called  for  a  welded  connection  between  the  load  distribution 
plates  and  the  sandwich  plates.  Because  the  materials  are  so  thin  (0.032”),  warping  of  the 
sandwich  plate  would  surely  have  resulted  from  welding.  Consequently,  the  short 
lever/distribution  plate  component  had  to  be  riveted  to  the  sandwich  plate.  Figure  6.9 
below  shows  the  rivet  patterns  for  the  23"  and  5.4"  plates. 
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5.4"  Sandwich  Rate 
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2.3"  Rate 


Figure  6.9  Distribution  Plate  Rivet  Patterns  with  Dimensions 


The  rivets  are  V32"  countersunk  solid  rivets.  They  formed  an  excellent  connection 
between  the  sandwich  plates  and  levers. 

Some  problems  arose  from  the  welds  being  slightly  out  of  square.  Once  the 
individual  parts  were  fabricated,  they  were  assembled  with  tape  for  the  purpose  of 
mapping  out  the  drill  holes  for  connection  to  the  truss.  Once  the  holes  were  drilled,  only 
assembly  remained  before  testing  began. 
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One  of  the  built-in  “fudge  factors”  of  the  design  was  the  oblong  hole  pattern  in  the 
fastening  plates.  The  slack  that  existed  was  very  helpful  during  the  actual  bolting  of  the 
device  to  the  truss.  Everything  was  set  up  easily,  leaving  only  the  machine  screws  on  the 
bottom  to  tighten.  The  device  was  able  to  be  secured  tightly  to  the  truss  while  the  rods 
rotated  with  a  minimum  amount  of  friction. 

Applying  the  viscoelastic  material  to  the  aluminum  required  patience  and  care 
because  it  is  so  sticky.  In  the  application  of  the  damping  material,  air  bubbles  became 
trapped  at  the  interface  (this  problem  is  very  difficult  to  avoid,  especially  for  wide  areas 
such  as  the  large  sandwich  beam).  Most  of  the  bubbles  were  removed  by  pricking  them 
repeatedly  with  a  pin.  It  was  virtually  impossible  to  remove  all  of  the  air.  However,  a 
good  bond  was  achieved,  one  that  was  fr^  of  delamination  at  even  the  most  highly 
stressed  regions.  One  of  the  most  reliable  indications  of  debonding  is  the  presence  of  an 
audible  “crack”  or  “pop”  upon  deformation  of  die  beam.  The  beams  did  not  make  this 
sound  even  when  strained  far  beyond  the  limits  of  the  testing. 

Experimental  Procedure 

The  logical  sequence  for  the  experiment  was  to  obtain  all  data  on  the  improved 
treatment  before  proceeding  to  test  the  conventional  treatment.  As  stated  in  Chapter  5, 
three  viscoelastic  thicknesses  were  chosen  for  testing  because  the  theoretical  values  of  t2 
were  suspect  The  three  thicknesses  and  the  reasons  for  their  choice  are  listed  below. 

t2  =  O.OOS"  fell  in  the  theoretical  range  where  the  shear  strain  is  largely  dependent 
on  the  term  ( 1  +  hlt2). 

t2  =  0.015"  was  the  theoretical  minimum  for  the  dominant  term  in  the  expression  for 
shear  strain  (discussed  earlier). 

t2  =0.030"  is  in  the  theoretical  range  where  the  shear  strain  term  (/i  +  ^2)^  starts 
becoming  significant.  Recall  Figure  5.13  for  reference. 

The  basic  procedure  for  testing  all  configurations  (including  the  conventional 
treatment)  was  as  follows: 

1. )  Set  up  the  Tektronix  TekStar  software. 

2. )  Apply  the  accelerometer  to  the  structure  at  a  point  of  interest. 

3. )  Impact  the  structure  with  the  hammer. 


67 


4. )  Record  the  transfer  function  of  accelerance/force  over  a  given  frequency  band 
(for  ail  of  the  tests,  the  frequency  band  was  either  0-20  Hz  or  0-50  Hz).  This  ratio  is 
actually  the  ratio  of  output  power  spectral  density  over  input  power  spectral  density. 

5. )  Repeat  step  4  for  the  desired  number  of  trials. 

6. )  Input  data  to  the  SMS-StarStruct  software. 

7. )  Use  the  polynomial  curve  fit  feature  to  record  resonant  frequencies  and 
damping  ratios. 

8. )  Repeat  steps  2-7  for  different  accelerometer  and  force  locations. 

For  each  viscoelastic  thickness,  this  procedure  was  repeated  over  twenty  times,  each 
time  using  a  different  accelerometer/impact  location  combination.  The  locations  were 
chosen  based  on  the  positions  of  the  nodes  and  antinodes  of  the  various  modes  of  interest 

For  the  low  modes,  collecting  repeatable  transfer  function  data  was  difficult  when  a 
sharp  impact  was  used.  Because  the  power  spectral  density  of  a  sharp  impact  is  relatively 
uniform  over  the  entire  frequency  range  of  interest  the  low  modes  simply  did  not  receive 
enough  input  power  to  cause  large  deflections.  For  this  reason,  an  alternative  impact 
method  was  developed  for  achieving  good  low  mode  data.  If  the  hammer  tip  was  barely 
in  contact  with  the  structure  at  the  beginning  of  an  impact  and  accelerated  rapidly,  the 
sharpness  of  the  impulse  was  reduced  because  the  impact  occurred  over  a  longer  time 
period  (about  100-120  msec  for  a  “slow”  impact  as  compared  with  15-20  msec  for  a 
sharp  impact).  The  resulting  power  spectrum  was  a  more-or-less  flat  curve  out  to  about 
10  Hz,  at  which  point  it  rolled  off  by  several  orders  of  magnitude  out  to  20  Hz.  Despite 
this  rolloff,  the  slower  impact  gave  smooth,  consistent  transfer  function  curves  in  the  low 
frequency  range  (0-20  Hz)  that,  in  turn,  gave  accurate  curve  fits. 


Results 

The  results  from  the  testing  make  up  five  sets  of  data,  which  are  the  following; 

Set  1  -  Undamped  Modal  Frequency  and  Loss  Factor  Data 

Set  2  -  Improved  Treatment  Frequency  and  Damping  Data,  t2  =  0.005" 

Set  3  -  Improved  Treatment  Frequency  and  Damping  Data,  t2  =  0.015" 

Set  4  -  Improved  Treatment  Frequency  and  Damping  Data,  /2  =  0.030" 

Set  5  -  Conventional  Treatment  Frequency  and  Damping  Data,  t2  =  0.005" 
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The  reduced  data  is  presented  here  for  convenience.  The  raw  data,  which  is 
presented  in  the  form  it  was  actually  taken,  appears  in  Appendix  I.  The  only  difference  is 
that  the  modal  damping  ratios  were  converted  to  loss  factors  by  multiplying  by  two. 

Some  of  the  information  in  the  tables  is  abbreviated.  The  following  is  an 
explanation  of  what  everything  means: 

Test  -  a  number  of  impacts  whose  frequency  response  functions  are  averaged  into  a 
single  frequency  response  function. 

Response  and  Excitation  -  the  location  of  the  accelerometer  and  impact, 
respectively,  on  the  truss.  Zero  is  defined  at  the  root  of  the  truss,  while  unity  is  at  the 
truss  tip.  The  majority  of  the  response  locations  were  on  the  cross  members,  placed  at 
intervals  of  0.2L  down  the  truss.  The  “R”  or  “L”  indicate  placement  on  the  right  or  left 
side  of  the  truss.  If  no  “R”  or  “L”  is  present,  the  accelerometer  is  positioned  in  the 
center  of  the  cross  member.  The  presence  of  a  “Y”  indicates  that  the  accelerometer's 
active  axis  is  in  the  plane  of  the  truss.  If  no  “Y**  is  present,  the  active  axis  is  orthogonal 
to  the  truss  plane.  A  negative  sign  indicates  that  the  truss  was  impacted  on  the  back  side. 
This  was  only  done  when  something  was  in  the  way.  Finally,  the  SB  denotes  placement 
on  the  sandwich  beam  device  itself.  BSBIY,  for  example,  denotes  “bottom  of  large 
sandwich  beam,  active  axis  in  the  y-direction.”  MSBI  denotes  “middle  of  sandwich  beam 
1,  active  axis  in  the  z*direction.”  .6R  denotes  “located  60%  of  the  way  down  the  truss  on 
the  right  side,  active  axis  in  the  z-direction.” 

f/  -  frequency  of  mode  i 

%i  -  damping  ratio  of  mode  /,  expressed  in  percent  (this  number  is  the  result  of  the 
polynomial  curve  fit) 

Details  -  the  best  way  to  describe  this  is  to  use  an  example.  “IK,  5,  slow,  0-20” 
means  “1024  instantaneous  time  samples  per  impact,  five  slow  impacts,  frequency  band 
0-20  Hz.” 

The  blank  spaces  in  the  tables  mean  that  no  peak  was  present  or  that  the  frequency 
band  did  not  include  that  particular  mode. 

Note  that  the  loss  factor  scale  is  different  for  each  graph. 
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Figure  6.10  Loss  Factors,  Undamped  Truss 
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Observations  of  Results 


Collecting  the  loss  factor  data  for  the  improved  treatment  was  difficult  because  the 
addition  of  the  damping  device  brought  quite  a  number  of  new  system  modes  into  the 
frequency  band  from  0-50  Hz.  The  data  is  very  consistent  with  regard  to  frequency,  but 
there  is  a  wide  scatter  of  loss  factor  values  for  most  of  the  modes  of  interest.  The  most 
prominent  new  modes  introduced  to  the  system  were  the  following: 

-  fundamental  sandwich  beam  bending  modes  (above  30  Hz) 

-  fundamental  “rocking”  modes  of  the  entire  sandwich  beam  device  (this 
motion  was  in  a  plane  parallel  to  the  truss  plane,  with  sandwich  beam  1 
resonating  at  about  14  Hz  and  sandwich  beam  2  resonating  at  about  27  Hz) 

-  in-plane  rotational  mode  of  the  entire  sandwich  beam  device,  where  the  axis  of 
rotation  (orthogonal  to  the  truss  plane)  occurred  at  the  middle  set  of  levers 

Some  of  the  impact  locations  excited  these  modes  more  than  the  modes  of  interest,  and 
some  of  the  loss  factor  data  points  are  grouped  at  two  or  more  distinct  values  r).  It  is 
believed  that,  for  some  data  points,  the  curve  fitter  found  these  modes  instead  of  the 
modes  of  interest.  For  example,  the  fourth  bending  mode  of  the  truss  (at  about  27  Hz) 
was  believed  to  be  one  of  the  more  highly  damped  modes,  yet  the  scatter  of  points  ranges 
from  low  to  high.  For  all  of  the  improved  treatment  graphs  (Figures  6.11-13),  the 
frequency  of  27  Hz  has  two  distinct  groupings.  A  look  at  the  data  sets  indicates  that  the 
low  damping  values  were  recorded  during  impacts  that  excited  torsion.  Intuitively,  the 
rocking  back  and  forth  of  the  damping  device  should  produce  negligible  damping. 
Therefore,  it  was  concluded  that  the  lowly  damped  groupings  are  attributable  to  the 
resonance  of  sandwich  beam  2.  For  a  similar  example,  consider  the  scatter  of  the  third 
truss  bending  mode  and  the  relatively  tight  grouping  of  a  few  points  well  below  the  rest 
of  the  mode  3  data.  They  too  are  the  result  of  torsion-inducing  impacts,  but  are 
attributable  to  the  rocking  of  sandwich  beam  1. 

Other  than  these  two  examples,  the  data  is  relatively  good.  An  exception  to  this 
conclusion  is  that  the  mode  1  data  is  highly  scattered  with  regard  to  loss  factor.  This 
scatter  is  caused  by  a  number  of  factors.  First,  the  frequency  band  was  from  0-20  Hz  for 
most  of  the  mode  1  data  points.  Noting  Table  6.4,  two  of  the  three  lowest  values  of  rji 
were  procured  from  tests  with  bandwidth  0-50  Hz.  1  Hz  is  on  the  low  edge  of  this  band, 
and  the  frequency  response  function  was  exceptionally  “ragged”  in  this  region  (a 
condition  which  made  for  an  unreliable  curve  fit).  Second,  the  accelerometer 
specification  for  the  lower  frequency  limit  is  0.5  Hz,  but  the  output  begins  to  roll  off 
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somewhere  between  0^  Hz  and  1  Hz.  The  effect  of  this  rolloff  is  probably  a  decrease  in 
the  signal-to-noise  ratio  (which  causes  a  random  distribution  of  damping  values),  but  it 
was  considered  to  be  small  because  the  slow  impact  technique  provided  ample  power  in 
the  low  frequency  region  to  obtain  a  strong  output  voltage. 

It  is  important  to  mention  here  that  the  data  points  should  not  be  averaged  into  a 
single  value.  For  any  given  mode,  each  data  point  represents  the  damping  value  recorded 
for  a  particular  combination  of  accelerometer  location  and  impact  location.  If  the 
accelerometer  or  impact  location  is  near  a  node,  the  output  voltage  will  be  small,  and  the 
corresponding  curve  fit  will  yield  a  less  reliable  damping  value  than  if  the  output  signal 
were  higher.  In  order  to  be  impartial,  as  many  data  points  were  recorded  as  could  be 
reasonably  obtained  for  any  given  frequency  response  function.  Curve  fitting  all  possible 
data  served  the  dual  function  of  helping  discern  sensitivities  in  the  testing  procedure. 
Hnally,  even  though  the  truss  was  relatively  massive  compared  to  the  accelerometer,  the 
sway  of  large  truss  deflections  caused  the  accelerometer  cable  to  influence  the  output. 
Appendix  C  addresses  this  particular  problem.  It  should  be  emphasized  that  this 
particular  problem,  although  annoying,  did  not  significantly  affect  the  data. 

The  data  from  the  conventional  treatment  (Figure  6.14)  was  very  good  for  all  modes 
except  the  “fourth  bending”  mode.  It  was  reputed  that  the  two  distinct  groupings  are 
actually  two  modes.  The  undamped  truss  data  showed  the  fourth  bending  mode  at  32  Hz, 
while  the  third  torsion  mode  was  determined  to  be  about  40  Hz.  The  effect  of  adding  the 
conventional  treatment  caused  the  fourth  bending  mode  to  shift  up  to  about  34  Hz  due  to 
a  stiffening  of  the  long  members  in  out-of-plane  bending.  Conversely,  because  the 
conventional  treatment  was  applied  at  the  furthest  distance  possible  from  the  axis  of 
torsion,  the  added  mass  caused  the  third  torsion  mode  to  shift  down  to  about  35  Hz. 

Conclusions  Based  on  the  Experimental  Results 

It  is  clear  from  viewing  the  loss  factor  graphs  that  the  damping  device  did  not 
achieve  the  levels  of  modal  damping  that  the  theory  predicted.  The  predicted  damping 
values  were  believed  to  be  too  high  before  any  of  the  data  was  collected,  but  they  were 
listed  at  the  end  of  Chapter  5  as  calculated.  It  is  obvious  from  reviewing  Table  5.1  that 
the  energies  in  the  viscoelastic  layer  predicted  by  the  Matlab  program  completely 
dominate  the  finite  element  prediction  of  truss  strain  energy.  Based  on  the  expected 
relative  deflections  of  the  pin  joint  rod  axes,  it  was  difficult  to  believe  that  such  a  small 
volume  of  material  could  possess  such  a  large  fraction  of  the  total  strain  energy. 
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Rod  Axis  Displacement  (In ) 


There  were  a  number  of  suspected  reasons  why  the  Matlab  computer  program  gave 
such  high  vales  of  strain  and  damping  energy.  Test  results  for  tn  =  0.015"  gave  damping 
ratios  that  were  much  lower  than  expected.  It  was  initially  presumed  that  the  high 
theoretical  values  for  strain  energy  were  the  result  of  an  incorrectly  determined  moment 
Ml).  The  moment,  taken  from  the  finite  element  static  solution  described  in  Chapter  5, 
was  supposedly  too  high  because  of  the  load  P  on  the  short  lever  tips.  The  only  way  to 
verify  this  suspicion  was  to  conduct  a  small  test  to  determine  the  actual 
force/displacement  relationship  between  the  lever  rods. 

A  test  apparatus  was  set  up  to  determine  the  behavior  of  the  sandwich  beams  in 
both  tension  and  compression.  While  the  beams  were  gradually  loaded  with  calibrated 
weights,  the  deflection  of  the  rod  axis  was  recorded  with  the  aid  of  a  traveling 
microscope.  The  microscope  is  a  precision  instrument  which  uses  an  eddy  current  sensor 
to  determine  the  displacement.  The  machine  gives  results  to  within  0.00005". 

With  the  t2  -  0.005"  sandwich  beams  fastened  to  the  truss,  a  number  of  readings 
were  taken  for  various  tip  deflections.  The  truss  was  loaded  by  a  point  load  at  the  center 
of  the  second  cross  member  from  the  bottom  (.8L)  to  obtain  the  desired  tip  deflections. 
Figure  6.15  shows  the  rod  axis  displacements.  This  data  gave  known  displacements  of 
the  lever  tips  for  a  known  truss  tip  displacement.  The  1  and  2  denote  different  trials. 
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Figure  6.15  Tip  Displacement  vs.  Lever  Rod  Axis  Vertical  Displacement 
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After  collecting  this  data,  the  sandwich  beams  were  removed  from  the  truss.  Then 
they  were  loaded  in  tension  and  compression  with  a  known  load  in  order  to  obtain  the 
force/displacement  relationship  between  the  lever  tips.  Figures  6.16-6.19  (next  page) 
show  graphically  the  results  of  these  tests.  From  these  graphs,  it  is  apparent  that,  indeed, 
the  sandwich  beams  behave  quite  differently  in  tension  than  compression.  The  theory 
also  predicts  a  divergence  of  the  cycle  halves  as  discussed  in  Chapter  5  (see  Figure  5.8 
for  the  theoretical  curves). 

The  real  information  to  deduce  from  analyzing  the  force/displacement  relationship 
is  noting  at  what  load  the  divergence  between  cycle  halves  becomes  significant.  It  seems 
that  the  appropriate  level  of  caution  was  taken  in  determining  the  geometric  parameters  h, 
s,  t,  and  the  cross-section. 

The  next  step  was  to  use  a  best  guess  linear  approximation  of  all  the  data  for  a  given 
configuration.  Both  tension  and  compression  data  was  taken  into  account  for  this  guess. 
The  resulting  slope  is  the  effective  spring  constant  between  the  rods. 

Once  the  spring  constant  was  determined,  finding  the  actual  load  P  between  the 
rods  for  a  given  truss  tip  displacement  was  a  simple  matter  of  multiplying  kgffhy  the 
actual  rod  displacement.  The  results  of  this  “mini-study”  are  tabulated  here  for  a  four- 
inch  truss  tip  displacement.  The  value  of  the  rod  axis  displacement  A  is  based  on  the 
actual  relative  rod  displacements  for  the  /2  =  0.005"  configuration.  It  was  believed  the 
actual  displacements  of  the  /2  =  0.030"  configuration  would  be  very  close.  This  belief  is 
based  on  the  fact  that  none  of  the  configurations  stiffened  up  the  truss  significantly  more 
than  any  other. 


f2  =  005" 

A(in) 

P 

t2  =  .030" 

A(in) 

P 

Sbi 

k  =  14  Ib/in 

.145" 

2  lb. 

k  s  12  Ib/in 

.145" 

1.71b. 

SB  2 

k  =  5.3  Ib/in 

.17" 

.9  lb. 

k  *  4.7  Ib/in 

AT 

.8  lb. 

Table  6.6  Sandwich  Beam  Force/Displacement  Experimental  Data 

After  the  actual  force/displacement  relationships  were  found,  a  modified  finite 
element  analysis  was  conducted  using  this  data.  The  sandwich  beam  and  short  lever 
components  in  each  bay  were  replaced  by  a  single  member  whose  “effective”  stiffnesses 
was  backed  out  by  using  the  relation 

AE  g 

F*  where  (6.1) 
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Force  between  axis  rods  (lb.) 


Figure  6.17  Force-Displacement  Relationship  -  Sandwich  Beam  2,  tz-  0.005" 
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Displacement  (in )  Displacement  (In.) 


Figure  6.18  Force-Displacement  Relationship  -  Sandwich  Beam  1 ,  =  0.030" 
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Figure  6.19  Force-Displacement  Relationship  -  Sandwich  Beam  2,  t2  =  0.030" 
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Table  6.7  shows  the  tabulated  results  of  two  static  models,  one  for  each  /2.  The  dynamic 
analyses  were  not  useful  for  displacements  because  the  maximum  displacement  (in  this 
case  the  truss  tip  elements)  is  default  set  to  one  meter.  However,  the  distribution  of  strain 
energy  in  the  dynamic  analyses  is  almost  identical  to  the  static  load  results  (±0.5%).  In 
addition,  the  percentage  of  total  strain  energy  contained  in  the  replacement  member 
conilguration  is  close  to  the  sum  of  the  strain  energy  in  the  sandwich  beam  and  short 
lever  elements  of  the  “accurate”  models. 


configuration 

A(in) 

P 

SB1.  f2  =  0.005" 

.124- 

1.71  lb. 

SB  2.  f2  =  0.005" 

.120- 

.579  lb. 

SB1.  f2  =  0.030" 

.131- 

1.51  lb. 

SB  2.  f2  =  0.030" 

.122- 

.520  lb. 

Table  6.7  Modified  Finite  Element  Model  Results 

Based  on  this  data,  it  was  concluded  that  the  moment  obtained  by  the  finite  element 
analysis  is  not  the  source  of  the  high  theoretical  shear  strains. 

Another  possible  reason  for  the  discrepancy  between  theory  and  experiment  is  the 
assumption  that  the  flexural  rigidity  D  is  independent  of  t2.  Clearly,  when  viscoelastic 
material  bonds  two  sandwich  plates  together,  the  effective  stiffness  of  the  composite 
structure  is  stiffer  in  bending  than  two  sandwich  plates  with  nothing  between  them  to 
resist  the  shear  forces. 

To  investigate  this  phenomenon,  another  small  test  was  conducted.  Two  sandwich 
beams,  whose  outer  plates  each  measured  18"  x  1.00"  x  0.032",  were  fabricated  in  order 
to  relate  the  change  in  stiffness  with  changes  in  the  viscoelastic  thickness.  The  two 
viscoelastic  thicknesses  tested  were  0.005"  and  0.015".  Each  beam  was  cantilevered 
(both  plates  clamped)  and  impacted  several  times  to  determine  the  fundamental 
frequency.  Since  the  fundamental  frequencies  were  close  to  one  another,  it  was  assumed 
that  the  shear  modulus  of  the  viscoelastic  was  constant.  With  an  experimentally 
determined  value  of  fo,  the  effective  moment  of  inertia  of  the  composite  beam  can  be 
found  by  using  the  frequency  equation 


where  =  1.875  for  a  cantilever  beam  vibrating  in  its  fundamental  mode  and  E  is  the 
elastic  .modulus  for  aluminum.  We  can  back  out  leff  as 
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(63) 


This  expression  was  used  to  calculate  the  effective  stiffnesses  displayed  below. 


fieq. 

(Hz) 

undamped 

n 

undamped 

lz(in») 

freq. 

(Hz) 

i|for 

O.OOS" 

lztor^  = 
0.005"  ( in*) 

freq. 

(Hz) 

Tjfor 
=  0.015" 

lzlbr^  = 
0.015"  (in*) 

6.8 

i.68E4J5 

4.147 

0.01282 

5.97E-06 

6.91 

0.215 

1.70862E-05 

6.605 

0.3096 

1.66E-05 

6.82 

0.211 

1.6644E-05 

6.604 

0.3072 

1.66E-05 

6.97 

0.192 

1.73842E-05 

6.685 

0.315 

1.7E-05 

6.82 

0.215 

1.6644E-05 

6.76 

0.3136 

1.73E-05 

6.89 

0.229 

1.69874E-05 

6.82 

0.222 

1.6644E-05 

6.89 

0.221 

1.69674E-05 

Table  6.8  Frequencies,  Loss  Factors  and  Effective 
Moments  of  Inertia  of  Test  Beams 


For  both  viscoelastic  thicknesses,  the  effective  area  moment  of  inertia  is  roughly 
tripled.  This  observation  raised  the  question  of  what  would  happen  if  these  values  of  / 
were  used  in  the  Matlab  program  in  place  of  the  flexural  rigidity  D.  This  is  a  straight 
substitution  because  the  test  beams  have  unit  width  and  the  same  thickness  as  the 
sandwich  plates.  Interestingly,  this  simple  change  reduced  the  strain  energy  to  what  can 
be  considered  a  reasonable  magnitude.  Table  6.9  lists  these  results. 


h 

D,  (Ib.in.) 

Pj  (Ib.in.) 

Ut+Msb 

Uvei 

tivte 

n 

- — 

cr25 

S!0B32 

0.^ 

I513S7 

.015* 

0.316 

0.117 

0.65 

0.0837 

0.0309 

0.0901 

.030" 

0.831 

0.307 

0.65 

0.22 

0.0814 

0.1904 

Table  6.9  Loss  Factor  Calculations  Using  Tripled 
Flexural  Rigidity  Dand  Old  FE  Model  Data 


Notice  that  the  strain  energy  value  in  the  previous  results  is  from  the  “accurate”  FE 
model.  In  addition,  the  old  FE  loads  (2.42  lb.  and  0.96  lb.  for  sandwich  beams  I  and  2, 
respectively)  were  used  in  the  Matlab  program.  When  the  modified  FE  strain  energy 
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calculation  (0.78  lb.  in.  for  the  replacement  member  model)  and  the  experimentally 
determined  loads  P  are  used,  the  results  become  the  following: 


h 

P»(ib.) 

0,  (b.in.) 

Psiti) 

Cb  (l>  in ) 

IMUse 

Uvei 

n 

17 

.m 

.8 

7Bg 

.015" 

1.82 

.1785 

.84 

.0893 

.776 

.0473 

.0237 

.050 

.030" 

2.0 

.568 

.9 

.270 

.778 

.1505 

.0715 

.132 

Table  6.10  Loss  Factor  Calculations  Using  Tripled  Flexural  Rigidity  D, 
Modified  FE  Model  Data,  and  Experimentally  Determined  Loads  P 

It  should  be  mentioned  that  the  loads  P  from  the  FE  model  (Table  6.7)  are  lower 
than  the  experimentally  determined  values  used  here.  The  effect  of  this  discrepancy  is 
that  the  tabulated  value  of  elastic  strain  energy  Uj-^  would  be  lower  for  the  given 
loads  P,  giving  a  still  lower  loss  factor.  In  any  case.  Tables  6.9  and  6.10  illustrate  the 
sensitivity  of  the  results  to  the  assumed  flexural  rigidity  D. 

The  development  of  the  sixth-order  equation  of  motion  (Appendix  D)  shows  how 
the  extension  of  the  constraining  layers  are  converted  from  expressions  in  terms  of  u  to 
expressions  in  terms  of  w.  It  is  difficult  to  find  the  real  effect  of  D  on  these  terms 
because  the  derivation's  intermediate  steps  convert  the  w-derivative  terms  into  expressions 
that  have  no  physical  meaning.  When  the  w-derivative  term  (Equation  D.23)  is 
substituted  into  Equation  (D.22),  a  D  shows  up  in  the  denominator  of  the  coefficient. 
Therefore,  even  though  the  flexural  rigidity  D  is  a  number  based  on  bending  phenomena, 
this  quantity  affects  the  value  of  extensional  terms. 

The  experimental  data  from  the  conventional  treatment  agrees  very  well  with  the 
theory.  The  expected  results  at  the  end  of  Chapter  4  list  only  the  first  two  modes.  In  fact, 
the  loss  factors  for  the  first  four  bending  modes  are  close  to  the  theoretical  curves.  Table 
6. 1 1  compares  the  theory  with  the  results  of  the  experiment. 


Mode  Number 

Theoretical 
Prediction  of  n 

Experimental 
Result  of  rj 

1 

?jaO.09 

2 

Tja014 

ua013 

3 

?ja0145 

?7s014 

4 

n«0135 

TjaOIS 

Table  6.1 1  Comparison  of  Theoretical  vs.  Experimental 
Damping  Values,  Conventional  Treatment 
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It  should  be  noted  that  the  theoretical  values  are  only  a  rough  estimate  based  on 
Rgure  6.14.  Considering  the  number  of  simplifying  assumptions  that  were  made 
regarding  this  analysis,  these  results  are  about  as  close  as  one  can  expect 
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VII .  Recommendations  for  Further  Research 


The  organization  and  nature  of  this  thesis  project  leaves  many  ‘^loose  ends”  that  can 
be  individually  investigated.  Several  of  the  conclusions  based  on  the  experimental  results 
from  the  previous  chapter  are  really  questions  because  the  author  ran  out  of  time.  For  a 
thesis,  the  presentation  of  the  results  is  rather  unorthodox.  However,  for  a  scientific 
research  endeavor,  the  issues  that  were  addressed  were  a  natural  consequence  of  the 
disparity  between  the  theory  and  the  experimental  results.  Truthfully,  some  of  the 
problem  issues  are  unanswered. 

Some  of  the  assumptions  made  in  the  analysis  are  suspect  and  warrant  further 
investigation.  The  assumed  shape  of  the  sandwich  beam  in  its  fundamental  bending 
mode  was  a  sine  wave.  This  problem  has  already  been  addressed  in  Chapter  5.  Another 
possible  source  for  error  is  in  the  assumed  boundary  conditions.  As  Cottle  [22:  18] 
observed,  a  three-layer  sandwich  beam  actually  has  eight  boundary  conditions,  two  of 
which  go  to  zero  because  the  difference  of  the  extensional  displacements  ui  is  considered 
instead  of  their  absolute  values.  However,  in  the  sandwich  beam  given  in  this  study, 
there  is  some  question  as  to  which  two  boundary  conditions  go  away.  As  discussed  in 
Chapter  5,  the  assumed  extensional  boundary  conditions  were  zero  strain  at  the  ends. 
Actually,  only  one  of  the  sandwich  plates  has  zero  longitudinal  strain  at  either  end.  What 
about  the  boundary  condition  on  the  other  plate  (the  plate  that  must  transfer  the  moment 
to  the  three-layer  interior)?  Recall  that  the  problem  was  divided  into  two  parts;  the  direct 
shear  static  case  and  the  total  (static  plus  dynamic)  bending  case.  It  seems  that  if  the 
tension  in  the  sandwich  plate  is  considered  in  one  part,  it  ought  to  be  disregarded  in  the 
other  part  so  as  not  to  contribute  to  the  solution  twice. 

Perhaps  the  most  neglected  area  of  this  project  was  the  procedure  for  finding  the 
dimensions  of  the  damping  device.  It  was  believed  that  this  part  of  the  project  was 
relatively  simple,  and  not  much  time  was  spent  on  it  in  consequence.  This  was  a  serious 
error  in  judgment  because  optimizing  the  **non-damping”  parts  would  have  allowed 
considerably  more  damping  material  to  be  applied.  However,  choosing  smaller 
dimensions  for  these  parts  might  have  invalidated  the  theory  because  of  the  asymmetric 
bending  phenomena.  Clearly,  this  is  a  more  difficult  problem  than  a  cursory  look  would 
indicate.  Using  the  finite  elements  turned  out  to  be  a  convenient  way  of  finding  strain 
energy  distributions  and  iterating  through  several  candidate  configurations.  However,  the 
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method  for  eliminating  parameters  was  not  at  all  an  optimal  technique.  Much  work  needs 
to  be  done  on  this  problem  alone. 

One  positive  observation  from  the  “improved”  treatment  is  that  a  measurable 
amount  of  damping  was  achieved  despite  the  relatively  small  area  of  coverage  (with 
respect  to  the  conventional  treatment).  The  area  of  viscoelastic  coverage  of  the  sandwich 
beams  is  about  96  in.-  as  compared  with  280  in.-  for  the  conventional  treatment.  This 
difference  is  significant  and  probably  has  a  lot  to  do  with  the  disparity  in  modal  loss 
factors  between  treatments.  For  some  existing  structures,  it  may  not  be  possible  to  cover 
as  much  area  as  the  conventional  treatment  does. 

Another  area  of  this  study  that  was  not  treated  with  appropriate  diligence  is  the 
transfer  of  strain  energy  with  structural  additions  like  the  damping  device.  Perhaps  the 
reason  for  this  is  that  the  idea  for  what  the  damping  device  would  look  like  was  decided 
early  in  the  project,  long  before  this  problem  was  discovered.  This  type  of  premature 
decision-making  is  not  a  very  sound  design  philosophy.  In  retrospect,  it  appears  that 
forcing  the  strain  energy  into  a  desired  region  by  adding  stiffness  is  just  as  important  as 
designing  the  actual  damping  device.  The  idea  of  using  added  stiffness  in  certain  regions 
as  part  of  the  damping  treatment  is  certainly  one  that  could  produce  fruitful  results.  For 
this  particular  project,  consider  what  would  happen  if  the  second  sandwich  beam  and 
levers  were  eliminated  and  all  of  their  weight  was  used  to  stiffen  up  the  four  truss  bays 
furthest  from  the  built-in  end.  This  concept  can  be  best  described  as  designing  in  a 
structural  “weak  link**  where  the  vast  majority  of  the  total  strain  energy  resides  (of 
course,  the  weak  section  would  still  have  to  be  fail-safe.  This  would  be  an  excellent 
study. 

A  design  technique  that  was  not  considered  was  the  partial  coverage  concept 
presented  in  Nokes  and  Nelson's  1968  paper  [23].  To  be  absolutely  fair,  this  idea  should 
have  been  employed  in  a  competing  configuration.  Certainly,  this  option  could  be  used 
in  concert  with  other  ideas  like  Plunkett  and  Lee’s  segmenting  technique. 

For  this  thesis,  ail  the  joints  were  modeled  and  constructed  as  ideal  connections. 
Any  design  that  would  actually  be  used  would  most  likely  involve  damped  joints  as  well 
as  damped  beams.  One  possible  study  might  be  to  compare  the  effect  of  adding  damping 
in  the  joints  as  well  as  in  the  main  part  of  a  damping  addition. 

Finally,  it  is  unfortunate  that  the  dynamic  terms  in  this  project's  analysis  were 
insignificant  compared  to  the  quasi-static  terms.  It  is  possible  to  conceive  a  structure 
where  the  dynamic  response  of  the  sandwich  beam  is  a  major  factor.  Verifying  the 
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dynamic  part  of  a  time-dependent  boundary  condition  solution  with  an  experiment  is  a 
problem  worth  studying. 
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V///.  Conclusion 


Experimental  results  show  that  the  sandwich  beam  damping  device  did  not  exceed 
the  conventional  damping  treatment  for  any  modes  of  truss  vibration.  Modal  loss  factors 
for  the  conventional  treatment  were  approximately  three  times  higher  than  the  sandwich 
beam  device.  The  conventional  treatment  was  well-designed  and  well-constructed,  and 
the  experimental  results  verify  that  Plunkett  and  Lee's  segmented  constraining  layer 
analysis  is  accurate. 

This  project  turned  out  to  have  complexities  (both  theoretical  and  experimental)  that 
were  not  anticipated  at  the  outset.  Even  though  the  sandwich  beam  damping  device 
appears  simple,  it  is  not.  A  number  of  separate  analyses  are  required  to  describe  the 
oscillatory  motion  of  the  device. 

The  expression  for  the  displacement  of  the  three-layer  sandwich  beam  is  quite 
complicated,  and  is  only  an  approximate  solution  since  part  of  the  solution  is  a 
comparison  function  that  satisfies  the  boundary  conditions  but  not  the  differential 
equation.  The  accuracy  of  the  solution  is  further  reduced  by  the  assumption  of  simple 
mode  shapes  that  made  finding  a  solution  in  series  form  possible.  Despite  its  complexity, 
expressions  were  obtained  for  the  shear  strain  of  the  viscoelastic  layer,  which  is  a 
function  of  the  derivatives  of  the  displacement  w.  For  the  given  sandwich  beam 
configuration,  a  single  term  dominated  the  expression  for  shear  strain.  This  term  caused 
problems  in  interpreting  the  results  of  a  computer  program  that  performed  numerical 
integration.  The  program  gave  relatively  high  values  of  viscoelastic  strain  energy  and 
damping  energy.  Although  this  problem  was  not  solved,  two  possible  reasons  for  the 
high  values  were  investigated  by  conducting  small-scale  experiments.  The  result  of  these 
investigations  were  that  the  problem  with  the  expression  for  shear  strain  was  not  found. 
One  possible  reason  why  this  problem  could  not  be  solved  is  that  Mindlin  and  Goodman's 
method  for  finding  displacement  solutions  to  problems  with  time-dependent  boundary 
conditions  is  mathematically  ill-suited  for  using  information  from  its  higher  derivatives. 
It  is  apparent  that  the  accuracy  of  the  comparison  functions  (which  are  really  truncated 
series  solutions  of  sinusoids)  degrades  considerably  as  more  and  more  derivatives  are 
taken. 
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Appendix  A  -  Characterization  of  the  Test  Article  and  Finite  Element  Model 


The  purpose  of  characterizing  the  undamped  response  of  the  test  article  was  threefold; 
the  system  natural  frequencies,  mode  types,  and  undamped  modal  loss  factors  had  to  be 
found.  These  three  sets  of  data  comprise  die  adequate  characterization  of  the  test  article. 

The  finite  element  model  described  in  Chapter  2  was  used  as  a  comparison.  Some 
modes  are  close  in  frequency,  and  the  model  helped  spot  frequency  regions  of  high  modal 
density. 

Several  trials  were  needed  for  the  purpose  of  distinguishing  mode  types.  For 
example,  if  the  accelerometer  is  placed  on  one  of  the  tip  comers,  the  transfer  function  will 
register  both  out-of-plane  bending  and  torsional  modes.  The  two  were  distinguished  from 
one  another  simply  by  performing  other  trials  where  the  accelerometer  was  placed  on  the 
longitudinal  axis  (in  the  center  of  a  cross  member).  Additional  trials  were  conducted  for 
determining  the  in-plane  bending  modes  by  flipjnng  the  active  direction  of  the  accelerometer 
and  impacting  the  truss  both  in  and  out  of  the  plane.  The  results  are  tabulated  below  for  the 
first  nine  modes.  The  modal  loss  factors  were  determined  only  for  the  first  six  modes.  The 
curve  fitter  in  the  SMS  StarStract™  software  package  was  employed  for  this  purpose.  Its 
output  was  in  percent  damping,  and  the  tabulated  results  shown  as  17  are  double  the  value  of 
percent  damping. 


Percent 

Mode  modal  FE  predicted  Difference 


Mode 

Freq.  (Hz) 

Type 

n 

frequency 

(1-FE/exp,)*100 

1 

1.05 

Bending  1 

0.008 

0.955 

9 

2 

6 

Bending  2 

0.003 

5.94 

1 

3 

7.4 

Torsion  1 

0.002 

7.98 

-7.3 

4 

16.5 

Bending  3 

0.0015 

16.5 

0 

5 

22 

In-plane  1 

0.0015 

24.2 

-10 

6 

22.5 

Torsion  2 

0.0005 

23.2 

-3.1 

7 

31.9 

Bending  4 

0.001 

31.7 

0.6 

8 

37.6 

Torsion  3 

41.3 

-9.8 

9 

50 

Bending  5 

49.3 

1.4 

Table  A.  1  Comparison  of  Finite  Element  Model  and  Experimental 

Modal  Frequencies 


The  modal  loss  factors  are  averaged  over  three  trials  where  the  force  was  applied  gradually 
over  about  100  ms.  The  maximum  truss  tip  deflection  was  about  four  inches. 
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Appendix  B  -  Analysis  of  a  Five  Layer  Sandwich  Beam  Using  the  Variational  Calculus 

The  following  is  taken  from  Torvik's  paper  [17].  Some  of  the  steps  omitted  in  his 
paper  are  included  for  clarity,  and  a  few  simplifications  are  made  in  the  application  of  this 
development  to  the  2D  truss  structure. 

Consider  the  following  slice  dx  of  a  layered  beam. 


dx 


Rgure  B.  1  Differential  Slice  c^of  a  Layered  Beam 

For  this  development,  assume  that  the  thicknesses  remain  constant  so  that  Mxf)  is  the  same 
for  all  layers.  Assume  the  shear  layers  do  not  carry  any  bending  or  axial  load  (because  G/ « 
£j).  Finally,  the  beam  is  assumed  to  be  of  unit  width. 


Hamilton's  principle  will  be  used  to  formulate  the  equations  of  motion  and  the 
boundary  conditions. 

The  elastic  strain  energy  due  to  bending  is 


(B.l) 
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where 


‘  12(1- v^) 

£«  Young's  Modulus 
V-  poisson’s  ratio 

w-  the  transverse  displacement  of  the  cross  -secticHi 
-  the  thickness  of  layer  i 


In  unabridged  form, 


Jo2[axl  Jo2\dxl  ^^)o2\dxl 


The  strain  energy  in  the  viscoelastic  layers  due  to  shear  is 

J  0^\  \  2^2 

The  kinetic  energy  of  a  slice  dx  of  the  composite  structure  is 

T'  -  (p,<l  +  Pnji  +  P/3  +PA  +  P/5)  ^ 


(B.2) 


(B3) 


(B.4) 


Since  we  are  using  this  analysis  to  study  the  2D  truss,  we  must  deal  with  the  cross 
members.  We  will  assume  they  act  as  rigid  bodies  (having  zero  strain  energy  for  all  time). 
For  the  frequency  range  of  interest,  this  is  a  good  assumption.  Therefore,  Uxm  =  0. 
However,  they  do  possess  kinetic  energy  equal  to 
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T 

x^m 


(5)’ 

jr-14.03 

n 

f— V 

[  dtj 

+  jPiA^, 

z.56.25 

f— VI 
\  */ 

*x 

b>42J6 


(B.5) 


I* -70.25 


where  -  the  length  of  each  cross  -  member  (13.0  in.) . 

These  discrete  masses  on  a  continuous  structure  will  prevent  us  from  being  able  to 
assume  an  orthogonal  modal  expansion,  but  we  can  simplify  to  orthogonalize  them  at  the 
expense  of  some  accuracy.  To  do  this,  a  modified  mass  density  p'  was  used.  Since  the 
cross-members  possess  roughly  half  of  the  truss's  mass,  the  modified  mass  density  p'  is 


p'=\^ 


(B.6) 


Since  the  loss  factor  T]s  will  be  shown  to  be  independent  of  p',  the  effect  on  the  accuracy  of 
the  solution  will  be  minimal.  A  modified  finite  element  model  was  constructed  to  compare 
the  result  of  using  p'  with  massless  cross-members  versus  the  “accurate”  FE  model. 
Indeed,  the  natural  frequencies  are  minimally  affected  for  the  first  five  nnodes  (±4%  between 
models).  Therefore,  this  was  considered  to  be  a  safe  assumption. 

Hamilton's  principle  is  the  following: 

a/^(L+ HO  t// -0  (B.7) 

where 

W  -  external  work  terms 

If  we  let 


D 

(PA) 


-  D,  +  +  Dj 

-  (P.  h  +  P-h  +P3^3  +  PaU  +  P5^) 


^4 


2/, 


2/4 
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substitute  p*  for  p,  and  consider  only  one  of  the  long  members  in  pure  out-of-plane 
bending  (no  torsion),  then 


6 


JoU  2  V^j  2  ”[dxj 

+  f  T )(■—]  ^  +  f +  Viv  +  F.m,  +  F3i^  +  F5U5] 

J  o2  \  St )  \  Sx  f 


(B.8) 


dt^Q 


Taking  the  variations  with  respect  to  the  ten  dependent  variables 

d^w  dw  dw  du.  du,  ,  du. 

Tpr.  — ’  T’  “t»  “j’  “5»  ”rr»  if* 
<?jc  dx  dt  dx  dx  Sx 


•‘s 

. 

-U 

« 

o\ 

+  G.,t^ 

f 


u,-u,  ,  dw 

- '■  +  d^  — 


\  ^2 


— 6mj  -  —  dw,  +  d,6 

/,  ~  Sx 


Sxfyt^  *2 


+G4/4 


,  dw 
-i - -  +  d, — 

\  ^4  ; 


dtc 


•  (Pi^c) — d~  fic  +  l  M6  —  -t-Vdw ^  Ffiu^  +Fy6u^  +  F^dw, 
Jo*’* 


dt 


dx 


(B.9) 


dt  =  0 


In  order  to  group  these  expressions  into  their  respective  spatial  variables,  we  must  integrate 
the  differential  variations 


«i!;,  «ia., 

dx-  dx  dt  dx  dx  dx 


by  parts.  Taking  each  term  separately. 
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(B.IO) 


-  D - r< 

Jo 


D — -6 — rdc-D — rd —  +  D — -6 — dx 

dx^  dx^  dxr  dbcL  Jo  o!*^  dx 


- 1  -  -E,/,  ^6u^  +  f  £;/,  -^  du,  <& 

\  "  dx  dx  <?JC  ^lo  Jo  dx-  ' 


■i: 

-}>•> 

-i: 


^  T„  Jo  ai- 


^6^ 

dx  dx 


^  du,  -  du.  , 

E,/,-^d— 

ac 


_  f  d^u.  .  , 

-^5^5 dx 

qX  b  J  0  M 


-1, 

-i: 


0  '‘I 


\/ 


/\‘2 


1  o  1  o  .  C 

-dMj  --d«,  +ti2d  — 


dx 


-G, 


1,  1.  1-  Ic  ,  dw  .  ,  dw  _ 

■-u,du,---u,6uj-—u^du^  +-M,6a,  +^4  — dw, — fin, 


\^2 


5x 


<?x 


i  J  j2.  dw\ 

+  d^u,o - Ojitd  —  +  a,/,  —  d —  iw: 

dx  dx  -  -  dx  dx) 


The  last  integral  has  itself  three  integrations  by  parts,  which  are 

-  f —dx  — G,4M,dwl^  +  f  G.d^^6wdx 

Jo  ~  -  dx  ~  ~  J  0  ~  -  dx 

[  G^dyUyd^dx  »Gn<LMiSw^  -  [  G^dy-^Sw  dx 
Jo*  dx  J  o  dx 

-  f  G,c/,'/,  ^6—dx  -  -Gydy\  ^dJ  +  f  GM.-ty  ^dwdx 
J  0  '  ~  ~  dx  dx  '  ~  '  dx  lo  Jo*  '  dx' 


The  other  shear  term  gives  a  similar  result 


(B.ll) 

(B.12) 

(B.13) 

(B.14) 


(B.15) 


% 


where 


-i: 


-GA  — fiw, -  ““Ml  &<5  +  +  ^4  ~  ^1  -  <^4  “SMs  (B.  16) 


^4 

+d.u.6 - dM.6 —  +  d  t, 

dx  *  ^  ax 


ax 


,  aw 

f* — o — 

*  ax  ax) 


dx 


-Vo^^u^d^dx  -  -G4ju,in^  +  [^G^d^  ^6w  dx 
Jo  ax  j  0  ax 

^G^dju^^dx  -  G^d^u^i^-jjj^d^^dw  dx 

f^^.,  aw .  aw  j  ^.5  j  -K 

-  G^d^\-—d-^dx~ -G^d/t^-—dw\  +  G^d^\—6wdx 

Jo  ax  ax  lo  J 0  <?jc 


(B.17) 


The  kinetic  energy  term  is  dealt  with  by  switching  the  temporal  and  spatial  integrals, 
performing  the  integration  by  parts  on  the  variation  of  the  temporal  derivative,  and 
prescribing  the  variation  on  w  at  times  ti  and  ^  to  be  zero.  There  results 


dx 


Tlf 


(B.18) 


{pf^)-^dwdx  dt 
0  dr 


Substituting  these  integrations  back  into  Hamilton's  principle  yields  an  integral  of  the  form 

(EOM  for  w  )dw  +  (EOM  for  u,  )da,  +  (EOM  for  u,  )du,  +  (EOM  for  u,  Ida, )  dx 
+(the  expressions  that  make  up  the  boundary  conditions  )\dt  »  0  (B.19) 


The  individual  equations  of  motion  are  then 


'  .  a'w  au,  au.  \ 

ax  ax  ax ) 


^  j( j  d'w  au.  au.]  ^  ^a'w 


(B.20) 
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-ri- - ~r^  +  -:^“s - ga -rr  +  *  ■  .“1  +  ^<4 't: " ® 


<?a:‘  ^4 


dx  t. 


_  <?“i^  ^  G^  r-  j  n 

£^3  — ^  — ^113  +  ~^u,  -  G^d^  —~  “  0 
dx^  dx 

_  d'^u.  G.  G.  ^.dw 

-- ;^“5+G4d4— -0 

ax  t^  ax 


(B.22) 


(B.23) 


with  corresponding  boundary  conditions 


(-'’S-451-'' 

D^-G:dJ^d/,^*u,-u^-G^^dJ,^  +  u,-u^i-V  6w  -0 

(-£,<,^  +  F,|ai.,|  -0 


+  F,  &i,  -0 


(B.24) 


We  can  condense  this  set  of  four  equations  to  three,  but  before  we  can  do  this,  we 
multiply  (B.22)  by  E\t\IE5f5,  (B.23)  by  Ei/i/Ea/s,  and  (B.20)  by  Vf).  ui"  stays  the  same. 

♦  ii±i)±L +^f,. + 

dx*  Dt^-\  2  }  dx^  Dt^\  2  I  dx  Dt^\  2  }  dx 

DU-\-  2  I  -  dx-  DU\-  2  }  dx  DuS'  2  )  dx  D  dt^ 


d\  G.  G4  r'A  ^  ^  r- A  ^  f\ 

— ^ ^u^  +  —^u,  -  G^d^ +  - ^“1  +  ® 


dx~  t 


dx  r. 


(B.25) 

(B.26) 


„  d^u^  E,ttG^  EitfGn  Ej,G^  (  t,+t-i\dw  ^ 

El — f-t-  » I-. -u  -  — u_2.  ;,  +  J — 2.  — _o 

dx-  E3V,  E^y^  E/jf,  \-  2  }  dx 


(B.27) 
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(B.28) 


Now  subtract  (B.26)  from  (B.28)  and  subtract  (B27)  from  (B26)  to  obtain 


(B31) 


The  mode  shape  of  the  cantilevered  composite  beam  is  a  combination  of  cos  and  cosh 
terms.  However,  as  the  mode  number  increases,  the  majority  of  the  beam  has  the  mode 
shape  of  a  sine  wave.  Therefore,  we  assume  a  bending  mode  shape  of  sin  and  will 
accept  the  error  in  the  low  modes  where  a  larger  percentage  of  the  beam  is  affected  by  the 
cosh  part  of  the  displacement  expression.  In  addition,  the  axial  displacements  of  the 
constraining  layers  are  assumed  to  have  a  mode  shape  of  cos  "**4,  Rnally,  assume  a 
temporal  dependence  of  e'"*'.  The  assumed  solutions  have  the  form 

*v  -  VVe'“*sin-^^  and  u- ^  U je"’"  (B32) 


Operating  on  eqs.(B.29-3 1),  the  new  equations  become 
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(B33) 


(B35) 


At  this  point,  we  define  quantities  that  help  amplify  the  algebraic  expressions.  Define 

G. 


(f) 

(^4  +  ^ 

IIl 
2  ; 

1  ^3 

(t) 

+  ^1- 

IIl 
2  > 

I 

i+-5i 

^7 

^3 

<h 

£;r,r4  \rm) 

.  _S_f_Lr 


D  Vwr; 


\n3i}  D 


(B36) 


to  help  simplify  the  arithmetic.  Making  these  substitutions  and  putting  the  amplitudes  of 
the  expressions  in  matrix  form. 


l+^3‘?5 

-^6 

(t/.-t',) 

‘72<?6‘?7 

IV 

■^3 

1+^4<76. 

(B37) 
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Setting  the  detenninant  of  this  matrix  equal  to  zero  and  solving  for  the  frequency  oP  yields 
the  frequency  equation  for  free  vibration  of  the  system.  Since  only  qg  contains  nP  (with  a 
constant  coefficient),  we  need  to  isolate  qg.  Since  the  design  for  the  improved  treatment  has 
been  prescribed  to  be  symmetric  and  the  structure  and  constraining  layers  will  have  the 
same  modulus,  £i  =  £3  =  £5  =  £ ,  (5  =  /3  ,  =  /2 » and  G4*  =  G2*  =  G*.  Therefore,  qi  - 

q2,q3=<i4,  and  qs  =  q^.  Consequently, 

1  +  ^345  -^5 

det  q,q^^  1+ “0  (B38) 

“^5  1+<?3^5. 


Solving  the  determinant  forqg. 


Letting 


and 


(B39) 


we  now  endeavor  to  find  the  system  loss  factor  ijiy.  To  do  this,  we: 
'  find  the  positive  imaginary  square  root  of  qg 
-  express  [qg]*^-  as  £  +  //and  assume  the  response 

c'*sin/8f 


-  compare  this  with  the  known  response  of  a  damped  linear  system 

e"^"'"*sin  (Ojt 

-  consequently, 

L  ^ 

R  (o,  ' 
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Assuming  the  loss  factor  ri  =  2C, 


The  following  procedure  was  used  to  calculate  I  and  R: 
Since  G*  is  a  complex  quantity,  g*  is  too.  Let 

g’  -  Re(g* ) + Im(g*)  -  g,  +  g, , 

,  and  r-y. 


qs  is  rewritten  as 


<78 


l+g,X^g,X+g^XY-^&XY 
l+SrX  +  giX 


Rationalizing  the  denominator  and  letting  g/ = irjvgn 


r}.8rXY 


(l+&X)^  +  (»7,g^)^ 


where  rjv  is  the  viscoelastic  material  loss  factor. 

We  then  find  the  magnitude  and  angle  of  the  complex  square  root  of  qs’ 

mag(^8)  -  [Re(^8)’  +  Im(%)']^ 
angle (q,)  -  tan"' 

Re(q,) 

niag[(qg)^I  -  ^magCq,) 

Finally, 
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R  -  Re[(<7g)^  I  -  cos|angle((<7g)^lJ 

/  -Im[(^g)^]-  mag[(<?g)^jsiii|angle((%)^jj 

This  process  allowed  for  efficient  parametric  study  on  a  computer  spreadsheet 
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Undamped  Loss  Factor  eta 


Appendix  C  -  Comparison  of  Experimental  and  Theoretical  Modal  Loss  Factors  (Modes  1 
and  2)  for  Various  Constrained  Layer  Segment  Lengths  on  a  6061-T6  Al  IT  Beam 


Three  11"  beams  with  cross-section  0.126*  x  0.7S”  were  fabricated  for  the  purposes 
of  verifying  the  accuracy  of  the  3M  data  and  determining  the  repeatability  of  accurate  data. 
Ail  damping  segments  were  prefabricated  by  3M.  The  material  properties  are  described 
below  forT=70T.* 


Material 

Density  p 

E 

G 

V 

Corstraining  Layer 

.1  lb/in’» 

10  Mpsi 

3.8  Mpsi 

6.03 

Damping  Layer,  Mode  1 

.04  Ib/in^ 

240  psi 

80  psi 

0.49 

0.9 

Damping  Layer,  Mode  2 

.04  Ib/in3 

500  psi 

170  psi 

0.49 

0.9 

Table  C.1  3M  ISO  112  Material  Property  Data 


♦The  values  of  G*  Jf,  and  t]  for  the  damping  layer  are  courtesy  of  3M.  All  other  values 
were  transcribed  from  Mantena,  Gibson,  and  Hwang's  paper  [18]. 

Figure  Cl  shows  the  undamped  beam  data.  Each  data  point  consists  of  the  curve 
fitted  value  of  percent  damping  (multiplied  by  two  to  get  i\)  based  on  a  set  of  twenty  impact 
trials  whose  frequency  response  functions  were  averaged  into  a  single  FRF. 


^  Trial  (Each  trial  is  one  set  of  iwanty  averaged  Impacts 

Figure  C.  1  Undamped  Experimental  Loss  Factors  for  1 1 "  Beam  Modes  1  and  2 


The  corresponding  data  is  tabulated  below. 


Trial 

U2 

1 

0.0094 

6.02 

2 

0.0068 

0.016 

3 

0.0083 

0.0207 

4 

0.007 

0.018 

5 

0.0075 

0.0175 

6 

0.0076 

0.0152 

Table  C.2  Undamped  Modal  Loss  Factors 
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There  was  a  significant  amount  of  60  Hz  noise  in  the  signal.  This  effect  was 
minimized  by  using  the  maximum  number  of  samples  (4096)  for  each  FRF  and  by 
choosing  the  minimum  frequency  bandwidth  needed  for  testing.  In  addition,  the  accuracy 
of  the  data  was  sensitive  to  the  following: 

-  location  of  the  accelerometer  (cleaner,  more  consistent  data  was  obtained  when  the 
accelerometer  was  near  the  clamped  end  of  the  truss) 

-  impact  location  (clean,  gently  sloping  force-time  histories  gave  more  consistent  data 
than  sharp-impulse  and  double-hit  impacts) 

-  machine  warmup 

-  battery  strength  in  the  force  transducer  power  supply 

-  bandwidth  chosen  for  the  curve  fit  (more  samples/trial  and  wider  curve  fit 
bandwidths  gave  more  consistent  damping  percentages) 

-  amount  of  tension  (or  slack)  in  the  accelerometer  cable.  This  particular  variable  had 
the  greatest  influence  on  accuracy  and  repeatalnlity.  It  was  determined  through  trial-and- 
error  that  a  loose  cable  of  medium  length  (20”)  with  the  accelerometer  placed  as  close  to  the 
cantilevered  end  as  possible  was  the  most  desirable  setup.  Fortunately,  the  effect  of  the 
cable  tension  is  small  when  collecting  data  on  the  truss  because  the  truss  is  much  more 
massive  than  the  1 1"  beams. 

The  outcome  of  the  undamped  beam  testing  was  the  undamped  modal  loss  factors  for 
modes  1  and  2  (0.008  and  0.018,  respectively).  These  values  were  then  subtracted  from  the 
experimental  damped  modal  loss  factors  to  obtain  the  actual  damped  moal  loss  factors. 

The  following  data  represents  the  damped  modal  loss  factors  for  given  segment 
lengths  of  damping  treatment.  The  treatment  was  applied  over  the  entire  length  on  both 
sides  of  the  beam.  Each  data  point  is  double  the  curve  fitted  value  of  percent  damping  based 
on  the  averaged  FRF  of  ten  impact  trials.  Remember,  the  value  of  the  reduced  rj  =  exper.  rf 
-  undamped  rj. 
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Loss  Factor  eta 


Length 

reduced 

elal 

reduced 

eta2 

0.37 

0.0076 

0.011 

0.37 

00073 

00125 

0.75 

0.0131 

00228 

0.75 

00136 

00212 

0.75 

0.0136 

0027 

0.75 

0014 

00248 

0.75 

00148 

00208 

0.75 

0.0149 

0.0192 

0.96 

002 

00302 

0.95 

002 

0.0268 

0.95 

0.0204 

0031 

1.2 

0.0212 

00302 

1.2 

00224 

00278 

1.2 

0022 

0031 

1.56 

0.0304 

0.0274 

1.56 

0.0304 

0.0322 

1.56 

0.0298 

0.0396 

1.97 

00362 

0.0434 

1.97 

0.0362 

0.0418 

1.97 

0.0368 

0.0462 

Length 

reduced 

etal 

reduced 

eta2 

25 

25 

0.0408 

0038 

25 

00404 

00384 

306 

0.0408 

00384 

3.06 

00406 

00306 

306 

0.0404 

0033 

4 

00328 

00262 

4 

00332 

0023 

4 

O032B 

0.0234 

5 

0.031 

00296 

5 

00356 

00314 

5 

0.0366 

00306 

813 

0.0306 

0.0296 

813 

003 

0028 

813 

00302 

00266 

8 

0.0276 

0025 

8 

00274 

00238 

8 

00278 

00204 

10 

00258 

00172 

10 

0.027 

00134 

10 

0.0256 

00111 

Table  C.3  Experimental  Data,  11”  Beam  Experiment 


Figure  (43)  is  displayed  again  for  convenience. 


Segment  Length  (inches) 

Rgure  4.3  Experimental  vs.  Theoretlcai  r}  for  Various  Segment  Lengths 
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Appendix  D  -  Derivation  of  the  Differendal  Equation  of  Motion  ofTransverse 
Displacement  of  a  Three-Layer  Sandwich  Beam 


The  following  development  is  mostly  taken  from  Mead  and  Marinis'  article  [8].  S<mie 
steps  that  were  omitted  in  the  article  are  developed  in  full  here. 

Consider  again  the  slice  dbc  of  a  three  layer  sandwich  beam  of  unit  width. 


dx 


Figure  D.1.  Differentiad  Slice  dx  of  a  TTvee- 
Layer  Sarxjwich  Beam 


The  following  assumptions  are  made  for  this  analysis: 

-  the  beam  is  syimnetric  (£i  =  £3  =  £ ,  /a  =  ri) 

-  shear  strains  in  the  constraining  layers  are  negligible 

-  direct  strain  in  the  viscoelastic  layer  is  negligible 

-  the  viscoelastic  layer  thickness  t'l  is  constant  for  all  x  and  time 

-  the  transverse  displacements  of  each  layer  are  equal  at  any  location  x 

-  the  effects  of  shear  deformation  and  rotary  inertia  are  negligible 

Unlike  the  direct  shear  case,  the  shear  strain  in  the  viscoelastic  layer  has  two  terms  and 
is  given  by 


du  dw 
dz  dx 


(D.l) 


It  is  seen  that 

du  t^  ^-t^  dw  Uy-Uj 
dz  2/,  dx 
r,  dw  ,  t<3-u, 
r,  dx  fj 


(D.2) 


where  u/  are  the  longitudinal  displacements  of  the  constraining  layers.  The  shear  strain  is 
then 
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/  <2*  ^2 


(D3) 


where  y  is  the  average  shear  strain  through  the  thickness  f2'  The  shear  stress  is  equal  to 


(D.4) 


The  flexural  rigidity  D  of  the  face  plates  is  given  by 


■  12(1- v^) 


D-J,D, 


(DJ) 


where  i  are  the  elastic  layers  1  and  3. 

Mead  and  Marinis  neglected  poisson  effects,  but  the  addition  of  the  ( 1-v^)  term  adds 
little  complexity.  Since  the  thicknesses  of  the  layers  are  small  compared  to  their  width  b, 
this  term  will  remain. 

The  total  shear  force  on  a  cross-section  is  the  summation  of  the  shear  forces  acting 
through  the  thickness  of  each  layer.  Hence,  the  shear  force  S  is 


where 


S  “  5|  +  ^2  ^3 


««  V  wr  ,  _  V  wr 


(D.6) 


(D.7) 


The  quantity  S2  is  given  by  the  shear  stress  r  times  the  thickness  over  which  it  acts. 
Therefore, 

■S2--tJi+^|  (D.8) 


where  the  minus  sign  indicates  that  the  positive  shear  force  direction  is  the  opposite  as  the 
positive  shear  stress  direction.  Making  the  appropriate  substitutions. 


5.D^-G/2lA  + 

dr*  •‘I  1,  M  tjdx  U 


(D.9) 


108 


The  distributed  loading  internal  to  the  boundaries  is  the  derivative  of  the  total  shear  force. 


dx 


(D.IO) 


Since  the  sandwich  beam  will  have  no  externally  applied  forces  in  its  interior  except  for  the 
ineitial  loading, 

^  dx*  \  fj[\ 


where 


3 

(pr)  *  mass  of  the  sandwich  beam  per  unit  length  . 


Next,  consider  die  axial  forces  Pj.  Assume  these  forces  act  along  the  neutral  surfaces  of  the 
individual  constraining  layers  and  have  magnitucfes 


Pi  •  o,A,  •Et,- 


(D.12) 


It  is  noted  here  that,  for  a  unit  width,  Ai  =  U,  so 


and  Pj-Pf,-^ 
'  '  dr  ^  '  dx 


(D.13) 


Summing  forces  on  each  constraining  layer  gives 


dP,  ^  ^ 

— '---T  and  — '•-T 
dr  dx 


(D.14) 


Therefore, 


_  d^u.  2G .  -  J.  ^|dw 


(D.15) 


^  d'u,  2G,  ^  _1  ,  t.  \dw 


/,  I  dw 


(D.16) 
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Remember  that  we  are  considering  the  deformation  of  the  sandwich  beam  to  be  the 
superposition  of  the  direct  shear  quasi-static  case  and  the  total  bending  case  (quasi-static 
plus  dynamic).  Since  the  effects  of  the  externa)  axial  forces  have  already  been  considered, 
the  axial  load  will  be  assumed  to  be  zero.  Consequently,  for  no  axial  loading  and  a 
symmetric  layup. 


di^  di^ 
dx  dx 


(D.17) 


Subtracting  (D.16)  from  (D.15)  yields 


d\  2G  -GltAdw 


dx^  Et^tJ^  *  £f,  (  r,J  dtc 


(D.18) 


Rearranging  (D.  1 1 )  gives 


dx*  h)  dx^  u]  dx  '  dt^ 


(D.19) 


We  can  use  (D.18)  and  (D.19)  to  solve  for  w  by  taking  derivatives  to  eliminate  u\. 
Differentiate  (D.  19)  once  and  (D.  18)  twice  to  obtain 


d\  2G  du^  G 

dx^  Et^U  dx  fir,  I  f,  j  dx' 


d  w  Gt,  L  d  w  2G  t  d^iL  /  V  <?  w 

— r - -  l+—  — r - 1+—  — 3  +(p0 — 

dx  D  uj  dx  D  y  r,/  '  dx^dt' 


(D.20) 


(D.21) 


Substitution  of  (D.20)  into  (D.21)  gives 


GtJ,  A  V  2G(,  r,"!  2G  du,  G  (,  t,  .  d*w  ^ 

— r- - -  1+-^  — T - 1+-!-  - ^ - 1  +  -! - r  +(pr) — r-7-0 

dx  D\  /,)  dx  D\  Etft,  dx  fit,  ^  ^2}  dx'di 


(D.22) 


We  can  determine  My'  in  terms  of  w  by  using  (D.19)  as 
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(D.23) 


ih.  2G  d^w  ^  2G^  ^  d’w  2G{pt)  ^w 

dx  Et^t^  dx*  dx^ 


Substitution  of  (D.23)  into  (D.22)  yields 


d^w 

9k 

dx^~ 

D 

Li)  (d^) 

/,j  dx*  Et^t^  dx*  dt^  D  dx^dt^ 


We  now  define  the  shear  and  geometric  parameters  to  be 


*  Er/j  2D 


which  simplify  (D.24)  to 


d^w 

dx^ 


ipt)(  d*w 
D  \dx^dt^ 


0 


(D.25) 


This  is  the  sixth-order,  homogeneous  differential  equation  of  motion  for  transverse 
displacement  of  a  three-layer  sandwich  beam. 
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Appendix  E  -  Development  of  the  Volume  Specific  Energy  Dissipation  Per  Cycle  D 
For  this  development,  we  define  the  cmnplex  quantity  r*  as 

T  -  G  y  (El) 

where 

G‘  -  G'  +  iG”  and  y’  -  (y '  +  ly  . 

The  volume  specific  energy  dissipation  per  cycle  D  is 

D^^’^'-xrdt  (E2) 

We  must  integrate  the  real  parts  of  each  term  over  one  cycle  cS  oscillaticm  (equal  to  Inlaid. 
Re{G*y*}Re{/<uyV“‘}  dlf 

■  Re|(G'  +  /G")(y '  +  /y  ")|Re|/<u(y'  +  iy"Xcostar  +  /sinm/)|  dt 

Re^G'y  'costur  +  /G'y '  sincu/  +  iCycostot  +  /‘G'y  "sintw/  +  iCy' cosax 
+  PG”Y'sinQJt  +  i~G"Y"  cosox  +  ^G"y"sin  ajif|Re|/<uy'  cosax  +  /^toy 'sin  ax 
+  r<wy "costor  +  /^a>y"sin<ur|  d( 

^J^’*“(G'y  'cosa?/  -  G'y"  sin  wt  -  G"y 'sinoft  -  G"y"  costor)(-o;y'  sin  tor  -  tuy  "cos  tor)  dt 

-  to[-G'y''  costo/sintor  +  G'y  'y"sin'  cor  +  G"y''sin-  ax  +  G"y'y"costorsin  ax 

-G'y'y"cos"  tor  +  G'y"*  costorsintor  +  G"y'y"  costorsin  ax  +  G"y"^  cos'torjdr 

-  J^"*”to[(G'y 'y"  +  G"y'^)sin*  tor  +  (-G'y'y"  +  G"y"^)cos'  art 

+  (-G'y'”  +G"y  y''  +  G'y"^  +  G"yV")costorsintorjtif 

“to(“)[G"y''  +  G"y"=] 

.;iG"(y'^  +  y"^) 

-  rtG'lyp  (E3) 

The  energy  dissipated  per  cycle  Ds  is 
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£>,-  S«G‘]rf‘n' 

Vatmi 


(E4) 
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Appendix  F  -  Determnaton  of  the  Shear  Strain  Due  to  Berating  yj, 

From  Appendix  D,  we  know  that 

(F.l) 

(F.2) 


Solving  for 


EtA3\ 

2g(  dx^  ^ 


“i  -“3 


(F3) 


Also  from  Appendix  D, 

£fi  ifi  g 

dx"'  dx  "  gEt,u(^*tJt^){dx*  ^  ax^~  D  at^ 

Differentiating  this  quantity  gives 

a\  a\  D  f  ^  Y  d^w  (pt)  a-w  aw 

dtc' “  ax^  "  gEt^t._(i+t,fu)[ax^  ^  ax^  D  ar  ax 


(F.4) 


(F.5) 


_  .  .  <?"«3  f 

Subtracting - —  from  — t, 

ax'  ax' 


Et,  2D  a^w  ipt)a-waw 

2G gEt^tJ^  1 + ^, A2 ) [ ^ ^  ^  at^  ax\ 


(F.6) 


Rewriting  the  displacement  w{x,t)  in  a  convenient  shorthand  form, 
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ll-l 


„  .  iwcc 
Csm - 


L  CO. 


cOy  sinco^j  f  -  CO sinco 


/J 


(F.7) 


+Afj  sincOyT^OiJC  + 


where 


-[(«i-.E|(g-.,i{i)'lg.|(i!C-.rl(i).(.!.i-«l(f 


<h 

<h 

a^m- 


l2L-gYe 

TAD 

_1_ 

ID 

gYL 

\2D 

ML 

TAD 


At  this  point,  we  are  ready  to  take  derivatives  and  find  an  expression  for  /£>. 


d  d‘w  d  ^ 

dx  dt^  dx 


C.sin- 


/  2  '  A 

wcc  Mj^ 

<»TS>n  tOjg  t  +  cOy-cOjB  sinctfyt 

t 

-  (Oy- Afj,  sin(07/(a,x  +  a,x^  +  a,x^  +  a^x*^ 


=  2 

n-l 


rm 


C.  cos 


ruvc  Aifc 


L  (o 


SB. 


-(Ogg  *COySin  (Ogg  t  +  (Or'^sg  S*** 


Wj?  -0)2 


/J 


(F.8) 


-iOfMg  sintOyr^a,  +  2a^  +  3a^x^  +  4a^^) 


d^w  ” 
dx^  «?i 


I  wry  /wt 

-It) 


rmx  M. 


<0. 


cOy  stnco5g^r  -  oi^g^  sincoy/ 


sn.  \  ~  ^SbI 


/J 


(F.9) 


+ AfjSin  €i>jt{6aj  +  24a4x) 
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77-.?[lTj  ^;ti — — ^)| 


(F.IO) 


Substitution  of  these  terms  into  (F.6)  gives 


2  Ceos 


2G\gEtyU(\  +  tJt^) 


rutx 

ptrut  Mg 

L 

D  L  (Ogg 

-(Ogg^  (Oj. sintOgg  t  +  (Oj?  fOg^  sinaij.t 


tt>r  -^sb: 


KrutV  JnnV]  Mg  fa>r  sincUj^/- to ja,  sin o)^ 
l]  ^^LjJtu,^  o>j?^o>ggl 

+A4sint07./|-gl'(6t^  +24a^)  -^la^x-^ZayK^  +  4t24jr’)j| 


(F.ll) 
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Appendix  G  -  Determination  of  Improved  Treatment  Parameters  h,  s,  t,  and  the  Cross-Section 


In  designing  the  improved  treatment,  one  of  the  limiting  constraints  was  that  "off-the- 
shelf*  parts  must  be  used.  This  requirement  limits  the  number  of  choices  of  constraining 
layer  thickness  t  and  the  cross-section  to  the  following  choices: 

t:  Since  aluminum  was  chosen,  the  reasonable  choices  for  thickness  were  0.025”, 
0.032",  and  0.040".  Thinner  sheets  required  special  ordering,  and  thicknesses  greater  than 
0.040”  possess  unacceptably  high  bending  stiffnesses. 

cross-section:  The  AfTT  shop  had  several  thick  sheet  sizes  and  a  beam  whose  cross- 
section  measured  0.12S"  x  0.75".  The  sheet  si^s  considered  to  be  reasonable  were  0.062", 
0.080",  and  0.090".  Of  course,  the  other  dimension  could  be  arbitrarily  chosen,  so  there 
were  actually  an  infinite  number  of  choices.  Since  the  task  of  choosing  a  cross-section  was 
considered  peripheral,  the  field  of  choices  was  arlritrarily  narrowed  to  the  following: 

Cross-Section  1:  0.062"  x  0.75" 

Cross-Section  2:  0.080"  x  0.75" 

Cross-Section  3:  0.090"  x  0.75" 

Cross-Section  4:  0.090"  x  1.0" 

Cross-Section  5:  0.125"  x  0.75". 

Some  rough  calculations  showed  that  the  optimal  range  of  values  for  h  and  s  were  4"- 
8"  and  1 .0"- 1 .5",  respectively.  The  choices  for  h  and  s  were  chosen  to  be  the  following: 

h:  5",  6",  and  7" 
r.  1.0",  1.25",  and  1.5" 

The  total  number  of  configurations,  then,  was  135.  The  total  width  of  the  sandwich  beams 
was  backed  out  by  using  the  constraint  of  maximum  weight  of  (levers  +  sandwich  beams)  = 
1.0  lb.  However,  the  proportion  of  the  width  in  each  bay  is  another  variable,  which 
complicates  matters  even  further. 

The  general  approach  of  choosing  the  ‘^besf*  configuration  was  the  elimination  of 
one  variable  at  a  time.  Rnite  element  modeling  was  chosen  as  the  tool  for  the  process. 
Basically,  one  parameter  was  varied  in  a  nominal  configuration  while  all  other  variables 
(except  hi)  were  held  constant.  The  particular  widths  bi  and  h2  were  assigned  widths  of 
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O.Tb/oio/and  QSbtouu^  respectively.  The  key  parameter  in  the  elimination  of  each  variable 
was  the  ratio  of 


^  elastic  strain  energy  in  the  sandwich  beam  elements 
^elastic  strain  energy  in  all  other  elements  (levers  -i-  truss) 

To  fix  parameter  i,  several  configurations  were  tested  where  the  only  difference 
between  them  was  the  value  of  parameter  i.  The  configuration  which  had  the  highest  value 
of  r  was  chosen,  and  the  process  was  repeated  for  next  parameter.  Parameter  i  was  then 
held  constant  for  the  remainder  of  the  study. 

One  of  the  peculiarities  of  the  1-DEAS  software  is  its  inability  to  allow  the  user  to 
bring  two  beam  elements  together  to  a  pin-jointed  connection.  The  use  of  multi-point 
degree-of-freedom  constraints  is  requited,  which  involves  the  creation  of  nearly  coincident 
nodes  that  are  constrained  to  have  the  same  values  for  the  desired  DOF.  The  software 
continued  to  show  a  non-zero  moment  at  the  end  of  a  simply  suf^ited  beam  (intuitively  an 
impossibility).  The  result  of  this  discrepancy  was  a  linearly  varying  value  of  strain  energy 
in  the  eight  elements  that  comprise  the  sandwich  beam.  This  makes  no  sense,  so  a  moment 
of  equal  magnitude  and  opposite  sense  had  to  be  administered  to  enforce  the  condition  of 
zero  moment  at  the  lever  tip.  This  correction  resulted  in  a  uniform  distribution  of  strain 
energy  over  the  length  of  the  sandwich  beam.  Although  the  actual  distribution  of  strain 
energy  is  non-uniform  (governed  by  the  beam-column  effect  of  the  sec  phenomenon), 
the  expected  result  from  the  model  is  uniform  because  the  solver  uses  linear  strength  of 
materials  theory  (which  assumes  circular  curvature  for  equal  moments  applied  at  both  ends 
of  an  element). 

Some  configuration  choices,  upon  closer  review,  were  ruled  out  immediately  or  were 
eliminated  because  of  too  high  a  value  of  ^/2  (which  produces  unacceptable  asynunetry  of 
cycle  halves).  A  discussion  of  the  ^/2  phencmienon  is  in  Chapter  5. 

The  following  is  a  list  of  judgments  that  eliminated  possible  configurations.  For  f,  the 
choice  was  rather  simple,  t  =  0.025"  caused  excessive  deflections  {^h  was  too  high),  while 
0.040"  was  too  stiff,  not  allowing  adequate  bending  deflections.  Hence,  t  =  0.032"  was 
chosen. 

In  choosing  the  cross-section,  it  was  determined  that  4  and  5  were  excessively  beefy 
and  merely  reduced  the  sandwich  beam  widths  because  of  their  weight.  Cross-section  1 
was  eliminated  because  it  is  dubiously  flimsy  in  bending  out  of  its  own  plane.  A  combined 
torsion-bending  motion  might  easily  cause  lateral  instability  of  such  a  thin  cross-section. 
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Cross-sections  2  and  3  were  tested  in  I-DEAS,  and  number  2  resulted  in  a  higher  value  of 
the  ratio  r. 

Since  the  choices  of  h  and  s  are  so  heavily  coupled,  they  were  varied  simultaneously. 
Two  test  case  finite  element  models  were  run  to  investigate  the  extreme  values  of  /i  =  5"  and 
5=1".  The  solution  indicated  that  choosing  h  =  5"  does  not  offer  enough  lever  tip 
displacement,  while  choosing  5=1"  results  in  unacceptably  high  values  of  ^/2  for  h  equal 
to  6"  or  7".  Therefore,  the  following  five  models  were  chosen  as  candidate  configurations 
based  on  the  individual  advantages  of  each.  The  models  were  the  following: 

Model  1: /I  =  5", 5=1", 

Model  2: /i  =  6",  5=1", 

Model  3: /i=  6",  5=  1.25", 

Model  4:  A  =  7"  ,  5  =  1.25",  and 
Model  5: /i  =  7",  5=  1.5". 

Model  4  produced  the  highest  ratio  r  and  was  thus  declared  the  winner.  Finally,  a  simple 
test  was  run  between  three  models  (configured  like  model  4)  whose  only  difference  was  the 
distribution  of  the  total  sandwich  beam  width  btot  between  bays  1  and  2.  The  three 
distributions  were  (65%  in  bay  1  and  35%  in  bay  2),  ^%o.  and  ^^/25-  The  ^%o 
configuration  produced  the  highest  ratio  r,  although  this  particular  test  produced  very  small 
differences  in  the  strain  energy  distribution  in  the  system.  The  final  configuration  is  then 

long  lever  length  ^-7in. 
short  lever  length  5  »  1.25  in. 
sandwich  plate  thickness  »  0.032  in. 

cross  -  section  -0.75  in.  x  0.080  in. 
wide  sandwich  beam  width  -  54  in. 
slim  sandwich  beam  width  -  23  in. 
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Appendix  H  -  MATLAB  Junction  file  LOSS 

For  this  program,  the  input  P  represents  the  tip  load  given  by  the  finite  element  model.  The 
FE  model  was  loEided  like  the  figure  shown  at  the  beginning  of  chapter  5. 

function  [sse,energy]=loss(z) 

%  lX)SS(z)  calculates  the  energy  dissipated  per  cycle  in  a  sandwich  beam  due  to 
%  induced  shear  and  moment  Fora  1  X  12  vector  z  whose  inputs  are  [^G*,tl,t2,s,P, 
%  omegsb,  omegt  L,  rhot,  b,  4-digit  model#],  LOSS(z)  returns  the  value  of  the  definite 
%  inte^al  of  tot^  (shear  strain)^!  over  the  area  of  the  SB  and  the  corresponding  energy 
%  dissipated  per  cycle. 

E=z(l); 

G=z(2); 

tl=z(3); 

t2=z(4); 

s=z(5); 

P=z(6); 

omegsb=z(7); 

omegt=z(8); 

L=z(9); 

rhot=z(10); 

b=z(ll); 

g=2*G/(E*tl*t2) 

D=E*tl''3/(6*.8911) 

Y=(tl+t2r2*E*tl/(2*D) 

al=(  12*L-g*Y*L''3)/(24*D) 
a2=-l/(2*D) 
a3=g*Y*L/(12*D) 
a4=-g*Y/(24*D) 

vnsq=2/(L*((pi/L)''2+g)) 


cape  1  =-g*L''2/pi  *a  1 

capc2=-a2*(L/pi*(g*L''244)+4*g*(Upi)''3) 

capc3=-a3*((6*L-g*L''3)*L/pi+6*g*L*(Upir3) 

capc4=-a4*((  12*L^2-g*LM)*L/pi+(  12*g*L''2-48)*(L/pi)''3-48*g... 

♦(UpD^^S) 

capc5=(  1  +Y)*L/pi*2*g''2*  Y/D 

capc=vnsq*(omegt''2*(capc  1  +capc2+capc3-«-capc4)+capc5) 

gamtcol=D*(P*s/b)/(G*g*(tl+t2)) 

gamtco2=((pi/L)''5-g*Y*(pi/L)''3)/(omegt^2-omegsb''2) 

gamtco3=rhot*omegt^2/D 

gamsco=(P/b)/(sqrt(g)*E*t  1  *t2)/((-tanh(sqrt(g)*L/2)... 
*cosh(sqrt(g)*L)+sinh(sqrt(g)*L))) 

hiterm=gamtco  1  ♦g*Y*6*a3 
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dx=.2; 

ii=:l; 

x=.l; 

iiintgrl=30; 

intgri=0; 


while  n<65 

gamt(ii)=gamtco  1  ♦(capc*gamtco2*cos(pi*x/L)-g*Y  ♦(6*a3+24*a4*x). 
+gamtco3*(capc*pi/(L*(omegt'^2-omegsb''2))*cos(pi*x/L)... 
-(al+2*a2*x+3*a3*x'^2+4*a4*x'^3))); 
gains(n)=:gainsco*(cosh(sqrt(g)*x)-taiih(sqit(g)*L/2)... 

*sinh(sqrt(g)*x)); 

gamtotal(ii)^gaint(n>t-gains(ii); 

uenrgy(ii)=^*(feal(gamtotal(n)))^2: 

senrgy(n)=dx*(abs(gaintotal(n)))^2; 

uint^=uiiitgrl+ueiirgy(ii); 

iiitgri=iiitgrl+seiirgy(n); 

unintot(n]^uintgrl; 

runtotal(n)=int^; 

x=x+dx; 

n=n+l; 

end 


enrgyl=uenrgy'; 
enrgy2=senrgy'; 
enrgyniat={enrgyl  enigy2]; 
urun=uruntot'; 
ran=runtotal'; 
roninat=[unin  run]; 

energy=pi*imag(G)*t2*intgrl 

benergy=b*energy 

uenergy=.5*real(G)*t2*uintgrl 

buenergy=b*uenergy 

gamtd) 

gamt(32) 

gaint(64) 

r=.l:.2:12.7; 

cig; 

plot(r,gamt) 

grid 

tiUe( ['Bending  Shear  Strain  (rad),  t2  =  ',nuni2str(z(4)),'  in.'|) 


121 


Appendix  I  -  Experimented  Results  of  Damping  Tests  -  Raw  Data 
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2 

1R 

1 
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5.798 
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3 

1R 

1 
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4 

1R 
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7.472 
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14.5 
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5 

1R 
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14.56 
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6 

1R 
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7 

1R 
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7.563 
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8 

1R 
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1R 

6R 
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14 
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IR 
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Table  1.2  Damping  Tests,  Improved  Treatment,  t2  =  0.015" 
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